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ABSTRACT 


The  objective  of  the  STAR  project  is  to  test  and  space 
qualify  a  new  continuous  cycle  cryogenic  refrigeration 
system  for  cooling  of  sensors  and  electronics  which  is  based 
upon  the  newly  discovered  thermoacoustic  heat  pumping 
effect.  The  new  refrigerator  has  no  sliding  seals,  a  cycle 
frequency  of  about  300  H2,  and  uses  acoustic  resonance  to 
enhance  the  overall  power  density  and  efficiency.  This 
thesis  is  concerned  specifically  with  the  design  and  testing 
of  the  electrodynamic  transducer  which  is  responsible  for 
the  electro-acoustic  power  conversion.  A  computer  model  of 
the  driver /resonator  system  is  presented  along  with  the 
techniques  for  measurement  of  the  electrical  and  mechanical 
parameters  used  as  input  for  the  model. 

A  final  driver  design  (including  dimensional  drawings) 
utilizing  a  modified  JBL  2450J  neodymium-iron-boron 
compression  driver  and  associated  leak-tight  electrical 
feed-throughs,  microphone,  accelerometer,  pressure  gage, 
pressure  housing,  and  resonator  interface  is  provided. 
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I .  INTRODUCTION 


A .  BACKGROUND 
1 .  History 

Thermoacoustics  can  generally  be  described  as  the 
study  of  the  interaction  between  heat  and  sound.  Scientific 
interest  in  this  interaction  is  not  new.  Lord  Rayleigh  (1878 
and  1945)  discussed  various  qualitative  examples  of  the 
production  of  sound  by  heat.  In  one  of  these  examples  he 
describes  investigations  by  Sondhauss  (1850)  of  an  effect 
observed  by  glassblowers :  heating  a  bulb  of  glass  at  the 
end  of  a  hollow  rod  occasionally  led  to  the  generation  of 
sound.  Faraday  (1818)  demonstrated  that  this  effect  occurred 
with  gases  other  than  air.  Sondhauss  proved  that  the  vibra¬ 
tion  of  the  glass  itself  did  not  generate  the  sound,  but  he 
offered  r.o  explanation  as  to  what  did.  In  his  description  of 
Sonhauss'  work,  Rayleigh  stressed  the  importance  of  the 
phase  difference  between  temperature  and  particle  motion. 
Rijke  (1859)  built  an  open  tube  with  a  wire  mesh  inside. 
When  the  mesh  was  heated ,  the  tube  produced  a  sound  of 
impressive  intensity.  The  functioning  of  this  oscillator  is 
slightly  different  since  it  requires  steady  gas  flow  for  its 
operation.  A  similar  sound  production  effect  was  observed  by 
Taconis ,  et  ai.  (1949)  in  hollow  tubes  immersed  in  liquid 
helium  (4.2«j^)_  The  Taconis  oscillations  were  addressed 
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quantitatively  by  Yazaki,  et  al.  (1980).  The  work  described 
above  deals  primarily  with  prime  movers — devices  that  con¬ 
vert  a  temperature  gradient  to  sound  energy.  We  call  this 
the  classical  era  of  thermoacoustics. 

Interest  in  thermoacoustics  was  renewed  when  the 
idea  of  the  use  of  acoustical  devices  as  refrigerators 
appeared.  Gifford  and  Longsworth  (1966)  described  the  pump¬ 
ing  of  heat  along  a  surface  caused  by  a  periodic  change  in 
the  pressure  of  the  adjoining  gas.  Such  a  change  can  be 
produced  by  the  oscillations  of  a  sound  wave.  In  their 
experiments  Merkli  and  Thomann  (1974)  explored  the  heating 
and  cooling  effects  on  the  wall  of  a  gas-filled  resonant 
tube.  They  found  that  heat  was  transported  from  a  region 
near  the  velocity  antinode  (or  maximum  value)  of  the  sound’s 
standing  wave  to  regions  near  the  adjacent  pressure  anti¬ 
nodes.  Thermoacoustic  prime  movers  have  also  received  recent 
attention.  Kempton  (1976)  investigated  the  excess  noise  of 
aeroengines  above  that  predicted  by  theory.  He  determined 
that  the  sound  was  produced  by  unsteady  heat  transfer.  Each 
of  these  latter  three  groups  of  exprimenters  used  some 
theory  for  comparison,  but  it  was  mostly  qualitative.  None 
provided  the  complete  quantitative  theory  that  would  explain 
their  experimental  observations . 

The  theoretical  breakthrough  in  thermoacoustics 
field  was  made  by  Nicklaus  Rott  (1969,  1974,  1975,  and 
1980).  His  theory  combined  basic  principles  from  physics. 


thermodynamics,  and  acoustics  to  quantitatively  describe  the 
effects  found  in  both  types  of  thermoacoustic  devices:  prime 
movers  and  heat  pumps.  Prime  movers  use  a  temperature  gra¬ 
dient  to  create  sound,  as  in  the  effects  discussed  by  Sond- 
hauss,  Rayleigh,  Taconis,  and  Kempton.  Keat  pumps,  or  refri¬ 
gerators,  use  the  oscillating  pressure  of  a  sound  wave  to 
produce  a  temperature  difference,  as  described  by  Gifford 
and  Longsworth  ana  Merkli  and  Thomann.  Rott  described  the 
effect  found  by  Merkli  and  Thomann  as  thermoacoustic  stream¬ 
ing.  It  is  this  effect  that  makes  thermoacoustic  refrigera¬ 
tion  possible. 

Inspired  by  Ceper ley's  (1979)  traveling  wave  Ster¬ 
ling  Cycle  heat  engine  and  Rott's  quantitative  theory,  the 
team  of  Wheatley,  Hofler,  Swift,  Migliori ,  and  Garrett 
(1982,  1983a,  1983b,  19e5,  and  1986)  developed  a  series  of 
thermoacoustic  experiments  at  Los  Alamos  National  Laboratory 
in  New  Mexico.  They  investigated  the  basic  thermoacous ”.ic 
effects  in  both  prime  movers  and  refrigerators  and  compared 
their  experimental  results  to  Rott's  theory. 

2 .  Thermodynamics 

We'll  digress  from  history  here  to  explain  the  ther¬ 
modynamic  distinction  between  prime  movers  and  refrigera¬ 
tors.  This  discussion  follows  Sears  and  Salinger  (1975).  A 
prime  mover  receives  heat  at  a  high  temperature,  TM  (the  hot 
reservoir),  does  work  on  its  surroundings,  and  rejects  heat 
at  a  lower  temperature,  Tc  (the  cold  reservoir) ,  as  shown  in 
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Figure  i-l.a.  The  first  law  of  thermodynamics  tells  us  that 
the  heat  flow  (Oh)  from  the  hot  reservoir  must  equal  the 
work  done  on  the  surroundings  (W)  plus  the  heat  flow  (Qe)  to 
the  cold  reservoir: 

Qh  =  W  +  Qc.  (1-1) 

The  second  law  of  thermodynamics  states  that  in  every  pro¬ 
cess  the  entropy  of  the  universe  must  either  remain  constant 
or  increase.  The  entropy  change  for  the  isothermal  heat 
transport  to/from  the  reservoirs  is  defined  as  the  heat  flow 
divded  by  the  reservoir  temperature.  For  a  prime  mover  the 
second  law  can  be  written: 

ASm  *  A Sc  or  Qm/Tm  i  Qc/Tc.  (1-2) 

The  thermal  efficiency  of  the  prime  mover  is  defined  as  the 

ratio  of  the  work  output  to  the  heat  input: 

n  *  W/Qm  *  ( 0*4  -  Qcs)/Qm  =  1  -  Qc/Qm.  (1-3) 

Comparing  equations  1-2  and  1-3  gives  us  a  limit  on  the 

efficiency  based  on  the  temperatures  of  the  reservoirs: 

nil-  tc/tw.  (1-4) 

This  limit  is  known  as  the  Carnot  efficiency. 

In  a  refrigerator  the  process  is  basically  reversed 
(see  Figure  I-l.b).  Heat  (Qc)  is  removed  from  the  low  tem¬ 
perature  reservoir,  work  (W)  is  received  from  the  surroun¬ 
dings,  and  heat  (Qw)  is  delivered  to  the  higher  temperature 
reservoir.  The  expression  for  the  first  law  is  the  same  as 
for  the  prime  mover,  Equation  1-1.  For  the  refrigerator  it 
means  that  the  heat  flow  into  the  hot  reservoir  must  equal 
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Figure  1-1.  The  Thermodynamics  of  (a)  Prime  Movers 
and  (b)  Refrigerators. 
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the  heat  flow  from  the  cold  reservoir  plus  the  work  done  to 
the  systen.  The  second  law  for  this  case  is: 

ASm  i  &So  or  Qm/Tm  A  Oe/To .  (1-5) 

The  efficiency  parameter  for  the  refrigerator  is  called  the 

coefficient  of  performance  (COP) .  It  is  defined  as  the  ratio 

of  the  heat  flow  from  the  cold  reservoir  to  the  work  input 

to  the  refrigerator,  or: 

COP  «  Qo/W  -  Qe/<Qn-Qo>  S  Tc'(Tm-To).  (1-6) 

The  limit  is  known  as  the  Carnot  coefficient  of  performance. 

The  work  in  this  thesis  deals  solely  with  refrigera¬ 
tors.  Even  though  prime  movers  are  mentioned  from  time  to 
time,  the  primary  focus  of  the  rest  of  our  discussion  will 
be  refrigerators. 

3.  Qygr.aU  BHislansz 

According  to  Hofler  (1986),  the  early  thermoacoustic 
work  at  Los  Alamos  focused  on  experimental  refrigerating 
engines.  The  performance  of  these  engines  fell  short  of 
expectations,  leading  to  simple  experiments  on  basic  thermo¬ 
acoustic  effects  and  a  proof -of -principle  experiment  on 
thermoacoustic  refrigeration.  Hofler  then  applied  the  Rott 
theory  to  the  experimental  systems  and  solved  the  resulting 
equations  numerically.  For  his  doctoral  dissertation  from 
the  University  of  California,  San  Diego,  Hofler  designed  and 
constructed  a  completely  functional  thermoacoustic  refriger¬ 
ator.  He  also  made  accurate  measurements  of  its  thermody¬ 
namic  efficiency,  and  used  this  efficiency  to  make  compari- 
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sons  to  the  Rott  theory.  After  receiving  his  doctorate, 
Hofler  case  to  the  Naval  Postgraduate  School  (NPS)  as  a 
post>doctoral  fellow,  and  brought  his  prototype  refrigerator 
with  his . 

The  purpose  of  this  thesis,  in  conjunction  with  the 
work  of  several  other  students,  is  to  modify  Hofler's  re¬ 
frigerator  design  in  order  to  improve  its  overall  efficien¬ 
cy  and  make  it  suitable  for  space  cryocooler  applications. 
The  basic  design  of  the  thermoacoustic  refrigerator  is  shown 
schematically  in  Figure  1-2.  The  driver  (A),  which  produces 
the  sound  waves,  is  coupled  via  bolts  to  the  resonator  (D) 
via  a  reducer  cone  (B)  and  oellows  (C>  Inside  the  resonator 
is  a  stack  of  plastic  plates  (E)  and  their  heat  exchangers 
(F)  which  allow  the  heat  to  be  removed  from  the  hot  end  and 
absorbed  by  the  cold  end.  It  is  the  interaction  between  the 
sound  waves  and  the  plastic  plates  that  produces  a  tempera¬ 
ture  difference  across  the  plates  and/or  pumps  heat  from  the 
cold  heat  exchanger  to  the  hot  heat  exchanger.  A  brief 
description  of  efficiency  would  be  useful  here.  Figure  1-3 
shows  a  diagram  of  the  chain  of  energy  transformations  which 
occurs  in  a  thermoacoustic  refrigerator.  There  are  three 
conjugate  variable  pairs  (V  and  I;  P  and  U;  and  AT  and  <3)c)  , 
and  two  transformations  (electrical-to-acoustical  and  acous- 
tical-to-thermal ) . 

We  start  with  electrical  input  power  to  the  driver, 
which  can  be  calculated  (using  r.m.s.  values)  with: 
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Figure  1-2.  Basic  Theraoacoustic  Refrigerator  Showing 

(A)  Driver,  (B)  Reducer  Cone,  (C)  Bellows, 
(D)  Resonator,  (E)  Stack,  and  (F)  Heat 
Exchangers . 
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Figure  1-3.  Chain  of  Energy  Transformations  in  a 
Thermoacoustic  Refrigerator. 
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whore  I  is  tho  current  into  the  electroecountic  driver,  V  it 
the  voltage  across  the  driver  terminals,  e  is  the  phase 
angle  between  I  and  V,  and  cos  e  is  known  as  the  "power 
factor."  An  electrodynaaic  driver  converts  the  electrical 
power  to  acoustic  power  (the  first  transf orsation ) .  Analo¬ 
gous  with  the  previous  definition  of  the  electrical  input 

power,  the  acoustic  power  is  given  by: 

fi  •  P  U  cos  #,  (1-8) 

where  P  is  the  acoustic  pressure,  U  is  the  voluaetric  velo¬ 
city  which  equals  the  product  of  particle  velocity  and 
resonator  cross-sectional  area,  and  #  is  the  phase  angle 
between  thes.  The  efficiency  of  this  electroacoustic  trans- 
f orsation  is  given  by: 

n«*  -  ft/Pg (1-9) 
In  the  second  transf orsation,  the  refrigerator  con¬ 
verts  the  acoustic  power  to  a  temperature  gradient  (vT)  and 
heat  flow  (G) .  As  discussed  previously,  the  efficiency  for 
this  transition  is  given  as  the  coefficient  of  performance: 

COP  -  0c/».  (I -10) 

The  overall  theraoacoustic  efficiency  is  therefore  the  pro¬ 
duct  of  tho  efficiencies  for  the  two  transformations: 

=  Hba  *  COP  *  ^c/Pbukc ■  (I— 11) 

Improvements  in  the  overall,  refrigerator  efficiency 

can  thus  te  divided  into  two  distinct,  though  interrelated. 

efforts.  The  first  is  the  suoject  of  this  thesis:  modifica- 
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tions  to  the  driver  to  optimize  the  electroacoustic 
efficiency.  Tha  second  affort  involves  tha  optimization  of 
tha  coafficlant  of  parforaanca  and  ia  addraaaad  by  LT 
Michaal  P.  Susalla,  USN,  in  hia  master's  thaaia. 

B.  SPACE  THERMOACOUSTIC  REFRIGERATOR  (STAR) 

1.  Motivation 

In  addition  to  improving  ovarall  efficiency,  our 
axpariaantal  refrigerator  will  be  designed  and  built  with 
tha  intantion  of  flying  it  on  tha  Spaca  Shuttla  as  part  of 
tha  National  Aaronautica  and  Spaca  Administration's  (NASA) 
Get  Away  Spacial  (GAS)  program.  As  of  3  February  1988,  thara 
is  a  Memorandum  of  Agreement  between  tha  Naval  Postgraduate 
School  (NPS)  and  tha  Air  Force  which  provides  funding  for 
tha  space  flight  and  which  assigns  the  NASA  payload,  number 
G-337  to  this  project.  Tha  reason  for  testing  the  thermo¬ 
acoustic  refrigerator  in  apace  is  that  the  immediate  appli¬ 
cations  of  the  STAR  are  expected  to  be  the  cooling  of  elec¬ 
tronics,  high  Tc  superconductors,  and  infrared  detectors  in 
space. 

There  are  currently  two  primary  methods  for  cooling 
satellite  systems  in  space:  (1)  evaporation  of  expendable 
cryogens  (liquid  helium,  nitrogen,  ammonia -methane,  and 
solid  hydrogen,  etc.),  and  (2)  closed  cycle  refrigerators 
(Stirling  cycle,  Vuilleumier  cycle,  etc.),  which  involve 
large  reciprocating  piston  masses  operating  at  low 
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frequency.  [Walker  (1903)  and  Salth,  at  al.  (1984)]  The 
disadvantages  of  thaaa  currant  cooling  aathods  ara  thair 
short  lifatiaaa  (axpandabla  cryogena)  and  high  vibration 
lavala  and  low  raliabllity  (cloaad  cycla  raf rigarators ) .  Tha 
advantagas  of  tha  tharaoacoustic  rafrlgarator  ovar  thaaa  and 
othar  ayataaa  liaa  in  its  aiaplicity  and  raliabllity.  Tha 
STAR  has  no  sliding  seals,  afficiant  performance,  low  vibra¬ 
tion  lavala,  and  (hopefully)  a  longer  Ufa  span. 

Tha  tharaoacouatic  rafrlgarator  naada  to  be  surroun¬ 
ded  by  a  vacuus  and  insulation  material  to  minimize  heat 
loss  due  to  tharsal  conduction  and  radiation.  The  vacuus  of 
space  will  provide  the  bast  insulation  to  determine  the 
refrigerator's  ultimate  efficiency.  Tha  absence  of  gravity 
removes  the  possibility  of  thermal  convection  as  an  addi¬ 
tional  nuisance  heat  transport  mechanism. 

2 .  Qet  Away  Special  (GAS)  Proaram 

NASA's  GAS  program  allows  small,  self-contained 
payloads  to  fly  on  the  Space  Shuttle  in  GAS  canisters  for 
relatively  low  cost  ($10,000)  (Get  Away  Special  Team,  1984). 
The  GAS  can  is  five  cubic  feet  in  volume  and  can  house  a 
payload  of  up  to  200  pounds .  Each  payload  must  contain  its 
own  electrical  power,  control,  data  acquisition  and  storage 
facilities,  etc.  The  Shuttle’s  astronauts  will  turn  a  switch 
on  or  off  at  designated  times  during  the  flight,  but  are 
otherwise  not  involved  with  the  experiment. 


3. 


Shand  Subsystems  fret  NASA  3-313 

Another  group  of  UPS  students  end  fsculty  (Boyd, 
et  si.,  19S7)  hsve  taken  advantage  of  the  QAS  prograa  to 
eeasure  the  resonant  acoustic  nodes  of  the  shuttle  payload 
bay  and  the  aabient  acoustic  environsent  produced  as  a 
result  of  sain  engine  and  booster  operation  during  launch. 
This  experiment  is  titled  "The  Space  Shuttlu  C'rgo  Bay 
Vibroacoustics  Experiment”  and  is  designated  by  NASA  as 
payload  G-313.  Several  subsystess  that  were  developed  for 
NASA  0-313  will  be  used  by  NASA  0-337  (STAR).  A  schenatic  of 
the  STAR  in  its  GAS  can  is  shown  in  Figure  1-4.  one  of  the 
systems  borrowed  from  NASA  G-313  is  the  electronics  system 
used  to  run  the  experiment  and  record  the  data  The  recorder 
system  consists  of  the  INTEL  model  BPK  5V7S  magnetic  bubble 
memory  module  and  an  NSC  800  microprocessor-based  control¬ 
ler.  Two  other  NPS  students,  LT  Charles  B.  Cameron,  USN  and 
CPT  Ronald  Byrnes,  USA,  will  be  designing  the  analog  elec¬ 
tronics  and  software  to  integrate  these  systess  into  the 
STAR  experiment  for  their  master's  thesis.  A  block  diagram 
of  the  electronics  is  shown  in  Figure  1-5. 

Another  borrowed  system  is  the  power  supply,  which 
consists  of  Gates  brand  lead  acid  battery  cells  (five 
ampere-hour,  two  volts  each).  These  gelled  electrolyte  bat¬ 
teries  are  ideal  for  the  STAR  due  to  their  high  power  den¬ 
sity,  low  cost,  and  the  absence  of  outgassing  during  dis¬ 
charge  cycles.  NASA  G-313  used  a  one  layer  battery  of  6© 
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Figure  1-4.  Schematic  of  the  Space  Thereoacoustic 
Refrigerator  in  its  Get  Away  Special 
Cannis ter . 
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Figure  1-5. 


Block  Diagram  of  the  Electronics  for  the 
Space  Thermoacoustic  Refrigerator. 
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cells  providing  680  watt-hours  of  energy  and  weighing  about 
80  pounds  (including  the  cells'  support  structure).  We  will 
be  using  two  battery  layers  with  as  many  cells  as  we  can,  up 
to  136.  This  will  give  us  a  minimum  of  680  watt-hours  and  a 
maximum  of  1,360  watt-hours  of  available  electrical  energy, 
and  a  battery  weight  between  80  and  160  pounds.  The  number 
of  battery  cells  we  can  use  will  depend  on  the  total  weight 
of  the  driver-resonator  assembly  and  its  auxiliary  equipment 
(vacuum  can,  gas  reservoir,  etc.)  and  electronics. 

The  use  of  the  GAS  can  imposes  certain  restrictions. 
Since  we  are  using  batteries  to  supply  the  power,  the  refri¬ 
gerator  has  to  be  energy  efficient.  Also,  the  GAS  can  setup 
requires  the  STAR  to  be  compact  and  lightweight.  Figure  1-6 
shows  a  photograph  of  Hofler’s  prototype  refrigerator.  This 
setup  is  approximately  six  feet  high.  In  comparison,  the 
maximum  payload  height  for  the  GAS  can  is  28.25  inches,  or 
less  than  2.5  feet.  These  considerations  played  a  maior  role 
in  the  choice  of  equipment  for  and  the  design  of  the  STAR. 


C .  SCOPE 

Chapter  II  discusses  basic  acoustical  and  thermodynamic 
theory  as  applied  to  STAR.  The  results  of  this  theory  deter¬ 
mine  the  driver  requirements.  The  theory  of  the  selected 
electrodynamic  driver  is  given  next.  Following  the  driver 
theory  is  a  presentation  and  discussion  of  a  computer 
program  that  models  the  coupled  system,  which  consists  of 
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Figure  I- 


Photograph  of  Hofler's  (1986) 
Prototype  Refrigerator. 
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the  electrodynamic  driver  and  the  resonator  containing  the 
gas  mixture  and  stack. 

Chapter  III  describes  the  technique  for  measuring  the 
electrodynamic  driver  parameters.  These  parameters  are  used 
as  input  for  the  driver  design,  the  computer  model,  and  the 
efficiency  calculations.  This  measurement  technique  was 
developed  using  an  Altec  290-16K  electrodynamic  driver  and 
was  applied  to  the  final  driver. 

Chapter  IV  describes  the  measurement  of  the  STaR  driver 
parameters.  The  driver  chosen  for  STAR  is  a  custom  unit, 
built  by  Harmon- JBL! t  which  uses  neodimium-iron-boron 
(NIB)  magnets.  The  STAR  driver  is  a  minor  modification  of  a 
new  line  of  compression  drivers  using  NIB  magnets  and  has 
characteristics  similar  to  the  JBL  2450 J. 

Chapter  V  discusses  the  construction  of  the  STAR  driver, 
including  the  housing  and  accessories  (reducer,  stiffener, 
bellows,  accelerometer,  microphone,  electrical  feed¬ 
throughs,  etc.)  necessary  to  integrate  it  with  the  resonator 
and  the  GAS  cannister  system.  Dimensional  shop  drawings  for 
the  parts  are  included  in  Appendix  C. 

Chapter  VI  gives  conclusions  and  recommendations  for 
further  development . 
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II.  THEORY 


A.  THERMOACOUSTIC  THEORY 

Thermoacoustic  theory  has  been  developed  in  detail  by 
Rott  (1969,  1974,  1975  and  1980)  and  adapted  to  the  thermo¬ 
acoustic  refrigerator  oy  Wheatley  et  al.  (1982,  1983a, 
1983b,  1985  and  1986),  Wheatley  and  Cox  (1935),  Hofler 
(1986)  and  Swift  (1989). 

The  space  thermoacoustic  refrigerator  (STAR)  basically 
consists  of  an  acoustic  driver  producing  sound  waves  in  a 
resonant  tube  (see  Figure  1-2).  This  resonator  is  filled 
with  a  mixture  of  helium-xenon  gas  (12.5%  xenon)  pressurised 
to  ten  atmospheres ,  and  contains  a  stack  of  plastic  plates . 
This  Chapter  will  present  a  qualitative  model  for  the  ther¬ 
moacoustic  heat  pumping  process  followed  by  a  quantitative 
development  for  heat  and  work  flow  at  a  plate.  The  Chapter 
concludes  with  the  efficiency  of  a  stack  of  plates  that  are 
very  much  shorter  than  one  quarter  of  a  wavelength. 

1 .  Qualitative  Picture 

Consider  a  parcel  of  jgas  that  moves  back  and  forth 
along  one  of  the  plates  at  the  acoustic  frequency  (see 
Figure  II-l).  As  it  moves,  the  gas  parcel  will  experience 
changes  in  temperature  and  volume.  Part  of  the  temperature 
changes  come  from  the  adiabatic  compression  and  expansion  of 
the  gas  by  the  sound  pressure,  and  part  as  a  consequence  of 
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Figure  II-l.  Gas  Parcel  Diagram.  (Note:  Displacement 
in  the  Vertical  Direction  is  Shown  for 
Clarity — Parcel  Moves  in  Horizontal 
Direction  Only.) 


the  local  temperature  of  the  plate  itself.  A  temperature 
fradient  may  develop  alone  the  plate  as  a  result  of  the 
operation  of  the  refrigerator.  The  temperature  and  volume 
changes  can  be  described  by  six  separate  steps  (the  changes 
in  volume  are  indicated  by  the  change  in  size  of  the  square 
parcel  in  the  figure) . 

Assume  the  plate  is  at  a  uniform  temperature,  T.  The 
parcel  starts  at  position  1  (X  ■  0)  with  a  temperature  of  T. 
The  sound  wave  moves  the  parcel  to  the  right  to  position  2. 
The  parcel  has  undergone  adiabatic  compression  and  its  tem¬ 
perature  is  now  T++ .  Since  the  temperature  of  the  gas  is 
higher  than  that  of  the  plate,  heat  (Q)  will  flow  from  the 
gas  parcel  to  the  plate.  The  parcel’s  volume  is  decreased 
and  its  temperature  lowered  to  T+.  This  heat  flow  also 
causes  the  plate's  temperature  to  increase  at  the  position 
X+.  Parcel  position  3  is  actually  in  the  same  location  as 
position  2,  but  it  is  displaced  vertically  in  the  figure  for 
clarity.  The  parcel  now  moves  to  position  4  and  expands 
adiabatically  to  a  new  temperature  of  T-  A  repeat  of  this 
process  puts  the  parcel  in  position  5  with  a  temperature  of 
T  .  Now  the  temperature  of  the  gas  is  lower  than  that  of 
the  plate  and  heat  (Q')  flows  from  the  plate  to  the  gas 
parcel,  expanding  the  parcel  and  raising  its  temperature  to 
T__.  This  flow  causes  the  plate's  temperature  to  decrease  at 
position  X-.  The  parcel  moves  to  the  right  again  under 
adiabatic  compression  and  we  are  back  where  we  started  in 
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position  1  with  temperature  T  ( adapted  from  Wheatley,  at 
ml.,  1985). 


This  css  parcel  cycle  is  repeated  at  the  resonator’s 
operating  frequency.  For  the  STAR  this  frequency  is  between 
approximately  250  and  600  cycles  per  second,  depending  on 
the  gas  mixture  (specifically  the  speed  of  sound  in  the  gas 
mixture)  and  the  length  and  shape  of  the  resonator. 

Notice  that  the  temperature  of  the  gas  parcel  at 
position  X  *  0  is  different  depending  on  which  direction  the 
gas  parcel  is  moving.  It  is  this  phase  shift  in  temperature 
relative  to  motion  that  produces  the  thermoacoustic  effect, 
as  we  will  show  in  the  calculations  that  follow. 

It  is  also  important  to  note  that  in  order  for  heat 
to  flow  between  the  gas  parcel  and  the  plate,  the  parcel 
must  be  vertically  located  within  about  a  thermal  penetra¬ 
tion  depth  of  the  plate.  The  thermal  penetration  depth  (8k) 
the  distance  that  heat  can  diffuse  through  the  fluid 
dv  ng  a  time  1/®,  where  ®  is  the  acoustic  angular 
frequency. 

Wheatley  et  al.  (1986)  describe  thermoacoustic  en- 
ginv  as  consisting  of  long  trains  of  these  gas  parcels,  all 
about  a  thermal  penetration  depth  from  the  plate.  The  par¬ 
cels  draw  heat  from  the  plate  at  one  extreme  of  their  oscil¬ 
latory  motion  and  deposit  heat  at  the  other  extreme.  Adja¬ 
cent  heat  flows  cancel  except  r_t  the  ends .  The  net  result  is 
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that  an  amount  of  heat  0  is  passed  from  one  end  of  the  plate 
to  the  other. 

2.  Slnilf  Pim 

Swift  (1989)  illustrated  the  basic  principles  of  the 
theraoacoustic  engine  by  examining  a  simple  example  of  a 
single  plate  in  a  gas-filled  tube  such  as  we  described 
above.  In  this  example  the  acoustic  and  thermodynamic  ef¬ 
fects  are  nearly  distinct.  In  the  absence  of  the  plastic 
p}.ate,  the  plane  standing  wave  sustained  in  the  tube  has 
well-known  acoustic  properties.  It’s  pressure  and  x-direc- 
tion  velocity  (along  the  length  of  the  tube)  are  given  as 
(using  peak  values): 

P*  =  p*  sin(x/fc)  (II-l) 

and  ux  •  i(p*/>„a)  cos(x/*)  (II-2) 

where  Pm  is  the  pressure  amplitude  at  the  pressure  antinode, 
the  reduced  wavelength  *  »  h/2 *  »  a/«,  is  the  wavelength, 
a  is  the  speed  of  sound  in  the  gas,  m  is  the  angular  fre¬ 
quency,  and  pm  is  the  mean  density  of  the  gas.  The  subscript 
1  indicates  the  first  order  time  oscillation,  and  the  sub¬ 
script  m  indicates  the  mean  value.  The  total  pressure  can  be 
given  as 

P  =  P-  +  Pie‘~V  ( II-3) 

A  similar  expression  represents  the  total  velocity.  Follow¬ 
ing  Swift,  we  will  assume  that  the  first  order  in  the  acous¬ 
tic  amplitude  is  sufficient  for  all  acoustic  and  thermo- 
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dynamic  variables  (pressure,  velocity,  tesperature,  density, 
and  entropy) . 

The  sound  wave  Is  adiabatic  and  has  an  oscillatory 
temperature  given  by: 

T*.  «  (T«.a/r.c^)  p»,  (H-4) 

where  g  *  -  (dr /dT )•»/?*»  is  the  isobaric  thermal  expansion 
coefficient  and  c»  is  the  isobaric  (constant  pressure)  heat 
capacity  per  unit  aass.  Note  that  T*  and  pa  are  in  phase. 
Por  ideal  gases,  theraodynaaics  tells  us  that 

T„0/f*,C_  -  ( g-I  )  tm/tPm,  (  II-5) 

where  %  «  the  ratio  of  isobaric  to  isochorlc  specific  heats 

(5/3  for  monatomic  gases,  smaller,  but  greater  than  one,  for 
other  gases).  Combining  Equations  II-4  and  II-5  gives  us: 

Ta/T«  -  [<*-l)/<]  p*/p„.  (II-6) 

The  introduction  of  the  plastic  plate  into  the 

standing  wave  modifies  the  original,  unperturbed  temperature 
oscillations.  This  modification  is  due  to  the  heat  flow 
between  the  gas  and  the  plate,  as  described  previously.  The 
temperature  is  modified  in  both  magnitude  and  phase  for  gas 
about  a  thermal  penetration  depth  away  from  the  plate. 
According  to  Swift,  this  results  in  two  important  effects: 
1)  a  time-average  heat  flow  near  the  surface  of  the  plate, 
along  the  direction  of  acoustic  vibration,  and  2)  the 
generation  or  absorption  of  real  acoustic  power  near  the 
surface  of  the  plate.  In  refrigerators  the  effect  is  an 
absorption  of  acoustic  power .  In  prime  movers  the  acoustic 
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power  is  generated  by  the  temperature  gradient  present  in 
the  plate. 

There  are  several  assusptions  we  will  make  to 
simplify  calculations,  again  following  Swift.  We  assume  that 
the  plate  is  short  enough  compared  to  a  reduced  wavelength 
<Ax  €  X)  and  far  enough  from  both  velocity  and  pressure 
nodes  that  px  and  ua  can  be  considered  uniform  over  the 
entire  plate.  We  assume  the  gam  has  sero  viscosity,  so  that 
ui  does  not  depend  on  y  (wp  had  already  made  this  assumption 
by  looking  at  u*  in  the  x-direction  only) .  We  assume  that 
the  plate  has  a  large  enough  heat  capacity  per  unit  area 
that  its  temperature  does  not  change  appreciably  at  the 
acoustic  frequency.  We  assume  that  the  plate  has  a  mean 
temperature  gradient  in  the  x-direction  vT*,.  Finally,  we 
neglect  the  plate's  and  gas*  thermal  conductivity  in  the  x- 
direction . 

Applying  these  assumptions,  we  see  that  the  mean 
fluid  temperature  (T„(x))  is  the  same  as  that  of  the  plate. 
Swift  calculated  the  oscillating  fluid  temperature  using 
the  general  equation  of  heat  transfer.  He  kept  only  first- 
order  terms,  neglected  thermal  conduction  along  x,  and 
applied  the  boundary  condition  Tt(0)=0  imposed  by  the  plate. 
The  resulting  equation  is: 

T*(y)*C  (T«e/p«.c^)p4-(vT»/®)u»]  Cl-exp[-(l+i)y/S]  ]  .  (II-7) 
To  interpret  this  equation,  we  look  at  it  in  the  limit  that 
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the  gas  la  far  enough  from  tha  Plata  (y*8x)  to  aaka  negli- 
glbla  tharaal  contact  with  tha  plata.  Thla  glvaa: 

T»(y)->t(T„g/Fll,c^)p»-(TT«/«)u»]  (11-8) 

Tha  first  tars  in  Equation  1 1-8  la  tha  taaparatura 

oscillation  dua  to  tha  adiabatic  coaprasslons  and  axpansions 

of  tha  gas  (saa  Equation  11-4).  Tha  sacond  tars  coaas  froa 

tha  saan  taaparatura  gradiant  in  tha  gas.  As  tha  gas 

osclllatas  along  x  with  displacaaant  aaplituda  u»/»,  tha 
taaparatura  at  a  givan  point  in  spaca  osclllatas  by  an 

aaount  vT«  ux/»  even  if  tha  taaparatura  of  a  givan  place  of 
fluid  raaains  constant.  Tha  actual  taaparatura  oscillations 
ara  a  linear  superposition  of  these  two  affects. 

The  y  dependant  part  of  aquation  II-7  is  complex.  It 
approaches  1  for  y  »  8k,  as  given  previously.  It  approaches 
zero  for  y  «  8k,  where  the  plate  imposes  tha  condition 
T*.  «  0.  Moat  importantly,  for  y  *  8k,  its  magnitude  is  still 
of  the  order  of  1,  but  it  has  a  substantial  imaginary  part. 
This  phase  shift  in  the  oscillating  temperature  of  the 
standing  wave  at  y  ar  8k,  due  to  the  thermal  presence  of  the 
plate,  is  an  important  result  because  it  leads  directly  to 
the  time-averaged  heat  flow  in  the  x-direction.  It  is  the 
same  phase  shift  we  emphasized  in  the  description  of  the  gas 
parcel  motion. 

a .  Heat  Flow 

Swift  argues  that  since  we  are  neglecting  ordi¬ 
nary  thermal  conductivity  in  the  x  direction,  the  only  way 
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heat  can  ba  transported  alone  x  is  by  the  hydrodynamic 
transport  of  entropy,  carried  by  the  oscillatory  velocity 

Ut : 

a*  ■  T»  Fai  STQT  •  1/2  P*Cp  IB  [T»  ]  U»  .  (II— 2) 

The  line  above  the  quantity  s*Ui  indicates  the  tise-averace 

of  the  first-order  entropy  and  velocity  product.  The  heat 
flow  is  a  second  order  quantity,  signified  by  the  subscript 
2. 

The  total  heat  flow  6*  along  the  plate,  in  the  x- 
direction,  is  found  by  integrating  q«  over  the  y-z  plane: 

da  »  rr  oX*q«  dy  ( 11-10) 

where  IT  is  the  periseter  of  the  plate  in  the  y-z  plane. 

Substituting  for  Tt  and  performing  the  integration  gives: 

O*  «  -1/4  ir«K(M)paU»(r-l)  (11-11) 

where  ITAk  is  the  thersodynasically  active  area  in  a  plane 
perpendicular  to  the  longitudinal  acoustic  aotlon,  T«J  is 
the  heat  parameter  of  the  gas,  and  T  is  defined  as  the  ratio 
of  actual  temperature  gradient  to  the  critial  temperature 

gradient  (vT/vT««-*%)  .  The  critical  temperature  gradient 
occurs  when  the  temperature  change  along  the  plate  just 
matches  the  temperature  change  due  to  the  adiabatic 
compression  of  the  gas,  and  no  heat  flows  between  the  gas 
and  the  plate.  It  is  the  boundary  between  the  refrigerator 
and  prime  mover  functions  of  the  thersoacoustic  engines. 

Equation  11-11  shows  that  when  r  <  1,  heat  flows 
up  the  temperature  gradient  from  cold  to  hot  and  work 
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( acoustic  power )  is  absorbed,  as  for  a  refrigerator .  When 
f  ■  l  there  is  no  hast  flow.  Whan  r  >  l ,  hast  flows  down  tha 

taaparatura  gradient  froa  hot  to  cold  and  acoustic  powar  is 
producad,  as  for  a  priaa  aovar. 

Nota  that  tha  total  haat  flow  is  proportional  to 
tha  araa  ITSk  ,  and  to  T*J  (*1  for  idaal  gaaas).  It  is  also 
proportional  to  tha  product  ptu* ,  and  so  aquals  zero  if  tha 
Plata  is  at  aithar  a  pressure  node  or  a  velocity  node  of  tha 
standing  wave.  Tha  aaxiaua  value  of  tha  product  occurs 
halfway  between  tha  nodes.  Finally,  tha  haat  flow  is 
proportional  to  tha  taaparatura  gradient  factor  r-1.  For 
▼Tm  >  vT«r-*.%.  f-l  >  0  and  tha  haat  flow  is  toward  tha  pres¬ 
sure  node,  while  for  *Tm  <  r-i  <  0  and  tha  haat  flow 

is  away  froa  tha  pressure  node.  If  suitable  haat  exchangers 
at  teaperatures  T«  and  T©  are  installed  at  the  ends  of  the 

plate  (with  -  T0  «  vT«  Ax),  this  haat  flow  carries  haat 

froa  one  exchanger  to  the  other. 

The  heat  flow  is  snail  under  ordinary 
circuastances.  However,  in  closed  resonators  it  is  possible 
to  achieve  sound  aaplitudes  many  orders  of  aagnituda  higher 
than  those  of  ordinary  conversation.  Since  0*  is 

proportional  to  p*»,  and  since  in  practical  acoustic  engines 
the  entire  cross  section  of  the  standing  wave  is  filled  with 

plates  (spaced  roughly  apart),  very  high  heat  flows 
and/or  a  large  TM  -  t©  aay  be  achieved. 
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b.  Work  Flow 


The  work  flow  (l.e.  acoustic  power)  is  given  by 
the  work  per  cycle  times  the  rate  at  which  that  work  occurs 
(the  acoustic  frequency  f ) .  From  thermodynamics,  the  average 
acoustic  power  produced  per  unit  volume  is : 

**  =  - ( «•/ Pm )  TpTFT  =  -  1/2  u>0  IratTx]  .  (11-12) 

The  gas  about  a  thermal  penetration  depth  away 

from  the  plate  "breathes,”  because  of  thermal  expansion  and 
contraction,  with  the  right  time  phase  with  respect  to 
oscillating  pressure  to  do  (or  absorb)  net  work.  This  is 
exactly  the  same  gas  that  we  have  seen  is  responsible  for 
the  heat  flow.  Gas  elsewhere  is  ineffective  in  doing  (or 
absorbing)  work.  The  density  oscillations  for  y  ^  Sk  and  for 

y  »  Sk  are  in  phase  with  the  pressure  oscillations,  and 
hence  do  (or  absorb)  no  net  work. 

The  total  acoustic  power  Ws*  produced  is  found  by 
integrating  fc*  over  all  space,  as  with  the  heat  flow: 

W*  =  1/4  TTSk  U-l)  <p**/pma)  (T-l)  Ax/A,  (11-13) 
where  all  the  terms  have  been  defined  previously. 

The  acoustic  power  is  proportional  to  the  volume 
ITSkAx  of  fluid  that  is  about  a  thermal  penetration  depth 
from  the  plate.  It  is  proportional  to  p4*,  and  so  is  quad¬ 
ratic  in  the  acoustic  amplitude  (as  was  the  heat  flow)  and 
vanishes  at  pressure  nodes.  Finally,  Wa  is  proportional  to 

(T-l) ,  the  same  temperature  gradient  factor  as  appeared  in 
Gz  (Equation  11-11).  When  vTm  =  vTcrl*,  (T-l)  =  0,  and  there 
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are  no  temperature  oscillations  in  the  fluid  other  than 
those  due  to  adiabatic  processes,  and  no  acoustic  power  is 

absorbed  or  generated.  For  9xm  >  vT,,.**,  r  -1  >  0  and  acous¬ 
tic  power  is  produced  near  the  plate.  Whether  this  power 
increases  the  amplitude  of  the  standing  wave,  is  radiated 
away  to  infinity,  is  simply  absorbed,  or  flows  through  an 
acoustic-to-electric  transducer  to  generate  electric  power 
depends  on  details  of  the  resonator,  not  on  the  plate  itself 

or  on  the  standing  wave  near  the  plate.  For  vT»  <  vTCi-t<»,  T- 
1  <  0  and  acoustic  power  is  absorbed  near  the  plate.  For  a 

tube  without  plates  which  has  a  diameter  less  than  the 
wavelength,  at  constant  temperature  (vT  *  0),  this  work  flow 
is  responsible  for  the  ordinary  thermal  attenuation  of  a 
sound  wave  (Kinsler,  et  al.,  1982). 
c.  Efficiency 

We  combine  Equations  11-11  and  11-13  to  get  the 
efficiency  of  the  plate  with  no  viscous  or  longitudinal 
conduction  losses : 

n  =  =  U-l)/(T,*e)  (®  Ax  p»)/(p„a*u*)  .  (11-14) 

Since  ux  =  u*»  sin  x/fc,  and  p»  *  pma  Uo  cos  x/A, 

n  *  <*-l)/(T«g)  (Ax/fc  tan  x/A) .  (11-15) 

For  x  «  n  =  U-l)/(T*,e)  (Ax/x),  (11-16) 

and  n  =  na.rnot/r  (11-17) 

This  efficiency  approaches  the  Carnot  efficiency 
as  the  power  output  and  heat  transfer  rates  approach  zero. 
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Me  can  make  a  similar  calculation  for  the 
refrigerator  mode  of  operation,  where  the  relevant 
efficiency  is  the  coefficient  of  performance,  COP  * 

Me  find  that 

COP  *  f  C0P.«^„o*.  (11-18) 

3.  Short  Stack 

After  developing  the  basic  principles  of  the  thermo¬ 
acoustic  engine  using  a  simplified  example,  Swift  developed 
a  more  realistic  model  by  including  viscosity,  longitudinal 
thermal  conductivity,  finite  (instead  of  infinite)  plate 
heat  capacity,  and  many  plates.  He  still  made  the  assump¬ 
tions  that  the  temperature  spanned  is  much  less  than  the 
absolute  temperature,  that  the  length  of  the  plates  is  much 
shorter  than  the  reduced  wavelength,  and  that  the  spacing 
between  the  plates  is  greater  than  a  few  penetration  depths. 

Using  the  equation  of  motion  for  the  gas,  the  boun¬ 
dary  condition  that  the  velocity  at  the  gas-plate  interface 
is  zero,  the  continuity  equation  of  the  gas,  the  heat  flows 
in  the  gas  and  plates,  and  basic  thermodynamic  properties, 
Swift  derived  a  wave  equation  for  pressure  p*(X)  in  terms  of 
dT»/dx,  material  properties  and  geometry.  Swift  also  derived 
an  equation  for  enthalpy  flow  along  x  in  terms  of  p*(x), 
T«(x),  geometry,  and  material  properties  (including  Prandtl 
number  of  the  gas  cr  =  c^u/k,  where  u  is  the  viscosity,  and  k 
is  the  thermal  conductivity  of  the  gas). 
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In  the  boundary  layer  approximation,  for  a  short 
(Ax  *  *)  stack  of  plates,  and  neglecting  viscosity,  the  heat 
flux  and  acoustic  power  are  given  by  expressions  very  simi¬ 
lar  to  the  simple  results  of  the  single  plate  example. 
Practical  engines  can  be  expected  to  have  all  available 
cross-sectional  area  filled  with  plates  spaced  about  48k 
apart.  The  only  effect  of  longitudinal  (i.e.,  along  x) 
thermal  conductivity  is  to  add  to  the  heat  flow  in  the  x 
direction  by  simple  conduction.  For  gases  this  effect  is 
normally  negligible. 

The  inclusion  of  viscosity  adds  considerable 
algebraic  complication,  and  a  little  conceptual  difficulty 
as  now  u>  is  a  function  of  y. 

Without  viscosity,  vT«r>&*  is  the  temperature 
gradient  for  which  the  temperature  oscillation  at  a  given  x- 
location  is  zero.  But  in  the  presence  of  viscosity,  Tx 
depends  on  x  as  well  as  y,  so  no  well-defined  vT«rl*  exists. 

To  the  lowest  order  in  VF  (where  a  is  the  fluid's 
Prandtl  number) ,  the  hydrodynamic  heat  flow  is  increased  by 

1/(1— /a).  This  factor  arises  because  the  mean  velocity  <u,> 
underestimates  the  velocity  with  which  entropy  is  convected 
along  the  stack.  As  we  saw  in  the  single  plate  calculations, 
the  convective  entropy  transport  occurs  mostly  at  a  distance 
6k  from  the  plate.  The  velocity  there  is  higher  than  the 
mean  velocity  by  1/d-Vtf). 
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Th«  lowest  order  effect  of  viscority  on  the  acoustic 
power  Is  the  power  dissipated  by  viscous  shear  in  the  fluid 
within  the  boundary  layer,  a  well-known  fluid  mechanical 
result . 

For  accurate  calculations,  the  full  general  equation 
must  be  used;  but  for  rough  estimates,  the  results  presented 
here  are  good  approximations,  and  are  much  simpler  to 
compute  with.  In  addition,  the  results  presented  here  are 
fairly  easy  to  interpret  physically.  The  expressions  for 
arbitrary  viscosity  are  complicated,  but  those  to  lowest 
order  in  viscosity  are  simple.  It  is  easy  to  see  the  single 
plate  expressions  for  heat  flow  and  acoustic  power  about  a 
thermal  penetration  depth  from  the  plates,  modified  as 
longitudinal  thermal  conductance  adds  to  the  heat  flow,  as 
viscous  shear  dissipates  acoustic  power,  and  as  non- ideal 
plate  properties  modify  the  thermal  boundary  condition  at 
the  fluid-plate  interface. 

4.  Design  Considerations 

One  of  the  important  considerations  in  the  design  of 
a  thermoacoustic  refrigerator  is  the  location  of  the  stack 
of  plates  with  respect  to  the  pressure  and  velocity  nodes 
and  antinodes  of  the  acoustic  standing  wave  in  the  resonant 
tube.  Experiments  performed  by  Wheatley,  et  al.  (1986) 
showed  that  the  heat  flow  as  a  function  of  position  fit  a 
simple  sine  curve  whose  spacial  period  is  half  the  wave¬ 
length  of  the  acoustic  standing  wave.  By  noting  how  the  sign 
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of  th«  temperature  difference  varied  with  respect  to  the 
plate's  location  in  the  sound  wave,  they  saw  that  heat 
always  flows  in  the  direction  of  the  closest  pressure  anti- 
node.  This  effect  is  expected  because  a  parcel  of  gas  moving 
in  the  direction  of  a  pressure  antinode  is  compressionally 
warmed  and  will  transfer  heat  to  the  plate.  A  parcel  moving 
toward  a  pressure  node  is  cooled  by  expansion  and  will  draw 
heat  from  the  plate.  At  both  the  pressure  nodes  (velocity 
antinodes)  and  antinodes  (velocity  nodes)  heat  flow  drops  to 
zero.  Thus  the  acoustic  heat  flow  depends  on  both  the  acous¬ 
tic  pressure  and  the  fluid  velocity. 

As  a  plate  or  stack  of  plates  is  moved  away  from  the 
pressure  antinode,  the  temperature  gradient  developed 
becomes  smaller.  At  a  quarter  of  a  wavelength,  no  gradient 
forms  (the  critical  temperature  gradient  equals  zero) .  This 
positioning  effect  is  important  in  the  design  of  a  refriger¬ 
ator,  because,  together  with  the  length  of  the  plates,  it 
places  an  upper  limit  on  the  maximum  temperature  drop 
possible  across  the  stack. 

B.  DRIVER  REQUIREMENTS 

In  order  to  drive  the  thermoacoustic  refrigerator  we 
needed  a  device,  a  transducer,  that  would  convert  the  elec¬ 
trical  energy  of  a  battery  pack  to  acoustic  energy,  as 
discussed  in  Subsection  I. A. 3.  We  wanted  as  high  a  pressure 
(or  force)  and  volume  velocity  (or  motion)  as  possible  to 
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get  the  aost  acoustic  power  provided  to  the  system.  There 
were  several  types  of  electroacoustic  transducers  to  choose 
from:  electro-static,  electret,  electrodynamic,  magneto- 
strictive,  electro-strictive,  and  piezoelectric.  Only  one  of 
these,  the  e.lectrodynamic  transducer,  can  provide  our  system 
with  both  the  large  forces  and  velocities  necessary.  This 
transducer  will  be  described  in  the  next  section.  The 
remainder  of  this  section  describes  the  rest  of  the  trans¬ 
ducers  and  explains  why  each  was  not  appropriate  for  use  in 
the  thermoacoustic  refrigerator. 

The  electrostatic  transducer  consists  of  a  pair  of 
charged  electrodes ,  or  capacitor  plates ,  one  of  which  is 
held  stationary  while  the  other,  the  diaphragm,  moves  in 
response  to  electrical  excitation.  The  electret  is  a  kind  of 
electrostatic  transducer  that  has  a  diaphragm  of  polarized 
plastic,  and  therefore  does  not  need  an  external  voltage 
supply  to  the  diaphragm  (Kinsler,  et  al,  1982,  p.  350).  If 
an  alternating  voltage  is  applied  across  the  plates,  the 
diaphragm  moves,  thereby  radiating  an  acoustic  wave. 
Unfortunately  there  is  very  little  force  (or  pressure)  and 
displacement  associated  with  this  motion.  This  type  of 
transducer  may  be  useful  as  a  microphone  but  is  not  an 
appropriate  choice  to  drive  the  thermoacoustic  refrigerator. 

In  a  magnetostrictive  transducer  a  change  in  the 
magnetic  polarization  in  the  material  causes  an  elastic 
strain  (Wilson,  1985,  p.  3.).  Although  the  forces  are  more 
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than  adequate  the  strain  produces  too  little  notion  to  be 
useful  to  us. 

In  electrostrictive  and  piezoelectric  transducers  an 
externally  applied  electric  field  causes  a  change  in  the 
dielectric  polarization  of  the  material,  which  in  turn 
causes  an  elastic  strain  (Wilson,  1985,  p.  3).  Once  again 
there  is  sufficient  force  but  the  displacement  due  to  the 
strain  is  too  snail  to  be  useful. 

Because  these  techniques  have  sufficient  force,  it  would 
be  possible  to  overcoae  the  displacement  limitation  using 
mechanical  leverage  or  other  forms  of  mechanical  impedance 
transformation,  but  this  was  not  considered  due  to  the 
associated  increase  in  device  complexity. 

C.  ELECTRODYNAMIC  THEORY 

i-  Physiss 

The  electrodynamic  driver  is  a  moving  coil 
transducer  ^Wich  converts  electrical  energy  to  kinetic  ener¬ 
gy  (motion) .  It  consists  of  a  diaphragm  attached  to  a  cylin¬ 
drical  coil  of  wire  (voice  coil)  that  is  suspended  in  a 
fixed  magnetic  field  (see  Figure  II-2) .  If  an  alternating 
current  is  supplied  to  the  coil,  the  interaction  of  the 
current  s--1  the  rr,.. ..tic  field  will  induce  a  force  on  the 
coil  so  that  the  diaphragm  moves.  The  magnitude  of  the  force 
is  F  =  B1I,  where  B  .is  the  magnitude  of  the  magnetic  induc¬ 
tion,  1  is  the  len>  .  of  the  coil  and  I  is  the  current.  The 
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Figure  II-2.  Schematic  of  an  Electrodynamic  Transducer 
Showing  (A)  Magnet,  (B)  Back  Plate, 

(C)  Pole  Piece,  (D)  Acoustic  Seal, 

(E)  Diaphragm,  (F)  Voice  Coil, 

(G)  Surround  (Suspension),  and 

(H)  Electrical  Lead  Connection. 


direction  of  the  force,  given  by  the  cross  product  of  the 
current  vector  and  the  magnetic  field  (fros  Ampere’s  law  £  « 
li.  X  g) ,  is  everywhere  orthogonal  to  the  current  and  the 
radial  eagnetic  field.  The  notion  of  the  diaphragm  produces 
sound  waves  with  a  frequency  equal  to  that  of  the  alterna¬ 
ting  current.  Another  important  property  of  the  electrodyna- 
nic  transducer  is  that  the  notion  of  the  voice  coil  induces 
an  electromotive  force  (emf)  which  equals  B1  times  the 
velocity  of  the  voice  coil  (emf  ■  Blv) .  This  ’’back  emf" 
increases  the  electrical  impedance  of  the  driver  and  can  be 
used  to  monitor  the  driver  motion  as  discussed  in  Section 
III.B. 

The  driver  is  attached  to  the  resonator  via  a 
reducer  and  bellows.  The  reducer  is  a  cone-shaped  piece  of 
aluminum  mechanically  connecting  the  driver  voice  coil  which 


produces 

the 

force  to 

the  bellows ,  which 

in 

turn  is 

in 

contact 

with 

the  gas 

in  the  resonator . 

The 

reducer 

is 

designed  to  be  lightweight,  rigid,  and  strong  for  a  direct 
transfer  of  motion  and  energy  from  the  moving  coil  to  the 
bellows.  The  bellows  is  a  lightweight  and  flexible  gastight 
seal  between  the  driver  and  the  tube.  An  aluminum  plate  is 
used  as  an  interface  between  the  reducer  and  the  bellows  to 
keep  the  bellows  face  rigid.  For  the  STAR  driver  the  reducer 
and  plate  are  incorporated  into  a  single  unit  machined  from 
aluminum  bar  stock. 
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Figure  I I -3  shows  s  mechanical  model  of  the 
electrodynamic  transducer  as  a  forced  simple  harmonic 
oscillator  (SHO),  which  consists  of  a  moving  mass  (n,) , 
stiffneau  (k),  mechanical  damping  or  resistance  (R„) ,  and 
driving  force  (F) .  Included  in  the  moving  mass  are  the 
masses  of  the  voice  coil,  reducer,  stiffener,  and  bellows. 
The  stiffness  includes  the  suspension  stiffness  of  the  voice 
coil  surround,  the  stiffness  of  the  gas  volume  trapped 
behind  the  diaphragm,  and  the  stiffness  of  the  bellows.  The 
equations  governing  the  SHO  can  thus  be  used  to  determine 
the  moving  mass  and  stiffness  of  the  driver.  The  primary 
relationship  used  is: 

*•*  »  k/m.,  (11-19) 
or,  the  resonant  angular  frequency  squared  equals  the 
stiffness  divided  by  the  moving  mass.  The  use  of  this  and 
other  relationships  will  be  discussed  further  in  Chapter 
III,  Driver  Parameter  Measurements. 

2.  gamHflians.  Sicsaii 

An  electrical  circuit  can  be  drawn  that  is 
equivalent  to  an  electromechanical  system  if  the  values 
given  for  Rm,  Lm,  and  CM  in  the  circuit  produce  the  same 
electrical  effects  as  the  electromechanical  system  itself. 
In  the  case  of  the  electrodynamic  driver,  the  equivalent 
circuit  is  shown  in  Figure  I I -4  with  the  values  givan  by 
Kinsler,  et  al.  (1982)  as: 
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Rm  -  (Bl)*/R«, 

(11-20) 

L«  -  (Bl)*/k, 

(II-2i) 

Cm  -  Mm/  ( BJ  )  * , 

(11-22) 

where  B1  is  the  transduction  coefficient, 

Rm  is  the 

mechanical  resistance,  and  k  and  a.  are  the  stiffness  and 
eovins  sass  as  described  previously. 

The  properties  of  the  equivalent  circuit  allow  us  to 
evaluate  the  performance  of  the  electrodynaeic  driver  as  if 
it  were  an  entirely  electrical  aystea.  The  amount  of  power 
the  driver  generates  will  depend  on  the  load  it  secs.  Since 
one  of  the  primary  purposes  of  this  thesis  is  to  maximize 
the  ratio  of  the  acoutic  power  delivered  to  the  electrical 
power  provided  (the  electroacoustic  efficiency),  it  is 
necessary  to  provide  a  load  to  the  system  in  order  to 
evaluate  driver  performance.  This  load  is  originally 
represented  as  a  complex  electrical  impedance  in  the 
equivalent  circuit.  The  actual  value  of  the  load  will  depend 
on  the  design  of  the  resonator  but  can  be  characterized  by 
an  acoustic  impedance  which  can  be  transformed  to  an 
equivalent  electrical  load  impedance. 

3.  Losses 

In  our  search  for  the  highest  electroacoustic 
efficiency  it  is  useful  to  understand  where  possible  losses 
may  occur  and  do  what  can  be  done  to  minimize  them.  Losses 
that  are  an  intrinsic  part  of  the  electrical  driver  are  eddy 
current  losses  in  the  magnetic  structure  generated  by  the 
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current  in  the  voice  coil,  Joule  heating  losses  due  to  the 
drive  current  dissipated  by  the  voice  coil  electrical 
resistance,  and  mechanical  losses  due  to  the  internal  fric¬ 
tion  of  the  surround  and  the  viscosity  of  the  gas  in  the 
"gap."  These  losses  need  to  be  taken  into  account.  Very 
little  can  be  done  about  the  eddy  current  losses  and  the 
viscosity  of  the  gas  which  is  necessary  to  improve  heat 
removal  from  the  voice  coil.  The  losses  due  to  surround 
(suspension)  mechanical  resistance  can  be  controlled  by 
choice  of  the  materials  and  removal  of  extraneous  suspension 
material  (see  Section  III.H). 

D .  COUPLED  SYSTEM 

1 .  Computer  Model 

Any  discussion  of  the  driver  section  of  the 
thermoacoustic  refrigerator  is  incomplete  without  some 
understanding  of  how  the  driver  parameters  are  coupled  to 
the  resonator  (acoustic  load).  The  equivalent  electro¬ 
mechanical  circuit  for  the  coupled  system  is  the  same  as 
that  shown  in  Figure  II-4  with  the  acoustic  load  impedance: 

a  — p«»csj  [  (l  +  j«/k)/(l  +  («/'k)*)  *  ( 11-23 ) 

( (cos (kL)sin(kL)  +  j  sinh(«L)cosh(aL) )  * 

(sina(iCL)COsh*(«L)  +  cos2( kL) sinha( <*L) )  )  ]  , 
where  P<e,  is  the  density  of  the  gas,  c  is  the  speed  of  sound 

(labeled  "a"  in  previous  discussions),  s  is  the  cross- 

sectional  area  of  the  resonator,  k  is  the  wave  number 
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(*  «/c),  «  *  ®«./(2Qc),  Q  is  the  quality  factor  of  the  reso¬ 
nator,  and  L  is  the  effective  length  of  the  resonator 
(Kinsler,  et  al.,  1982).  This  circuit  combines  the  equiva¬ 
lent  circuit  of  the  electrodynamic  transducer  with  the  com¬ 
plex  impedance  load  representing  the  resonator.  One  way  to 
examine  this  coupling  is  through  computer  modeling. 

LT  Michael  P.  Susalla,  USN,  developed  a  computer 
program  which  models  the  thermoacoustic  refrigerator  as  the 
equivalent  electrical  circuit.  The  program  models  each  of 
the  mechanical  and  acoustical  components  as  equivalent 
electrical  components,  as  discussed  above.  Using  the 
electrical  equation  from  the  combined  circuit  we  determine 
the  system  performance  as  different  parameters  are  varied. 

a.  Input 

There  are  two  sets  of  input  parameters  needed 
-for'  the  model:  (1)  driver  parameters  and  (2)  resonator 
parameters.  The  driver  electrical  parameters  include  the 
voice  coil  inductance  (LE)  and  resistance  (RE) ,  and  the 
current  (I)  supplied  to  the  voice  coil  (which  is  held  at  one 
ampere  for  all  test  cases).  The  driver  mechanical  parameters 
used  in  the  model  are  the  same  equivalent  resistance,  induc¬ 
tance,  and  capacitance  as  those  developed  for  the  equivalent 
circuit.  Therefore  the  parameter  values  needed  for  input  are 
the  transduction  coefficient  (BL),  the  mechanical  resistance 
(RM),  the  moving  mass  (M) ,  the  suspension  stiffness  (SSUP) , 
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the  volume  of  the  gas  behind  the  diaphragm  (VOL),  and  the 
cross-sectional  area  of  this  volume  (A) . 

The  resonator  parameters  needed  include  the  tube 
quality  factor  (Q) ,  the  speed  of  sound  (O  ,  the  resonance 
frequency  (FO) ,  the  density  of  the  gas  in  the  tube  (RHO) , 
the  resonator  cross-sectional  area  (S),  and  the  equivalent 
length  of  the  resonator  (L).  For  this  model  the  piston  area 
is  set  equal  to  the  resonator  cross-section  (S) .  For  the 
tests  reported  here  they  were  nearly  equal  but  could  differ 
sustantially  in  other  designs . 
b .  Output 

Four  output  plots  are  generated  by  the  graph 
program  that  accompanies  the  model:  (l)  voltage  across  the 
driver  vs.  frequency  (since  we  put  unit  current  into  the 
driver,  this  graph  can  also  be  read  as  the  input  electrical 
impedance),  (2)  piston  velocity  vs.  frequency,  (3)  acoustic 
power  delivered  to  the  load  (resonator)  vs.  frequency,  and 
(4)  electroacoustic  efficiency  vs.  frequency.  Each  of  these 
plots  should  produce  peaks  at  both  the  driver  and  the  tube 
resonant  frequencies. 

The  model  uses  the  electrical  equation  for  the 
coupled  equivalent  circuit  to  calculate  voltage  directly 
(holding  the  current  constant  at  one  ampere)  for  the  first 
output  plot.  It  also  directly  calculates  the  force  on  the 
moving  mass,  F  =  B1I.  From  this  force,  the  piston  velocity 
for  the  second  plot  is  found  by  dividing  the  force  by  the 
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mechanical  impedance  (Z«Mh)  .  Pressure  is  calculated  as  tube 
impedance  times  volume  velocity,  and  acoustic  power  is 
pressure  times  volume  velocity  times  the  cosine  of  their 
phase  difference  (Equation  1-8)  for  the  third  plot.  For  the 
final  plot  the  electrical  power  is  calculated  as  voltage 
times  current  times  the  cosine  of  their  phase  difference 
(Equation  1-7) ,  and  the  model  takes  a  ratio  of  the  electri¬ 
cal  to  the  acoustical  powers  for  the  electroacoustic  effi¬ 
ciency.  Since  the  components  of  these  powers  are  complex  and 
have  phase  differences,  the  actual  powers  are  calculated  as 
one  half  the  real  part  of  [current  times  complex  conjugate 
of  voltage] ,  and  one  half  the  real  part  of  [pressure  times 
the  complex  conjugate  of  volume  velocity] .  Each  of  the  above 
mentioned  plots  is  the  discussed  output  parameter  versus 
frequency. 

2.  Assumed  Typical  Parameters 

Table  II-l  lists  the  input  parameters  with  sample 
values  for  the  Altec  290- 16K  driver  in  air  with  a  long 
straight  resonant  tube  attached  (see  Chapter  III).  Sample 
outputs  of  the  four  plots  using  these  input  values  are  given 
in  Figures  II-5  through  II-8. 

Figure  II-S  plots  the  driver  voltage  as  a  function 
of  frequency.  Since  the  current  is  held  constant  at  one 
ampere,  this  plot  actually  gives  us  information  on  driver 
electrical  input  impedance  vs.  frequency.  Notice  that  there 
are  two  peaks  in  this  curve.  The  first  peak  occurs  at  a 
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TABLE  II-l.  MODEL  INPUT  VALUES— EXAMPLE  DATA 


Symbol 

Meaning 

Units 

Value 

RM 

Mechanical  resistance 

kg/s 

1.33 

M 

Moving  mass 

kg 

0.0151 

SSUP 

Spring  constant 

N/m 

9.2X103 

RHO 

Density  of  gas 

kg/ma 

1.2 

C 

Sound  speed 

m/s 

342.0 

S 

Resonator  X-sec .  area 

m5* 

9.35X10 

BL 

Driver  transduction  factor 

N/A 

21.5 

Q 

Resonator  quality  factor 

60.0 

FO 

Resonator  resonance  freq. 

Hz 

240.0 

A 

Back  volume  X-sec .  area 

ms 

7.92X10 

VOL 

Back  volume 

m* 

1.0X10" 

LE 

Voice  coil  inductance 

H 

1.0X10- 

RE 

Voice  coil  DC  resistance 

ft 

10.6 

I 

Driver  current 

A 

1.0 

L 

Resonator  length 

m 

0.712 

VOLTAGE  VS.  FREQUENCY 
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Figure  I 1-6.  Computer  Model  Sample  Output — Piston 
Velocity  vs.  Frequency  for  One  Ampere 
into  Driver. 
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ACOUSTIC  POWER  VS.  FR 
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ELECTROACOUSTIC  EFFICIENCY  VS. 


Figure  I I -8.  Computer  Model  Sample  Output — 
Electroacoustic  Efficiency  vs. 
Frequency . 
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lower  frequency  end  has  a  such  higher  amplitude  than  the 
second.  This  higher  amplitude  peak  occurs  near  the  driver 
resonant  frequency  <f.  .  (1/2*) [k/s.3 *  124  Hz).  It 
dominates  the  curve  because  it  is  the  driver  iapedance  that 
we  are  plotting,  and  the  voice  coil  notion,  and  consequently 
the  "back  eaf,"  is  the  largest  at  that  frequency.  The  fact 
that  there  is  a  second  peak  near  the  tube  resonance  (FO  * 
240  Hz)  shows  the  effect  of  the  coupling  of  the  two  parts  of 
the  system. 

In  the  plot  of  piston  velocity  vs.  frequency  (Figure 
I I -6),  the  dominant  peak  again  occurs  near  the  driver 
resonance  for  the  reasons  presented  above.  The  smaller 
coupled  peak  at  the  tube  resonance  is  due  to  the  change  in 
the  load  impedance  which  occurs  near  the  tube  resonance. 

Figure  II-7  shows  the  plot  of  the  acoustic  power  vs. 
frequency .  Now  the  dominant  peak  occurs  at  the  tube 
resonance  because  acoustic  power  is  primarily  a  function  of 
the  tube  parameters.  Here  the  secondary  peak  at  the  driver 
resonance  is  a  consequence  of  the  greatly  increased  piston 
velocity  at  the  driver  resonance  putting  energy  into  the 
smaller  non-resonant  acoustic  load  of  the  tube. 

In  the  plot  of  electroacoustic  efficiency  vs. 
frequency.  Figure  I I -8,  the  primary  peak  at  the  tube 
resonance  is  sharp  and  the  secondary  peak  at  the  driver 
resonance  is  small  but  visible.  It  is  important  to  note  that 
the  electroacoustic  efficiency  peaks  at  a  value  of  nearly 
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60S  •van  though  the  two  individual  resonances  are  separated 
by  nearly  an  octave  in  frequency.  The  peak  acoustic  power 

for  one  asp  as  shown  in  Figure  II-7  is  10.3  W.  The  losses 
(40%)  are  due  primarily  to  Joule  heating  in  the  voice  coil 
DC  resistance  (P«  *  1/2  (I*R*o)  *  S.3  W) .  This  accounts  for 
34%.  The  reaaining  loss  can  probably  be  attributed  to 
mechanical  resistance  losses  in  the  suspension. 

Appendix  A  contains  a  printout  of  the  calculation 
program.  Appendix  B  contains  a  printout  of  the  graph 
program . 
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III.  DRIVER  PARAMETER  MEASUREMENTS 


A.  INTRODUCTION 

In  Sections  II.B  and  II.C  we  discussed  the  choice  of  the 
electrodynssic  transducer  to  drive  the  theraoacoustic  refri¬ 
gerator.  Given  this  choice,  we  needed  to  develop  a  system 
for  measuring  the  mechanical  parameters  of  this  type  of 
driver  to  be  used  in  determining  the  electroacoustic  effi¬ 
ciency.  The  relevant  parameters  (as  discussed  in  Section 
II. D)  are:  (1)  moving  mass,  (2)  suspension  stiffness,  (3) 
transduction  coefficient  (SI),  and  (4)  mechanical  resis¬ 
tance  . 

The  driver  we  selected  for  the  preliminary  measurements 
is  the  Altec  model  290-16K  moving-coil  loudspeaker.  There 
are  several  pieces  of  equipment  needed  to  interface  the 
driver  with  the  resonator.  Three  of  these  parts — the  reducer 
cone,  stiffener,  and  dynamic  bellows — become  an  integral 
part  of  the  driver,  and  their  effects  on  the  driver  parame¬ 
ters  must  be  taken  into  account.  We  first  measured  the 
parameters  of  the  driver  alone,  then  added  the  other  parts 
and  determined  the  new  parameters  for  the  entire  system.  The 
driver-bellows  combination  was  then  fitted  to  one  end  of  a 
resonant  tube.  The  actual  resonator  used  in  the  thermo¬ 
acoustic  refrigerator  will  be  a  complex  structure  that  is 
not  easily  modeled.  For  the  purposes  of  the  driver  parameter 
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measurements  we  used  a  resonant  tuba  of  unifora  cross - 
•action  to  allow  easy  coapariaon  of  axpariaantal  valuas  to 
thaory.  Tha  non -driver  and  of  tha  closed  tuba  housed  several 
■icrophonas.  Tha  electroacoustic  efficiency  of  tha  systea  at 
resonance  was  datarainad  froa  aeasureaents  of  electric  power 
delivered  to  tha  driver  and  the  acoustic  pressure  at  tha 
microphone.  The  acoustic  power  dissipation  was  datarainad 
froa  tha  quality  factor  of  tha  tuba  at  resonance  and  tha 
aicrophone  sensitivity  (tha  relationship  between  the  pres¬ 
sure  that  tha  aicrophone  senses  and  its  voltage  output ) . 

B.  IMPEDANCE  ANALYZER 

Tha  priaary  instruaant  used  in  the  testing  and  evalua¬ 
tion  of  the  electrodynaaic  driver  was  the  Hewlett  Packard 
Model  4192A  LF  Iapedance  Analyzer.  According  to  the  manufac¬ 
turers  Operation  and  Service  Manual,  the  4192A  is  a  fully 
automatic,  high  perforaance  test  instrument  designed  to 
measure  a  wide  range  of  iapedance  parameters  as  well  as 
gain,  phase,  and  group  delay.  Its  two  measurement  display 
sections  provide  direct  readout  of  the  selected  measurement 
parameters  with  4.S  digit  resolution  along  with  the  appro¬ 
priate  units.  In  noraal  mode  operation,  the  4192A  performs 
approximately  five  aeasureaents  per  second.  It  also  provides 
an  average  measurement  node  (one  measurement  per  second)  for 
higher  resolution  and  noise  rejection.  A  high  speed  mode  can 
perform  ten  measurements  per  second. 
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The  4192A  can  provide  measuring  frequency,  OSC  voltage 
level,  and  DC  bias  voltage  ( impedance  measurements  only) 
equivalent  to  actual  operating  conditions.  The  sweep  capa¬ 
bility  of  the  built-in  frequency  synthesiser  and  DC  bias 
source  permits  quick,  convenient,  and  accurate  measurements . 
Measuring  frequency  or  DC  bias  voltage  can  be  autoaatically 
or  manually  swept  in  either  direction.  OSC  level  can  be 
manually  swept  in  either  direction  in  lmV  increments.  Actual 
test  voltage  across  or  test  signal  current  through  the 
device  under  test  is  also  measured. 

In  amplitude/phase  measurements,  the  4192A  can  measure 
four  transmission  parameters  -  gain/loss,  level,  phase,  and 
group  delay.  In  impedance  measurements,  the  4192A  can 
measure  eleven  impedance  parameters  -  absolute  value  of 
impedance,  absolute  value  of  admittance,  phase  angle,  resis¬ 
tance,  reactance,  conductance,  susceptance.  Inductance, 
capacitance,  dissipation  factor,  and  quality  factor.  All 
iapedence  measurements  have  a  basic  accuracy  of  0.1%  and  a 
resolution  of  4.S  digits.  Also,  the  circuitry  of  the  4192A 
provides  direct  and  accurate  impedance  measurements  of  both 
grounded  and  floated  devices. 

C .  MECHANICAL  PARAMETERS 

The  most  important  parameters  for  a  moving  coil  trans¬ 
ducer  are  the  moving  mass,  suspension  stiffness,  and  the  B1 
(transduction)  coefficient.  The  moving  mass  of  the  driver 
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includes  the  mass  of  the  voice  coil,  the  diaphragm,  and  the 
hydrodynamic  mass  of  the  air  being  pushed  by  the  “diaphragm. 
The  masses  of  the  reducer  and  stiffener  are  also  important, 
as  are  the  moving  mass  and  suspension  stiffness  of  the 
bellows.  The  moving  mass  and  stiffness  of  the  combined 
driver-reducer-bellows  system  is  also  important.  The  values 
for  the  combined  system  should  equal  the  sum  of  the  values 
of  the  individual  parts.  The  stiffness,  mass,  and  B1  coeffi¬ 
cient  give  us  the  information  we  need  to  construct  the 
equivalent  electrical  circuit,  as  discussed  in  Subsection 
II. C. 2.  Also,  we  have  to  insure  that  the  stiffness  of  the 
driver  system  is  just  large  enough  to  produce  a  mechanical 
resonance  at  the  tube  resonant  frequency,  at  which  point  the 
driver  piston  looks  like  a  soft  (or  open)  tube  end.  This  may 
appear  to  be  contradictory  to  having  a  velocity  node  of  the 
standing  wave  at  the  piston,  but  it  is  not.  It  is  how  we  get 
high  efficiency. 

The  important  design  parameters  of  the  tube  are  its 
resonance  frequency  and  quality  factor.  The  length  of  the 
tube  and  the  speed  of  sound  in  the  contained  gas  determine 
the  resonance  frequency. 

1 •  Dynamic  Mass  Loading  of  Driver 

We  used  the  driver's  similarity  to  a  simple  harmonic 
oscillator  to  find  the  moving  mass  and  suspension  stiffness 
of  the  driver  by  taking  advantage  of  the  fact  that  the 
resonance  frequency  squared  equals  the  stiffness  divided  by 
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the  mas* .  If  we  add  mass  to  the  driver,  it  should  linearly 
change  the  inverse  of  the  frequency  (or  period)  squared: 

T*  .  i/f*  ■  4if*(a.  ♦  m)/k,  (III-l) 

where  is  the  moving  mass,  a  is  the  added  mass,  and  k  is 

the  suspension  stiffness.  The  instrument  we  used  for  this 
measurement  was  the  HP  4192A  Low  Frequency  Impedance 
Analyzer,  described  in  section  B.  He  added  various  masses  of 
putty  (one  to  nine  grams)  to  the  diaphragm  and  found  the 
resonance  frequency  of  the  driver  for  each  mass.  The  putty 
mass  was  measured  on  the  Sartorius  model  2403  mass  balance 
(serial  number  1287) .  This  balance  has  an  accuracy  of  ± 
0.00005  g.  However,  the  putty  could  be  expected  to  gain  or 
lose  more  mass  than  that  due  to  handling.  Based  on  mass 
measurements  of  the  putty  after  use,  we  expect  an 
experimental  accuracy  of  ±  0.005  g  for  the  masses.  A  more 
significant  uncertainty  involving  the  masses  arises  from  the 
inability  to  place  the  masses  evenly  on  the  voice  coil  or  to 
produce  a  perfectly  rigid  coupling  between  the  masses  and 
the  voice  coil.  It  is  not  possible  to  assign  a  specific 
numerical  percentage  error  to  this  uncertainty,  but  we  must 
take  it  into  account  when  evaluating  the  data  associated 
with  mass  loading.  A  typical  plot  of  the  period  squared  vs. 
added  mass  in  air  is  shown  in  Figure  III-l.  We  used  the 
least  squares  method  to  fit  a  straight  line  through  the 
data.  We  determined  the  suspension  stiffness  (k  =  47t3/slope) 
to  be: 
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(PERIOD  SQUARED  (xlO-6  S*S) 


ALTEC  290-16K  MASS  LOADING 


ADDED  MASS  (g) 


Figure  III-l.  Least  Squares  Fit  of  Period  Squared  vs. 

Added  Mass  to  Determine  Stiffness  and 
Moving  Mass  of  the  Altec  290-16K  Driver. 
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k  »  2.7  ±  0.1  X  10*  N/m,  (III-2) 
and  the  aovinf  aass  vm.  .  intercept /slope) : 

a.  =  5.8  ±  0.1  g,  ( III-3) 
for  this  driver  in  air. 

The  moving  aass  of  the  driver  in  vacuum  is  expected 
to  be  less  than  the  value  in  air  because  in  vacuum  the 
driver  is  not  pushing  a  hydrodynaaic  aass.  However,  we  did 
not  take  any  measurements  of  this  driver  in  vacuum. 

2.  Bellows  pimmtfn 

The  Servometer  custom-designed  dynamic  bellows 
(shown  in  Figure  II I -2  with  test  stand)  is  used  to  provide  a 
flexible  gas  seal  between  the  driver  and  the  resonator. 

The  stiffness  of  the  bellows  was  measured  using  two 
independent  methods  -  static  and  dynamic.  For  the  static 
method  we  used  Hooke's  law:  k  »  -F/x.  The  force  was  deter¬ 
mined  by  pulling  a  partial  vacuum  on  the  bellows  and  multi¬ 
plying  the  vacuum  pressure  (read  on  a  manometer)  by  the 
surface  area  of  the  bellows  (A  =  1.1  +  o.i  X  10-®  ma) .  We 
used  a  G.  L.  Collins  Corporation  model  LMS-119V22-01  Linear 
Variable  Differential  Transformer  (LVDT)  (serial  274674) 
connected  to  an  EG&G  128A  Lock-in  Analyzer  to  determine  the 
displacement  of  the  bellows  under  a  partial  vacuum.  Figure 
III-3  shows  the  LVDT  in  its  calibration  stand.  The  calibra¬ 
tion  curve  for  the  LVDT  is  shown  in  Figure  III-4.  The  least 
squares  fit  slope  of  this  curve  gives  us  the  calibration 
constant  for  the  LVDT: 
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Figure  III -2. *  Photograph  or  Servometer  Custom-Designed 

Dynamic  Bellows — With  and  Without 
Test  Stand. 
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Figure  III-3 .  Photograph  of  Linear  Variable 

Differential  Transformer  ( LVDT ) 
in  its  Calibration  Stand. 
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PERIOD  SQUARED  (xl  0*6  S‘S)  LOCK-IN  OUTPUT  VOLTAGE  (mV  dc) 


LVDT  CALIBRATION 


LVDT  DISPLACEMENT  (INCHES) 

Figure  III-4.  Calibration  Data  for  LVDT  Voltage 

Output  vs.  Displacement . 


BELLOWS  DYNAMIC  MASS  LOADING 


ADDED  MASS  (g) 

Figure  II 1-5.  Least  Squares  Fit  of  Period  Squared 

vs.  Added  Mass  for  Dynamic  Loading 
of  Bellows . 
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dV/dx  «  -0.378  ±  0.000S  V/in  •  -14.9  t  0.02  V/m.  (III-4) 
Table  III-l  gives  the  static  tsst  data  for  ths  stiffness 
measurement.  By  averaging  tha  valuas  wa  obtained: 

k  «  4.3  t  0.2  X  10»  N/s.  ( III-5) 

TABLE  III-l.  BELLOWS  STIFFNESS— STATIC 
MEASUREMENT  RESULTS 


LVDT  Voltage 

Pressure 

Pressure 

Force  («P*a) 

Stiffness 

(mV) 

(cm  Hm) 

(X10*»Pa> 

ilLL. 

4.39 

(X10»N/m> 

14.84 

2.9 

3.87 

4.40 

15.62 

3.1 

4.13 

4.69 

4.47 

12.89 

2.3 

3.07 

3.48 

4.02 

Tha  dynamic  measurement  was  a  mass  loading  technique 
similar  to  that  used  for  the  driver.  In  this  case,  we 
measured  the  period  as  a  function  of  added  mass  using  an 
Endevco  model  2225M2  accelerometer  (serial  FA80)  and  a 
Nicolet  model  3091  digital  storage  oscilloscope  to  record 
the  free  decay  of  the  bellows  after  it  had  been  disturbed  by 
the  blunt  wooden  end  of  a  cotton  swab.  We  then  plotted 
period  squared  versus  added  mass  and  found  the  stiffness  and 
moving  mass  from  the  slope  and  intercept  of  the  least 
squares  fitted  line  (see  Figure  I I 1-5)  as  discussed  in  the 
previous  Section.  This  gave  us  values  of: 
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(III-6) 


k  -  4.58  ±  0.05  X  10®  N/m 
•nd  m_  .  5.84  t  0.05  g.  (III-7) 

This  value  for  Moving  mass  includes  tha  accelerometer  cable 
•nd  is  not  representative  of  tha  actual  moving  mass  of  tha 
hallows.  Subtracting  tha  waight  of  tha  acceleroaete:  givas 
us  a  valua  of  approximately : 

*•  «  1.5  ±  0.1  g.  ( III-8 ) 

for  tha  hallows.  We  know  tha  moving  mass  has  to  be  lass  than 
2.3  g,  since  that  is  tha  bellows'  total  mass — moving  and 
static.  Tha  values  obtained  for  bellows  stiffness  using  the 
two  different  methods  ware  within  experimental  error.  We 
therefore  felt  confidant  using  a  valua  of  ks  ■  4.5  i  0.1  X 
10®  N/a  in  our  later  calculations. 

3.  Combined  System 

After  we  glued  the  reducer  (a  ■  4.8  g) ,  stiffener 
plate  (■  ■  4.4  g)  and  bellows  to  the  driver,  we  wanted  to 
find  the  stiffness  and  moving  mass  of  the  combined  driver- 
reducer-bellows  system.  We  first  tried  a  dynamic  mass 
loading  measurement  (same  procedure  as  with  the  driver 
alone) ,  but  found  that  the  period  squared  versus  added  mass 
plot  (see  Figure  III-6)  gave  unrealistic  results.  So  instead 
we  used  the  relationship: 

F  »  -kx  »  B1I  ( III-9 ) 

or  k  =  -BlI/x  =  B1  (di/dV) (dV/dx)  (III-10) 
and  made  a  static  measurement  using  the  LVDT .  We  expected 
the  additional  equipment  (reducer,  stiffener,  and  bellows) 
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MASS  LOADING  DRIVER-CONE-BELLOWS 


Raductr-Btllom  Systta. 


ALTEC  290-16K  STATIC  DISPLACEMENT 


CURRENT  IN  VOICE  COlL(mAdc) 

Figure  III-7 .  Driver  Displaceaent  vs.  Voice  Coil 

Current  as  Measured  with  LVDT/ 
Lock-in  Analyzer. 
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to  change  the  aoving  mass  and  suspension  stiffness  of  tha 
ayataa  but  not  tha  tranaduction  coafflciant  (Bl).  Wa  appliad 
a  DC  currant  to  tha  drivar  and  racordad  tha  LVDT  voltata 
output,  than  changed  tha  currant  and  racordad  tha  new  vol¬ 
tage.  Thia  procaaa  producad  tha  values:  AI  *  283  *  o.S  aA 
and  gv  ■  2.58  t  0.005  aV.  Wa  aubatituta  thaaa  valuaa  and  tha 
valua  of  Bl  ■  21.5  ±  0.5  N/A  (axpariaantal  aaaauraaant  of 
thia  quantity  la  diacuaaad  in  tha  naxt  section)  into  Kqua- 
tion  III-10  for  a  coabinad  atiffnaaa  valua  of: 

k  «  3.5  ±  0.1  X  10*  N/a.  (III-ll) 
If  wa  add  tha  atiffnaaa  valuaa  for  tha  drivar  and  ballowa 
together  a  total  atiffnaaa  of  k  •  3.2  ^  q  x  x  10*  N/a  for 
tha  ayataa  ia  obtained.  Thua  tha  total  axperlaental  stiff- 
nasa  for  tha  coabinad  ayataa  diffara  by  about  10%  froa  tha 
aua  of  tha  individual  atlffnaaaas.  Using  tha  relationship 
a*  *  k/(2*f)a.  and  a  aaasurad  rasonanca  frequency  of 
235  t  0.5  Hz,  tha  coabinad  aoving  aass  is: 

a.  -  16.2  ±  0.5  g, 

coaparad  to  m*  *  is. 9  ±  0.5  g  for  tha 
added  together.  The  two  types  cf  aass 
within  2%. 

D.  TRANSDUCTION  COEFFICIENT  (BL) 

We  determined  the  Bl  coefficient  of  the  driver  using 
three  independent  methods.  The  first  method  was  to  remove 
the  voice  coil  from  the  magnet  and  directly  measure  the 


(III-12) 
individual  masses 
measurement  agree 
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magnetic  induction  us ins  •  Dowdy  RFL  Industriss  sodsl  912 
Guassaeter  (serial  808)  usins  s  Hsll  offset  probs  sodsl 
912015  (ssrlsl  21122)  (B  «  2.2  ±  o.05  Tssls).  Us  thsn 
ssssursd  1  (ths  lsngth  of  ths  wirs)  by  count ins  ths  nuabsr 
of  turns  of  wirs  in  tbs  coll  (undsr  a  alcroscops)  and 
aultiplyins  this  nuabsr  (44)  by  ths  circusfsrsncs  of  ths 
coil  (1  •  44  X  s  X  dlasstsr  ■  9.98  i  0.005  a).  Coablnins 
thsss  two  asasursssnts  savs  us  ths  valus: 

81  •  22.0  t  0.7  N/A.  ( III-13 ) 
In  ths  sscond  asthod  ws  ussd  ths  rslationshlp : 

F  «  B1I  •  -kx  or  B1  *  -kx/I  «  -k(dx/dV) (dV/dl ) .  (III-14) 
Us  prsviously  dstsrainsd  k  for  ths  drivsr  alons  (Equation 
III-2) .  and  now  used  ths  LVDT  and  lock-in  analyxsr  to  find 
dx/di.  Ths  LVDT  calibration  was  discusssd  in  ths  prsvious 
Section,  with  a  ca leu la tad  valus  of  dV/dx  ■  -14.9±0.05  V/a. 
Us  thsn  asasursd  ths  diaphragm  displacsasnt  as  a  function  of 
DC-current  through  ths  voice  coil.  This  gave  us  a  valus  of 
dV/dl  «  11.75  +  o.OS  aV/A  (sss  Figure  III-7) .  Using  ths 
relation  in  Equation  III-14  ws  obtain  the  value: 

B1  ■  21.2  1  o.l  N/A.  ( III-15) 
In  the  third  method  we  used  the  relationship: 

F  »  Bli  «  mg  or  81  »  ag/i.  (III-16) 
For  this  we  seasured  the  change  in  DC-current  through  the 
coil  as  a  function  of  added  mass,  using  the  LVDT  as  a 
monitor  to  keep  the  displacement  of  the  diaphragm  constant. 
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ThrM  different  aui/currtnt  coabinationa  (m«  Table  III-2) 
produced  an  avar a|a  value  of: 

Bi  •  21.4  *  0.5  N/A.  (III-17) 

Averaging  tha  raaulta  of  these  thraa  aathoda  gava  ua  a 
valua  of  Bl  »  21.5  t  0.5  N/A  with  25  davlatlon  aaong  tha 
thraa  aathoda. 

TABLE  I 11-2.  STATIC  MEASUREMENT  OF  ALTEC  290-16K 
TRANSDUCTION  COEFFICIENT  (BL) 


Al_Ul 

AI  (aA) 

Bl  ■  1/1  (N/A) 

132.5 

64.4 

20.2 

157.5 

74.0 

20.9 

257.5 

108.2 

23.3 

E.  MICROPHONE  CALIBRATION 

Microphona  calibration  ia  an  iaportant  part  of  tho 
process  leading  to  the  neasureaent  of  the  electroacoustic 
efficiency.  A  aicrophone  is  used  to  aeasure  the  acoustic 
pressure  inside  the  resonator  by  converting  the  acoustic 
pressure  to  an  output  voltage.  The  ratio  of  the  open  circuit 
voltage  to  the  pressure  is  called  the  open  circuit 
aicrophone  sensitivity. 

Me  used  three  aicrophones  for  the  calibration:  (1)  a 
General  Radio  (GR)  aodel  1962-9611  half-inch  electret 
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condenser  microphone  (serial  liiil  >  with  OR  preamplifier 
model  1560-P42  (aerial  4061);  (2)  an  Kndevco  model  8510B-5 
piexoreslstive  bridge  microphone  (serial  73NJ);  and  (3)  an 
1/R  inch  BAR  model  4136  microphone  (serial  19663) .  We  first 
calibrated  the  OR  using  a  OR  type  1562- A  sound- level  cali¬ 
brator  (serial  19663,  H4dB  SPL,  250  Hs) ,  and  obtained  a  OR 
sensitivity  of  M*  .  7.42  aV/Pa.  This  compares  within  about 
IX  of  the  calibration  certificate  provided  by  the  manufac¬ 
turer  of  7.50  aV/Pa.  Additional  checks  were  made  with  a  BAK 
plstonphone  calibrator  which  substantiated  the  OR  calibra¬ 
tion  results. 

To  calibrate  the  other  two  microphones  we  attached  the 
resonator  tube  to  the  290-16K  driver.  We  used  an  HP  332SA 
Synthesizer /Function  Generator  to  excite  the  driver  at  240 
Hz,  the  resonant  frequency  of  the  tube.  We  recorded  the 
output  voltages  for  all  three  microphones  and  used  the 
comparison  relationship: 

M»/V*  -  M./V,  (HI-18) 
to  determine  the  sensitivities  of  the  BAK  and  tndevco.  The 

BAK  sensitivity  was  determined  to  be  (b  *  0.463  aV/Pa, 
comparing  within  2X  of  the  manufacturer's  specification 
value  of  0.470  aV/Pa  (adjusted  to  include  the  capacitance  of 
the  BAK  model  VA0160  1/8"  to  1/4”  adapter  with  an  expected  4 
dB  "line  loss").  The  sensitivity  of  the  Endevco  microphone 
was  calculated  to  be  8.17  X  10 v/Pa,  comparing  within  IX 
of  the  manufacturer's  specification  of  8.09  X  10-«  v/Pa.  For 
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all  three  microphones  we  measured  sensitivity  values  within 
2%  of  the  manufacturers'  specifications. 

F.  ELECTROACOUSTIC  EFFICIENCY 

The  electroacoustic  efficiency  is  the  ratio  of  the 
acoustic  power  dissipated  in  the  resonator  to  the  electric 
power  delivered  to  the  driver.  The  electric  power  was 
determined  using  the  relationship: 

=  V  I  cos  e,  ( III-19) 
where  V  is  the  voltage  across  the  driver  terminals ,  I  is  the 
current  flowing  in  the  driver,  and  e  ±s  the  phase  angle 
between  the  voltage  and  the  current.  The  values  for  voltage, 
current  and  phase  angle  delivered  to  the  driver  were 
obtained  using  the  HP  4192A. 

Calculation  of  the  acoustic  power  was  a  little  more 
complex.  We  started  with  the  basic  relationship  for  the 
quality  factor  of  the  tube: 

Q  =  27fEs-r/Ei_e ,  (IIT-20) 
where  E»T  is  the  energy  stored  in  the  tube  and  Ei_c  is  the 
energy  lost  per  cycle.  The  acoustic  power  generated  in  the 
tube  is  the  energy  lost  per  cycle  times  the  number  of  cycles 
per  second,  or  frequency .- 

P<m=  =  Eucf  =  2rtfE®-r/Q.  ( III-21 ) 
Following  Burmaster  (1985)  the  energy  stored  in  a  plane 

standing  sound  wave  is  equal  to  the  energy  density 
’ "tegrated  over  the  volume  of  gas  in  the  tube  (VOL): 
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E.t  *  (Prmm*/2T(P^)  *  VOL  ,  (III-22) 

where  P„m<B  iS  the  root  mean  square  acoustic  pressure  in  the 
tube,  ^  i 3  the  ratio  of  the  specific  heats  of  the  gas  in  the 
tube  (1.4  for  air),  and  P„  is  the  ambient  pressure  (one 
atmosphere  or  1.013  X  10®  Pa  for  these  measurements).  The 
r.m.s.  acoustic  pressure  is  the  r.m.s.  microphone  output 
voltage  divided  its  sensitivity  (V/M)  as  discussed  in  the 
previous  Section.  The  expression  for  the  acoustic  power 
becomes : 

P-c  =  (JffV»VOL)/(QMa#Po>  .  ( II 1-23 ) 

The  only  remaining  quantity  in  Equation  II 1-23  that 

needs  to  be  calculated  separately  from  the  efficiency 
measurement  is  the  quality  factor  of  the  tube  (Q) .  To 
determine  Q  we  removed  the  290- 16K  driver  from  the  resonator 
and  replaced  it  with  an  electret  transducer  to  drive  the 
tube.  This  was  necessary  because  there  were  losses  due  to 
the  290- 16K  that  made  the  tube  Q  appear  to  be  lower  than  it 
actually  is  (we  discovered  this  problem  while  investgating 
the  reasons  for  the  calculated  efficiency  to  be  greater  than 
one).  We  used  the  B&K  microphone  for  the  pressure  (converted 
to  voltage)  readings  at  the  other  end  of  the  tube.  We  plot¬ 
ted  voltage  as  a  function  of  frequency,  and  determined  the 
resonant  frequency  (f0)  and  frequencies  of  the  half  power 
points  (f^  and  f»).  Using  the  relationship: 

0  =  f«/(f„-fi)  ( III-24 ) 

we  obtained  a  value  of  0  =  60  +  l  for  the  tube. 
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To  determine  the  electroacoustic  efficiency  we  returned 
the  tube  to  the  290-  16K  driver  (see  Figure  III-8)  and  drove 
the  tube  at  resonance.  We  recorded  the  voltage,  current,  and 
phase  angle  for  the  driver  and  the  B&K  output  voltage.  The 
electroacoustic  efficiency  was  then  calculated  using: 

=  P*c/P«x«m=.  ( III-25) 
Table  III-3  lists  the  electric  power,  acoustic  power,  and 
efficiency  for  five  test  runs  at  driver  currents  ranging 
from  3-38  mA  (with  corresponding  voltages  0.15-1.43  V).  The 
electroacoustic  efficiency  ranged  from  84.2  for  the  lowest 
driver  current  to  75.9  for  the  highest. 

TABLE  III-3 .  ALTEC  290-16K  ELECTROACOUSTIC 
EFFICIENCY  MEASUREMENTS 


(mW> 

P~=  (mW) 

(X) 

0.56 

0.47 

84.2 

3.41 

2.70 

79.2 

10.0 

7.79 

77.9 

17.0 

13.1 

77.2 

54.2 

41.2 

75.9 

G.  COMPARISON  TO  COMPUTER  MODEL 

The  final  task  in  the  evaluation  procedure  for  the 
electrodynamic  driver  is  to  compare  the  measured  electro¬ 
acoustic  efficiency  to  that  predicted  by  the  computer  model 
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Figure  III-8 .  Photograph  of  Altec  290-16K  Driver 

Assembly  with  Resonator  Tube  and 
Endevco  Microphone. 
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discussed  in  Section  II.D.  We  have  already  measured  or 
calculated  all  the  input  values  for  the  model  except  for  the 
driver's  mechanical  resistance  <Rm). 


l.  Mechanical  Resistance 

We  can  determine  Rm  by  using  the  relationship: 

R«  a  2rtf «*m0/Q ,  ( III-26 ) 

where  f0  is  the  resonant  frequency,  su  Is  the  moving  mass 
and  Q  is  the  quality  factor  of  the  driver.  There  are  several 
methods  available  for  determining  the  driver’s  quality 
factor.  One  is  to  plot  the  magnitude  of  the  modified  driver 
impedance  vs.  frequency  and  pick  the  half  power  points  off 
the  graph.  Then,  as  with  the  tube,  Q  is  given  by: 

Q  =  fo/<fM-f*),  ( HI-27 ) 

where  f„  and  fi  are  the  frequencies  associated  with  the 
upper  and  lower  half  power  points  (located  at  [Zl^^/VZ, 
where  [Z]  is  the  magnitude  of  the  driver  impedance  which  is 
modified  by  subtracting  the  direct  current  (DC)  resistance 
from  the  measured  complex  impadance) . 

Another  method  is  to  use  the  relationship  [Hofler, 


1986]  : 


O  *  J-  f.  d$/df , 


( III-28 ) 


where  0  is  the  phase  angle  between  driver  resistance  (R)  and 
reactance  (X),  the  two  components  of  modified  impedance  (Z  = 
R  +  iX,  where  R  is  the  measured  resistance  minus  the  DC 
(direct  current)  resistance).  This  me. ified  impedance  is 
proportional  to  the  voice  coil  velocity.  The  quantity  d^/df 


75 


is  the  slop*  of  the  phase  angle  versus  frequency  curve 
evaluated  at  the  resonant  frequency  (where  the  curve  is  a 
straight  line) .  We  used  the  HP4192A  to  obtain  the  data  for 

these  measurements. 

Figure  III-9  shows  a  plot  of  modified  impedance  vs. 
frequency  for  the  Altec  290-16K  in  air.  From  this  graph  we 

find  f..235.2  Hz,  f«»249.6  Hz  and  f» -222.0  Hz.  The  Q  for 
these  values  is  8.S2.  Figure  I XX- 10  shows  a  plot  of  phase 
angle  vs.  frequency  for  the  same  data  set.  Over  this  large 
range  of  frequencies  the  phase  angle  curve  approximates  an 
arctangent.  At  resonance,  however,  where  d$/df  is  evaluated, 
the  curve  is  linear  (see  Figure  III-ll).  Equation  III-28  is 
the  slope  of  the  arctangent  at  resonance  (f«).  The  slope  of 
the  plotted  data  is  71.8X10-®  sec  (dd  was  converted  from 
degrees  to  radians) .  The  Q  calculated  from  this  method  is 
8.44.  These  values  for  Q  agree  within  1%.  From  these  values 
we  get  R^-2.77  kg/s. 

We  tried  yet  another  method  for  determining  driver 
Q.  This  third  method  measures  the  decay  constant,  %,  using  a 
digital  storage  oscilloscope  to  record  the  "ring  down"  of 
the  driver.  This  technique  takes  less  time  to  run  the  exper¬ 
iment  and  reduce  the  data  using  the  relationship: 

Q  *  ®t/2  or  x  =  2Q/m  =  QT/*  ( III-29) 
where  ®  is  the  resonant  angular  frequency  in  radians/sec,  T 
is  the  period  in  seconds,  and  t  is  the  free  decay  time 
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ALTEC  IMPEDENCE  PHASE  ANGLE 


Figure  III-10.  Altec  290-16K  Modified  Impedance 

Phase  Angle  vs.  Frequency — Full 
Frequency  Range. 


ALTEC  PHASE  ANGLE  -  LINEAR  SEGMENT 


Figure  III-ll.  Least  Squares  Fit  of  Altec  290-16K 

Modified  Ispedance  Phase  Angle  as  a 
Function  of  Frequency  Near  Resonance. 
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constant  in  seconds.  This  sathod  will  be  discussed  in 
greater  detail  in  Chapter  IV. 

2.  Bfidftl  Rtmltl 

Table  II 1-4  lists  the  coaputer  model  input  values 
for  the  Altec  290-16K.  Figures  III-12  through  III-15  show 
the  output  plots  for  these  values  as  discussed  in  Section 
II. D.  In  this  case  the  resonance  peaks  are  much  closer 
together,  and  the  magnitude  of  tha  secondary  peak  is  closer 
in  value  to  the  primary  peak.  Once  again  the  driver  reso¬ 
nance  is  at  a  lower  frequency  than  the  resonator  resonance. 
It  is  interesting  to  note  that  the  modeled  efficiency  (about 
60ft)  is  15-20ft  less  than  the  measured  efficiency  reported  in 
the  previous  Section.  We  believe  that  this  results  from  the 
fact  that  for  the  coaputer  model,  a  high  tube  0  means  that 
the  driver  losses,  represented  by  R„,  play  a  more  important 
role  in  the  efficiency  calculation,  and  the  results  are  very 
sensitive  to  small  errors  in  R*,.  The  primary  ways  for  the 
driver  energy  to  be  dissipated  are:  (1)  Joule  heating  due  to 
DC  resistance  in  voice  coil  and  eddy  currents  induced  in  the 
magnet  structure  due  to  the  time- varying  magnetic  fields 
generated  by  the  AC  voice  coil  current  (small  but  not  negli¬ 
gible)  ;  (2)  mechanical  resistance  of  driver  (R„);  and  (3) 

energy  transferred  to  acoustic  load.  Another  possible  reason 
for  the  different  efficiencies  is  that  the  efficiency  is  not 
independent  of  current  hence  the  real  physical  sys  .am  has 
some  non-linearity,  as  Table  III-3  clearly  shows.  The 
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TABLE 

111*4.  COMPUTER  MODEL  INPUT- 

-ALTEC 

290-16K  DRIVER 

gyalrel 

Mamin* 

unit* 

RM 

Mechanical  resistance 

kg/s 

2.77 

M 

Moving  easa 

kg 

0.0162 

.SUP 

Spring  constant 

N/a 

3 . 52X10* 

RHO 

Density  of  gas 

kg/e® 

1.2 

C 

Sound  speed 

a/s 

342.0 

S 

Resonator  X-sec .  area 

a* 

9. 35X10-* 

BL 

Driver  trans .  factor 

N/A 

21.5 

Q 

Resonator  quality  factor 

60.0 

FO 

Resonator  resonance  freq. 

HZ 

240.0 

A 

Back  voluee  X-sec.  area 

a* 

4. 13X10-® 

VOL 

Back  voluee 

a® 

2.8X10-* 

LE 

Voice  coil  inductance 

H 

1.0X10-® 

RE 

Voice  coil  DC  resistance 

ft 

10.6 

I 

Driver  current 

A 

1.0 

L 

Resonator  length 

a 

0.712 
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VOLTAGE  VS.  FREQUENCY 


MTiTTTTTTl 

(gffO  juqoua  mum 


Figure  III- 13.  Computer  Model  Output  for 

Altec  290-16K — Piston  Velocity 
vs.  Frequency  for  One  Amp. 
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Figure  III -14.  Computer  Model  Output  for 

Altec  290-16K — Acoustic  Power 
vs.  Frequency  for  One  Amp. 


ELECTROACOUSTIC  EFFICIENCY  VS. 


computer  model  la  baaed  on  a  llnaar  theory  and  will  have  an 
efficiency  which  la  Independent  of  drive  current. 

H.  DRIVER  RESONANCE  ALTERATION 

The  refrigerator  must  operate  at  the  resonant  frequency 
of  the  reaonator  to  obtain  the  hifheat  electroacoustic  effi¬ 
ciency  and  tc  aalntain  stack  position  with  respect  to  the 
standing  acoustic  wave  (discussed  in  Sub-section  II. A. 4).  It 
appears  from  the  computer  model  results  and  our  experiments 
so  far  that  the  operating  frequency  should*  also  be  near  the 
driver’s  resonance  to  get  the  most  power  out  of  the  driver 
The  resonant  frequency  of  the  resonator  is  determined  by  the 
effective  length  of  the  resonator  and  the  speed  of  sound  in 
the  contained  gas.  The  effective  length  for  a  tube  of  uni¬ 
form  cross-section  and  closed  end  is  the  actual  measured 
length.  For  the  STAR  resonator,  however,  the  effective 
length  is  the  length  that  would  produce  the  measured  reso¬ 
nant  frequency  for  the  given  sound  speed  if  the  resonator 
were  a  uniform  cross-section  closed-enied  tube  (see  Figure 
1-2  for  approximate  shape  of  resonator).  The  sound  speed  can 
be  altered  by  changing  the  percentage  of  xenon  in  the 
helium-xenon  mixture,  since  the  sound  speed  depends  on  the 
density  of  the  mixture,  and  the  density  depends  on  the 
percentages  of  the  individual  gases  present  (see  Equations 
IV-24  and  IV-26) .  The  gas  mixture  and  effective  resonator 
length  will  be  chosen  to  maximize  the  coefficient  of 


performance  of  the  refrigerator.  Therefore  if  we  want  to 
natch  the  resonances <wewi 11  have  to  alter  the  resonance  of 
the  driver  once  the  resonant  frequency  of  the  resonator  has 
been  determined. 

The  resonant  frequency  of  the  driver  is  proportional  to 
the  square  root  of  the  stiffness  divided  by  the  moving  mass. 
If  we  want  to  raise  the  resonant  frequency  of  the  driver  we 
can  either  increase  the  stiffness  or  decrease  the  mass. 
Convaraely,  to  lower  the  resonant  frequency  we  can  either 
decrease  the  stiffness  or  increase  the  mass.  Generally  we 
prefer  to  decrease  stiffness  ar.d  aass  to  get  a  larger 
velocity.  Reduction  of  stiffness  usually  also  reduces  the 
mechanical  resistance. 

For  the  Altec  290-10K  attached  to  the  uniform  •  cross- 
section  resonant  tube  in  these  preliminary  experiments,  the 
frequencies  were  already  well  matched  (note  the  closeness  of 
the  poaks  in  Figures  III-12  through  III-14)  and  changing  the 
ariver  frequency  was  not  required.  However,  we  know  that  in 
future  applications  we  say  need  to  alter,  most  probably 
lower,  the  driver  resonant  frequency.  We  developed  an 
experimental  procedure  to  do  this. 

Adding  moving  mass  to  the  driver  would  appear  to  be  the 
easiest  way  to  lower  the  driver  resonance.  Unfortunately  the 
added  mass  would  also  reduce  the  velocity  by  increasing  the 
mecahnical  impedance  above  resonance.  We  therefore  chose  to 
experiment  with  decreasing  the  driver  stiffness.  Recall  that 


thm  driver -reducer-bellows  tyitM  stiffness  includes  the 
stifffness  of  the  bellows  (4.5  X  10®  n/b)  and  the  suspension 
stiffness  of  the  driver  phenolic  surround  (29  X  10®  n/b). 
The  driver  suspension  was  the  obvious  choice  for 
experimental  alteration  for  two  reasons:  (1)  it  is  over  six 
tlaes  as  larfe  as  the  bellows  stiffness,  and  (2)  it  would  be 
very  difficult  for  us  to  change  the  bellows  stiffness 
without  affecting  its  perforaance. 

Since  the  290- 16K  voice  coil  surround  is  made  of  a 
phenolic  material,  it  was  relatively  easy  to  cut  through  it 
with  an  ordinary  scalpel.  We  made  a  series  of  successive 
cuts  in  the  surround  material  and  ran  measurements  of  the 
driver  resonant  frequency  after  each  cut.  The  results  of 
this  experiment  are  given  in  Table  III-5.  We  plotted  the 
phase  angle  between  the  driver's  electromotive  force  (emf) 
and  the  voice  coil  displacement  as  a  function  of  frequency. 
For  the  emf  we  measured  the  voltage  across  the  driver  termi¬ 
nals  in  channel  A  of  the  HP  4192A  Impedance  Analyzer  and 
assumed  that  the  emf  was  in  phase  with  the  driver  current 
and  hence  force.  We  used  the  photonic  sensor  to  measure  the 
displacement  by  connecting  its  voltage  output  to  channel  B. 
The  4192A  then  displayed  the  phase  angle  between  the  two 
voltages  representing  force  and  displacement.  The  driver  is 
at  resonance  when  this  phase  angle  equals  90  degrees .  The 
quality  factor  (Q)  is  equal  to  one-half  the  resonant  fre¬ 
quency  times  the  slope  of  the  phase  angle  vs.  frequency 
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curve  evaluated  at  resonance  (see  Equation  HI-28) .  Using 
the  previously  determined  valuas  tor  moving  mass  (16  g)  and 
ballovs  stiffness,  we  calculstad  tha  suspansion  stiffnass: 

*****  -  k%-%  -  ktolUM  (III-30) 
uhora  k%«%  *  Wa  calculstad  tha  machanical 
rasistanca  (.Km)  using  aquation  111-26.  Wa  stoppad  cutting 

tha  surround  whan  its  stiffnass  was  approxiaataly  aqual  to 
that  of  tha  ballows  and  tharafora  tha  total  stiffnass  was  no 
longar  suspansion  doainatad. 

TABLE  III -5.  RESULTS  OF 
DRIVER  RESONANT  FREQUENCY  ALTERATION 


Percent  Surround 

Resonant 

Quality 

Suspension 

Mechanical 

Remaining 

Frequency 

Factor 

Stiffness 

Resistance 

.(.eppeml  ?. 

(H&L- 

LKg/gl 

100 

241.0 

8.55 

32.2 

2.83 

100  (six  slits) 

229.5 

8.09 

28.8 

2.85 

80 

213.5 

5.17 

24.3 

4.15 

50 

183.5 

8.65 

16.8 

2.13 

20 

119.5 

8.56 

4.52 

1.40 

Despite  the 

anoaoly  at 

the  80% 

data  point, 

we  feel  that 

this  experiment 

shows  that 

partial 

removal  of 

the  surround 

can  reduce  stiffness,  resonant  frequency,  and  mechanical 
resistance.  It  appears  from  the  above  data  that  the  flexing 


Xf  _  -i  - 


of  tho  phono lie  ourround  is  rosponmiblo  for  tho  majority  of 
tho  mocha ni cal  loomom  in  tho  driVor. 


89 


XV.  STAR  DRIVIR  PARAMETIRS 

A.  INTRODUCTION 

This  chapter  discusses  the  selection  of  the  electro- 
dynamic  driver ,  describes  the  parameter  measurements  on  the 
chosen  driver,  end  estimates  the  electroacoustic  efficiency 
of  the  STAR  using  the  computer  sodel  described  in  Section 
II.  D. 

l.  interface  with  NASA  QAS  Cannister 

In  addition  to  the  aaxiaua  weight  of  200  pounds 
discussed  in  Section  I.B,  restrictions  due  to  the  Get  Away 
Special  (GAS)  cannister  dimensions  give  us  a  aaxiaua  payload 
height  (28.25  Inches)  and  diaaeter  (19.75  inches).  In  our 
allowed  height  we  must  fit  two  battery  layers,  the  driver 
with  its  housing,  the  resonator  assembly,  and  its  associated 
vacuus  space  and  theraal  insulation.  The  battery  layers  are 
each  about  four  inches  high,  for  a  total  battery  height  of 
eight  inches.  This  leaves  about  20  inches  for  the  driver  and 
resonator  asseablies  (see  Figure  1-4  for  approxiaate 
configuration) .  The  electronics  package  and  gas  reservoir 
should  fit  in  the  area  radially  surrounding  the  resonator, 
and  will  therefore  not  add  to  the  height.  Although  the 
diaaeter  of  the  cannister  voluae  is  19.75  inches,  the  driver 
housing  aust  attach  to  the  standard  GAS  can  lid  at  a  9.5 
inch  diaaeter  bolt  circle,  otherwise  we  would  need  a  lid 


specially  designed  for  our  prelect.  As  far  as  weight  is 
conosrned .  approximately  every  pound  we  save  elsewhere 
allows  us  to  carry  an  extra  battery  cell ,  which  results  in 
an  extra  ten  watt-hours  of  available  energy. 

2.  Maximizing  Kleetroacoustic  IfUfillDSac, 

in  addition  to  the  space  and  weight  restrictions,  we 
want  to  maximize  the  electroacoustic  efficiency  of  the 
driver.  He  can  do  this  by  obtaining  the  aaximus  value  for 
the  cosbination  of  parameters  Bl/R^m*.  A  larger  B1  will 
give  us  a  larger  driving  force  (P*Bli),  which  in  turn  gives 
us  a  larger  acoustic  vcluse  velocity,  one  of  the  components 
of  acoustic  power.  For  a  given  value  of  magnetic  induction 
(B),  a  longer  voice  coil  (achieved  by  having  a  larger  voice 
coil  diameter  with  the  same  number  of  turns,  or  the  same 
diameter  and  more  turns)  will  result  in  a  larger  Bl.  This 
longer  1  would  also  result  in  a  greater  voice  coil 
electrical  resistance,  so  the  two  effects  cancel.  We  also 
want  to  minimize  the  moving  mass  (a.)  to  allow  greater 
velocities  and  minimize  the  driver's  mechanical  resistance 
(R»)  in  order  to  minimize  mechanical  losses  in  the  driver. 


3.  Pclvflr  gfllfls&laa 

We  had  two  off-the-shelf  drivers  available  for  con¬ 
sideration.  The  Altec  260-16K  was  discussed  in  Chapter  III. 
The  other  driver  was  a  JBL  2445J.  Both  drivers  were  designed 
for  applications  somewhat  different  from  the  STAR.  Figure 
IV-i  shows  a  side  view  of  the  244SJ  driver.  The  hatched 
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Figure  iv-i.  Schematic  Drawing  of  JBL  244SJ  Driver. 

Hatched  Area  Shows  Unnecessary  Throat 
Section. 
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area,  known  aa  the  throat  of  the  drivar,  la  conaidarad 
unnecessary  for  our  application,  and  adda  a  lot  of  axtra 
waifht  to  tha  drivar.  Tha  Altac  290-16K  haa  a  aiailar  axcasa 
throat  volume.  Ma  contactad  JBL  and  aakad  if  they  would 
ouatoa-daaicn  a  drivar  which  would  aliainata  thia  axcaaa 
weight.  Portunately  for  our  projact,  JBL  waa  in  tha  procaaa 
of  designing  a  naw  drivar,  tha  2450 J,  which  uaaa  neodymium- 
iron-boron  aagnata  inataad  of  tha  traditional  caraaic.  Thay 
agreed  to  custom-design  a  vergion  of  tha  2450 J  for  our  uaa. 
Tha  raaultlng  drivar  (hereafter  r afar rad  to  aa  NIB)  la  ahown 
in  Figure  IV-2.  Tha  NIB  drivar  accapta  tha  aaaa  voice  coils 
(four  inch  diameter,  aluainua  wire,  and  titanium  suspen- 
sion/diaphraga)  aa  tha  2445 J  (ceramic  magnet)  drivar.  JBL 
claims  that  it  should  also  have  tha  same  properties , 
including  Bl. 

In  tha  remainder  of  thia  chapter  we  describe  the 
parameter  teats  on  the  NIB  driver.  We  then  compare  tha  NIB 
driver  to  the  290- 16K  driver.  The  results  of  the  2445 J 
parameter  measurements  are  given  for  comparison.  Finally  we 
present  results  generated  by  the  computer  model  for  the  NIB 
driver  given  various  acoustic  load  impedances. 


B.  NIB  DRIVER  PARAMETERS 

For  the  measurement  of  the  NIB  driver  parameters,  we 
used  the  techniques  we  developed  in  Chapter  III  for  the 
measurement  of  the  290-16K  driver  parameters . 
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Figure  IV-2.  Photographs  of  JBL  Neodynium- 
Iron-Boron  Driver. 


i.  tttvlM  ttaiii  *n&  SuBaanilsn  stiff nm 

The  ,  moving  mass  amt  suspension  stiffness  wara 

ealoulatad  using  tha  mass  loadinc  technique.  Wa  placed  aaall 
aaaaaa  (0.3  to  3.3  grams)  of  black  putty  on  top  of  tha 
volca  coil  in  half-graa  increments .  for  aach  added  aaaa  tha 
raaonant  frequency  of  tha  drivar  was  found  using  tha  HP 
4192A  LP  lapadanca  Analyser  (saa  Saction  IXX.B  for  a 
dascription  of  this  instruaant).  Figura  IV-3  ar.owa  a  aaapla 
plot  of  period  squarad  vs.  addad  aaaa  for  tha  NIB  drivar  in 
air.  Tha  stiff naas  was  found  by  calculating  four  pi  squarad 
dividad  by  tha  slopa.  Tha  aoving  aaaa  is  aqual  to  tha 
intarcapt  dividad  by  tha  slopa.  For  tha  NIB  in  air  thasa 


valuas  ara: 

k  *  9.8  ±  0.1  X  10*  kg/s*  (IV-1) 

and  a«  «  7.8  ±  0.1  g,  (IV-2) 

with  a  resonant  frequency  of  564  Hz.  Wa  repeatad  tha 
aeasureaent  in  a  vacuus  chamber  and  found: 

k  -  8.6  *  o.l  X  10*  kg/s*  (IV-3) 

and  au  «  3.2  ±  0.1  g.  (IV-4) 


with  a  resonant  frequency  of  833  Hz  (see  Figure  IV-4). 

The  aoving  mass  is  greater  in  air  because  of  the 
driver  having  to  move  tha  air  itself.  Tha  discrepancy  in  the 
stiffness  is  not  as  easy  to  explain.  To  determine  the  actual 
value  of  k,  we  used  the  static  mass  loading  method.  This 
gave  us  a  value  of: 


k  -  6.4  ±  0.2  X  10*  kg/s*.  ( IV-5) 

(The  static  loading  measurements  will  ba  discussed  in  tha 
next  Sub-section . ) 

This  last  value  for  suspension  stiffness  differs  by 
about  40%  from  the  first  two.  We  did  not  have  sufficient 
time  to  test  this  driver  until  we  were  sure  we  had  the 
"correct”  value  for  the  stiffness.  We  recommend  that  the 
follow-on  students  do  so.  Lacking  the  correct  value,  we  ran 
the  computer  model  for  two  sets  of  data:  one  set  represen¬ 
ting  the  "high"  k  (9.0  X  10*  kg/s*)  and  one  using  the  "low" 
k  (6.4  X  10*  kg/s*) . 

2.  Transduction  Coefficient  (Bl) 

The  transduction  coefficient  (Bl)  was  calculated 
using  two  methods.  The  first  method  was  to  measure  B  (magne¬ 
tic  induction  in  the  voice  coil  gap)  and  1  (length  of  the 
voice  coil)  directly.  We  measured  B  with  the  Dowty  RFL 
Industries  model  912  Guassmeter  (serial  803)  using  a  Hall 
effect  probe  model  912015  (serial  21122) .  The  maximum 
reading  using  this  probe  was  18.7  kGauss  =  1.87  Tesla  (T) . 
This  value  differs  by  less  than  one  percent  from  the  value 
stamped  on  the  driver  by  the  manufacturer:  1.858  T.  We  will 
therfore  use  a  value  of  B=1.86  ±  0.01  T  for  this  driver.  We 
then  counted  the  turns  of  aluminum  wire  on  the  voice  coil  (n. 
=  30  ±  l  turns)  and  calculated  1  using  the  relationship: 

1  =  n  d  n.  ( IV-6 ) 
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For  the  9.98  *  0.005  cm  diameter  coil  we  calculated  a  value 
of  1  *  9.4  ±  o.3m.  Multiplying  these  two  values  together 
gives  us  for  the  transduction  coefficient: 

B1  -  17.5  ±  0.6  N/A.  ( IV-7) 

The  other  method  used  to  determine  B1  used  the 
relationship: 

F  »  B1I  *  -kx  or  B1  ■  -kx/I,  (IV-8) 

where  I  is  the  direct  current  applied  to  the  coil,  k  is  the 
spring  constant  measured  earlier,  and  x  is  the  displacement 
of  the  voice  coil  produced  by  the  current.  For  this  measure¬ 
ment  on  the  290- 16K  driver  we  had  used  a  Linear  Variable 
Differential  Transformer  (LVDT) .  For  the  NIB  driver, 
however,  we  used  Mechanical  Technology  Incorporated's  MTI- 
1000  photonic  sensor  (serial  38010818)  with  plug-in 
cartridge  model  MTI-3812  (serial  1366)  and  probe  model  125R. 
Figure  IV-5  shows  the  photonic  sensor  with  the  NIB  driver. 
One  major  advantage  of  the  photonic  sensor  over  the  LVDT  is 
that  the  photonic  does  not  need  to  make  physical  contact 
with  the  surface  it's  measuring  (it  uses  light)  whereas  the 
LVDT  does.  With  the  Altec  290-16K's  rigid  phenolic  dia¬ 
phragm,  the  contact  with  the  LVDT  probably  did  not  affect 
the  measurement.  The  NIB  diaphragm,  however,  is  made  of  thin 
-  titanium  and  is  quite  flexible.  It  would  have  been  impos¬ 
sible  to  get  an  accurate  reading  of  the  displacement  of  this 
diaphragm  using  the  LVDT.  We  calibrated  the  photonic  sensor 
using  calibration  fixture  KD-CH-III  A.  The  output  voltage 
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from  the  sensor  was  read  on  a  HP  3478A  multimeter  (serial 
A110136) .  We  then  plotted  DC  voltage  vs.  displacement  in 
millimeters.  The  resulting  graph  is  shown  in  Figure  IV-6. 
The  slope  of  this  line  gives  us  the  calibration  constant: 

dV/dx  *  -2.086  ±  0.001  V/mm.  (IV-9) 

This  value  is  within  two  percent  of  the  calibration  value 
provided  by  the  manufacturer  of  dV/dx  »  2.045  V/mm  (back 
slope) . 

The  next  step  was  to  run  DC  current  through  the  coil 
and  record  the  displacement  of  the  diaphragm  (via  a  voltage 
reading  from  the  photonic  sensor).  For  the  voltage  supply  we 
used  the  HP  467A  power  amplifier  (serial  444-00243) .  We 
found  the  current  by  measuring  the  voltage  (with  the  HP 
3476A  digital  multimeter,  serial  292162)  across  a  ten  ohm 

*108)  resistor  placed  in  series  with  the  circuit,  and 
dividing  the  voltage  by  ten  ohms.  The  resulting  graph  is 
shown  in  Figure  IV-7.  The  slope  of  the  line  gives  us: 

dl/dV  «  2.10  ±  0.005  A/V.  (IV-10) 

Combining  the  values  in  Equations  IV-9  and  IV-10  with  the 
previously  calculated  high  end  value  for  the  suspension 
stiffness  k=9.2±0.6X  10^  kg/s3  gives: 

B1  =  k(dx/dV) (dV/dl )  =  21.5  ±  1.4  N/A.  (IV-11) 
This  15%  difference  between  the  values  of  B1  using 
the  two  independent  methods  is  what  prompted  us  to  test  for 
B1  and  k  using  the  static  loading  method .  Our  primary 
concern  is  that  we  would  expect  the  direct  measurements  of  B 
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CURRENT  (A)  DC  VOLTAGE  (V) 


•  IV-6 .  Calibration  Data  for  Photonic  Sensor 
Voltage  Output  vs.  Displacement. 


NIB  -  PHOTONIC  CURRENT  CALIBRATION 


VOLTAGE  (V) 

Figure  IV-7.  Least  Squares  Fit  of  DC  Current  Through 
Voice  Coil  vs.  Photonic  Voltage  Output. 
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and  1  to  give  us  the  highest  possible  value  for  Bl.  This  is 
because  the  maximum  aagnatic  field  nay  not  be  acting  on  the 
full  length  of  the  coll.  We  did  not  expect  the  value  of  Bl 
calculated  using  the  spring  constant  to  be  1S7,  higher  than 
the  direct  measurement  value. 

We  were  able  to  obtain  static  loading  values  for 
both  Bl  and  k  using  the  same  equipment  setup  and  the 
formulas : 

Bl  *  Amg/AI  and  k  ■  -  Amg/Ax  (IV-12) 
For  these  measurements  the  voice  coil  is  driven  with  a  DC 
current.  We  then  add  previously  measured  masses  (putty  and 
soldering  wire)  to  the  voice  coil,  and  measure  the  resulting 
displacement  of  the  diaphragm  using  the  photonic  sensor .  To 
measure  Bl  we  determine  how  much  current  is  needed  to  bring 
the  voice  coil  back  to  its  null  displacement  (original 
output  voltage  reading  fron  the  photonic  sensor)  after  the 
mass  is  added.  We  measure  k  by  keeping  the  current  constant 
and  measuring  the  difference  in  photonic  output  voltage 
(multiplied  by  the  calibration  constant  to  get  displacement) 
resulting  from  the  added  mass.  Table  IV-l  lists  the  results 
of  these  measurements. 

The  bracketed  value  in  each  k  and  Bl  pair  was 
calculated  using  the  measured  value  in  the  pair  and  the 
relationship: 

B1I  =  -kx.  ( IV-13 ) 
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From  these  Maiurtaanti  we  get  the  value  k  •  6.^  ±0.2  kg/s* 
as  reported  aarliar,  and  tha  transduction  coefficient: 

Bl  •  15.0  ±  0.2  N/A.  (IV-14) 
For  tha  nodal  runs  we  usad  a  valua  of  Bl  *  19.0  for  tha  high 
k  and  Bl  *  1S.0  for  tha  low  k. 


TABLE  IV-1.  STATIC  LOADING  OF  NIB  DRIVER  TO  DETERMINE 
SUSPENSION  STIFFNESS  AND  TRANSDUCTION  COEFFICIENT 


An  (s) 

AI  <SA) 

ay  <yl 

k  (X10*  ks/s*) 

Bl  (N/A) 

17.1 

5.2 

6.7 

[15.3] 

17.1 

— 

5.5 

6.4 

[14.6] 

34.7 

23.3 

— 

C6.4] 

14.6 

34.7 

22.3 

C6.7] 

15.3 

34.7 

— 

.  288 

6.1 

— 

(front  slope) 

3. 

Resist  ftn<?9 

The  next 

parameter  we 

measured  was  the 

mechanical 

resistance  (R*,)  using  the  relationship: 

R-  -  2*f.n./Q.  (IV-15) 
Me  attempted  to  find  the  driver  quality  factor  (Q)  using  the 
same  modified  impedance  technique  described  in  Section 
III.G.  As  can  be  seen  in  Figure  IV-8,  the  graph  of  the 
modified  impedance  magnitude  vs.  frequency  is  not  a  smooth 
resonance  curve,  as  it  was  for  the  290- 16K  (see  Figure  III- 
9).  In  fact,  it  looks  like  the  superposition  of  two 
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different  impedancw  curves ,  the  first  on*  shorter  and  fatter 
and  the  second  one  taller  and  thinnar.  It  was  iaposslbla  to 
accurately  aea&ure  half -power  points  froa  this  curve,  so  we 
chose  the  phase  angle  aethod.  (It  is  interesting  to  note 
that  the  JBL  2445 J  (ceraaic)  driver  did  have  a  smooth  modi¬ 
fied  impedance  curve  similar  to  that  of  the  290-16K. )  Figure 
XV-9  shows  the  graph  of  phase  angle  vs.  frequency  over  a 
wide  range  of  frequencies.  Although  this  curve  generally 
approximates  an  arctangent,  the  segaent  of  the  line  immedi¬ 
ately  adjacent  to  the  resonant  frequency  (aaauaed  to  be  the 
frequency  where  d-0,  in  thia  case  534.5  ±  0.5  Ha)  is  linear 
(for  a  frequency  range  of  about  ten  Hz).  This  segment  is 
plotted  in  Figure  IV-10.  The  slope  of  the  line  is  dd/df  »  - 

42.0  ±  0.5  X  10-®s  (after  converting  degrees  to  radians).  As 

in  Section  XI.G,  Q  is  given  by: 

Q  «  - ( f o/2 ) ( dd/df )  ■  11.2  ±  0.2.  (IV-16) 

This  gives  a  value  for  R„  (using  the  moving  mass  in  air 
«  7.a  ±  0.1  g) : 

R«.  «  2.34  ±  0.05  kg/s.  ( IV-17) 

We  decided  that  a  more  accurate  measurement  of  Q  and 

R«,  could  be  found  by  making  the  G  measurement  in  vacuum.  For 
this  we  used  the  open  circuit  free  decay  aethod.  This  aethod 
consists  of  sending  a  signal  to  the  driver  at  resonance  with 
an  HP  3314A  function  generator  (serial  A342044) ,  then 
simultaneously  switching  off  the  signal  and  triggering  a 
Nicolet  model  310  storage  oscilloscope  (serial  88004641)  to 
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PHASE  ANGLE  (DEG) 


NIB  IMPEDANCE  PHASE  ANGLE 


480  500  520  540  560  580 

FREQUENCY  (HZ) 

Figure  XV-tf.  NIB  Modified  Impedance  Phase  Angle 
a*  a  Function  of  Frequency — Full 
Frequency  Range. 


NIB  PHASE  ANGLE  -  LINEAR  SEGMENT 


520  530  540  550 


FREQUENCY  (HZ) 

Figure  iv-io.  NIB  Modified  Iapedance  Phase  Angle 

as  a  Function  of  Frequency — Linear 
Segaent  Near  Resonance. 
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record  the  free  decay.  A  aaeple  of  the  output  fros  this 
process  is  shown  In  Figure  iv-li.  The  equation  for  a  slsple 
harmonic  oscillator  experiencing  free  decay  Is: 

A(t)  •  A«  exp  <-t/*)  sin(»*t+#),  (IV-18) 

where  A(t)  Is  the  amplitude  as  a  function  of  time,  A«  is  the 
Initial  (maximum)  amplitude,  t  is  the  time  (*  nT,  or  the 
number  of  periods  times  the  period ) ,  x  is  the  free  decay 
constant  (■  GT /*,  see  Equation  111-29) ,  is  the  angular 
resonant  frequency,  and  #  is  the  phase  angle.  We  recorded 
the  peak  amplitude  as  a  function  of  period  number,  and 
plotted  the  natural  logarithm  of  the  peak  amplitude  vs. 
period  number.  The  equation  representing  this  line  is: 

In  A ( T ) -  In  A.  -  nT/t.  (IV-19) 

\ 

The  slope  of  this  line  is  t/x  *  nT/*t  *  nx/Q.  Therefore  Q  can 

be  found  by  dividing  *  by  the  slope  of  the  line.  Figure  IV- 

12  shows  a  semi-log  plot  for  the  NIB  driver  in  air.  The 
value  of  Q  for  this  plot  is  21.5,  which  is  almost  twice  the 
value  we  found  using  the  modified  impedance  phase  angle 
measurements  (Q  *  11.2).  Obviously  the  air  strongly  affects 
the  measurements.  The  scattering  of  the  data  points  in  this 
figure  is  probably  due  to  the  same  phenomena  that  causes  the 

double  peaks  in  the  impedance  curve.  Figure  IV-13  shows  the 

semi-log  plot  of  the  amplitude  vs.  period  number  for  vacuum. 
This  curve  is  much  smoother  and,  we  assume,  more  reliable. 
The  Q  for  this  data  is : 

Q  =  X/slope  ■-=  19.8  ±  0.05.  ( IV-20 ) 
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Figure  IV-n.  Sample  of  Free  Decay  Output  from 

Nlcolet  310  Storage  Oscilloscope. 


4 


HIP  FREE  DECAY- AIR 

y  •  &S792  *  0.1464X  R-1.00 


Figure  IV-12.  Semi-log  Plot  of  Amplitude 


vs.  Period  Number  for  NIB 
Free  Decey  in  Air. 


Period  • 

Fifure  IV-13.  Semi-log  Plot  of  Aeplitude 

vs.  Period  Number  for  NIB 
Free  Decay  in  Vacuum. 
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«e  can  also  find  tha  period  froa  the  free  decay  curve  (and 
thus  the  reaonant  frequency  aa  the  lnverae  of  the  period). 
He  uae  the  oedlloacope  to  eeaeure  the  tine  for  ten  cycles 
and  divide  this  nuaber  by  ten.  For  the  vacuus  run  the 
resonant  frequency  was  839  ±  0.5  Hs.  These  values  (and  the 
aovinf  sass  *  3.2  f  in  vacuus)  give  us  the  sechanical  resis¬ 
tance  in  vacuus: 

EU  -  2*f«s«/0  a  0.83  i  0.01  kf/s.  (IV-21) 
This  is  the  value  we  used  for  the  aodel. 

Table  IV- 2  fives  a  cospar ison  of  applicable 
paraaeters  for  the  three  drivers.  Notice  that  the  paraaeter 
Bl/R„„— is  alaost  three  tiaea  greater  for  the  NIB  driver 
than  the  290-16K.  This  is  the  paraaeter  that  we  said  (Sub¬ 
section  IV. A. 2)  we  wanted  to  saxiaise  to  get  the  highest 
theraoacouatic  efficiency.  Thus  the  driver  that  gives  us  the 
best  height  and  weight  characteristics  (see  photograph  of 
the  thi  drivers  in  Figure  IV-14)  also  gives  us  the  best 
electroacoustic  efficiency . 

C.  MODELED  EFFICIENCY 

In  order  to  aodel  the  electroacoustic  efficiency  of  the 
driver,  we  need  to  know  what  acoustic  load  (ispedance)  it 
will  see.  The  final  resonator  and  stack  coaposition  will  not 
have  been  deterained  by  the  tiae  this  thesis  is  coapleted. 
He  therefore  assuaed  a  range  of  possible  values  for  the 
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TABLE  IV-2. 

COMPARISON  OF  PARAMETERS  FOR 
THE  THREE  DRIVERS 

Property  (units) 

Altec  290 

JBL  2445J 

IBkHIB 

Moving  mass,  m.  (kg) 

0.0151 

0.015 

0.015 

Spring  Constant  (N/m)  3.52X10“ 

•*6.5X10“ 

•*9.0X10“ 

Magnetic  Induction, 

P  (Tesla) 

2.1 

1.97 

1.87 

Voice  Coil 

Diameter  ( cm ) 

7.22 

9.98 

9.98 

B1  (N/A) 

21.5 

16.5 

17.5 

Voice  coil  DC 
resistance  (a) 

10.6 

7.1 

7.6 

Mechanical 

resistance  (kg/s) 

2.77 

•*1.0 

0.83 

Voice  coil 

inductance  (H) 

(at  20  kHz) 

1.0X10-“ 

1.14X10"“ 

1.14X10- 

Mass  (kg) 

14.5 

13.8 

3.1 

Height  (in) 

6.0 

5.5 

2.5 

Magnetic  energy 
density  (J/cm3) 

**0.5 

•*0.5 

3.2 

Bl/R _ ;hn0 

(m/Coul-kg) 

510 

1230 

1410 

Ill 


Figure  IV-14.  Photographic  Comparison  of  the 

Three  Drivers — Altec  290- 16K, 
JBL  2445 J ,  and  JBL  NIB. 


acoustic  load  based  on  the  experimental  data  acquired  to 
date.  J 

The  model  parameters  that  relate  to  the  acoustic 
load  are  the  density  of  the  gae  (RHO) ,  the  speed  of  sound  in 
the  gas  (C) ,  the  resonator  quality  factor  (Q),  the  resonator 
resonant  frequency  (FO),  and  the  resonator  length  (L).  The 
model  was  designed  tor  a  closed -ended  tude  of  uniform  cross- 
section.  The  STAR  resonator  has  an  open  tube  (with  a  bulb 

attached)  and  nonuniform  cross-section.  The  problem  of 

>  ■ 

applying  the  existing  computer  model  to  the  STAR  resonator 
was  solved  by  scaling  the  STAR  parameters  to  make  it  look 
like  a  closed-ended  uniform  cross-section  tube, 
l.  Incut  Values  for  Pure  Helium 

For  pure  helium,  the  experimental  data  available  was 
the  range  of  resonator  quality  factor  values  (10-20),  and 
the  speed  of  sound  and  gas  density  at  room  temperature 
(about  300K)  and  ten  atmospheres  (c„T  =  1024  m/s, 

*  1.618  kg/m*) .  Also  available  were  four  resonant  fre¬ 
quencies,  two  for  a  short  resonator  at  warm  (580  Hz)  and 
cold  (520  Hz)  temperatures,  and  two  for  a  long  resonator  at 
warm  (520  Hz)  and  cold  (470  Hz)  temperatures.  Using  this 
information,  we  first  calculated  the  effective  lengths  of 
the  resonators  at  room  temperature  as  if  they  were  closed- 
ended,  uniform  cross-section  tubes.  For  this  we  used  the 
relationship: 

*  -  2L  «  c/f,  or  L.*.*  =  c/2f ,  (IV-22) 
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where  \  is  the  wavelength,  L  is  the  length  of  the  tube,  c  is 
the  speed  of  sound,  and  f  is  the  resonant  frequency.  This 
calculation  gave  us  effective  lengths  for  the  resonators: 
L.  a  0.883m  and  Li  »  0.985m.  For  the  computer  runs  we  use  a 
convenient  value  of  *  0.90m.  We  then  calculated  the 
speed  of  sound  for  the  operating  temperature  of  approxi¬ 
mately  250K  using  the  relationship: 

C(T)  -  C«r  tT/T^]1'*.  (IV-23) 
This  gives  us  a  speed  of  sound  cme  «  935  m/s.  If  we  use 
this  sound  speed  and  the  effective  resonator  lengths  to 
estimate  the  resonant  frequencies  for  250K,  we  get  f,  *  529 
Hz  and  =  474  Hz,  which  agree  with  the  measured  values.  We 
calculated  the  gas  density  at  the  operating  temperature 
using  the  relationship: 

P  »  <P/C»  «1.94  kg/m®,  ( IV-24) 
where  <  is  the  ratio  of  the  specific  heats  (1.667  for 
helium)  and  P  is  the  pressure. 

Table  IV-3  lists  the  input  parameters  for  the  pure 
helium  computer  runs.  Note  that  for  the  computer  model,  the 
moving  mass  must  take  the  mass  of  the  reducer  cone  and 
bellows  into  account  (see  Figure  IV-15).  The  total  moving 
ma3S  is  expected  to  be  approximately  15  g. 

2.  Input  Values  fo£  Helium-Xenon 

There  was  somewhat  less  data  available  for  the 
helium-xenon  mixture  (12.5%  Xe) .  We  had  the  sound  speed  at 
300K  and  ten  atmospheres  (c„x  =  472  m/s).  We  also  had 
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TABLE  IV-3 .  COMPUTER  MODEL  INPUT- 
NIB  DRIVER  WITH  HELIUM  AT  TEN  ATMOSPHERES 
AND  TEMPERATURE  *  250K 


;-V.  ; 

Symbol 

MsuUni 

Units 

NIB 

RM 

Mechanical  resistance 

kg/s 

0 . 83 

M 

Moving  sass 

kg 

0 . 015 

SSUP 

Spring  constant 

N/m 

6.6,  9.0  X10** 

_ 

RHO 

Density  of  gas 

kg/m® 

1.94 

C 

Sound  speed 

m/s 

935.0 

S 

Resonator  X-sec .  area 

a* 

1. 14X10“® 

BL 

Driver  transduction  factor 

N/A 

15.0,  19.0 

Q 

Resonator  quality  factor 

10.15,20 

FO 

Resonator  resonance  freq. 

Hz 

519.4 

♦ 

A 

Back  volume  X-sec.  area 

na 

4. 42X10“® 

•* 

VOL 

Back  volume 

m® 

3.28X10“** 

LE 

Voice  coll  inductance 

H 

1.14X10“** 

RE 

Voice  coil  DC  resistance 

ft 

7.6 

I 

Driver  current 

A 

1.0 

L 

Resonator  length 

m 

0.90 
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rai onant  frequencies  for  the  saae  four  conditions  as  for  ths 
purs  hslius:  short-wars  (265  Hz),  short-cold  (230  Hz),  long- 
wars  (235  Hz),  and  long-cold  (215  Hz).  Us  ussd  Equation 
XV-22  to  find  ths  sffsctivs  rssonator  lsngths  ( Lm  .  0.891a 
and  Lx  »  1.004s).  Howsvsr,  ones  again  ws  ussd  a  convsnisnt 
sffsctivs  lsngth  of  0.90  a  for  ths  cosputsr  runs.  Us  ussd 
Equation  IV-23  to  find  ths  sound  spssd  for  Hs-Xs  at  the 
operating  teaperature,  cOTO  .  431  a/s.  Ths  estiaatsd  reso¬ 
nant  frequencies  for  th«  cold  teaperature  were  f_*242  Hz  and 
f i*215  Hz,  ones  again  coaparing  favorably  with  the  aeasured 
values.  To  find  the  density  of  the  aixture  at  the  operating 
teaperature  we  used  a  fora  of  ths  ideal  gas  law: 

e  *  PM/RT,  ( IV-25) 
where  M  is  the  atoaic  aass  of  the  aixture.  R  is  the 
universal  gas  constant,  and  P  and  T  are  the  pressure  and 
temperature.  Ue  calculated  the  atomic  mass  of  the  mixture 
using  the  equation: 

=  Mh«Xh«  +  ,  (IV-26) 
where  the  M’s  are  atomic  masses  (4  for  helium  and  131.3  for 

xenon)  and  the  x’s  are  percentages  of  the  gases  in  the 
mixtures.  This  equation  gives  us  an  atomic  aass  for  the 
12. 5X  Xe  mixture  of  19.9  g/aol,  and  a  density  of  9.70  kg/m® 
at  250K  and  ten  atmospheres.  Table  IV-4  lists  the  computer 
input  values  for  He-Xe,  including  a  range  of  Q's  from  5  to 
30.  (We  are  not  sure  what  to  expect  for  the  actual  value  of 
Q,  since  a  new  stack  will  also  be  used.) 


TABLE  IV-4 .  COMPUTER  MODEL  INPUT-- 
NIB  DRIVER  WITH  HELIUM-XENON  AT  TEN  ATMOSPHERES 
AND  TEMPERATURE  *  250K 


Symbol 

tteaajjaa 

ynUs 

NIB 

RM 

Mechanical  resistance 

kg/s 

0.83 

M 

Moving  aass 

kg 

0.015 

SSUP 

Spring  constant 

N/a 

6.6,  9.0  X10* 

RHO 

Density  of  gas 

kg/a® 

9.70 

C 

Sound  speed 

a/s 

431.0 

S 

Resonator  X-sec .  area 

a* 

1. 14X10“® 

BL 

Driver  trans .  factor 

N/A 

15.0,  19.0 

Q 

Resonator  quality  factor 

5,10,15,20,25,30 

FO 

Resonator  resonance  freq. 

Hz 

239.4 

A 

Back  volume  X-sec.  area 

m* 

4. 42X10“® 

V'L 

Back  volume 

a® 

3.28X10“* 

LE 

Voice  coil  inductance 

H 

1.14X1C’* 

!2 

Voice  coil  DC  resistance 

ft 

7.6 

Driver  current 

A 

1.0 

L 

Resonator  length 

a 

0.90 

3-  Computer  Model  Results 

Figures  iv-16  through  ZV-19  ere  sample  output  plots 
for  pure  heliua  with  81  ■  15  N/A,  k  *  6.4  X10*  kg/s*.  end 
Q  ■  20.  Notice  thet  the  peek  velues  for  the  driver  end 
resonator  ere  elsost  matched  (the  driver  resonance  Is  at  the 
lower  frequency) .  The  thermoecoustic  efficiency  for  this  set 
of  velues  Is  about  70. A%.  Figures  IV-20  through  IV-23  show 
the  corresponding  plots  for  the  helium-xenon  mixture.  This 
time  the  peeks  ere  not  matched  (the  driver  peek  is  off  the 
scale)  end  the  electroacoustic  efficiency  is  about  16. 6%. 
The  primary  reason  for  the  difference  between  the  two  gases 
is  the  lower  sound  speed,  and  therefore  lower  resonant  fre¬ 
quency,  in  the  He-Xe  mixture. 

Table  IV-5  lists  the  maximum  value  of 
electroacoustic  efficiency  for  each  of  the  test  runs.  The 
consistently  low  efficiencies  for  the  He-Xe  emphasize  the 
importance  of  matching  the  two  resonant  peaks.  He  recommend 
that  the  driver  frequency  be  adjusted  to  match  the  resonator 
operating  frequency.  This  can  best  be  done  by  reducing  the 
stiffness  of  the  titanium  surround  (see  Section  III.H). 


VOLTAGE  VS. 


Figure  IV- 16.  Computer  Model  Output  for  NIB  with 

Helium  (B1  «  15,  Q  *  20) — Driver 
Voltage  vs.  Frequency  for  One  Amp. 
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PISTON  VELOCITY  VS  FREQUENCY 


Figure  IV-17.  Computer  Model  Output  for  NIB  with 

Helium  (B1  »  15,  O  »  20) — Piston 
Velocity  vs.  Frequency  for  One  Amp. 
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ACOUSTIC  POWER  VS. 


Figure  IV -18.  Computer  Model  Output  for  NIB  with 

Helium  (B1  «  15,  0  «  20) --Acoustic 
Power  vs.  Frequency  for  One  Amp. 
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ELECTROACOUSTIC 


Figure  IV- 19.  Computer  Model  Output  for  NIB  with  Helium 

<B1  «  15,  O  *  20) — Electroacoustic 
Efficiency  vs.  Frequency. 
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VOLTAGE  VS.  FREQUENCY 


Figure  IV-20.  Computer  Model  Output  for  NIB  with  Helium 

Xenon  (B1  »  15,  Q  «  20) — Driver  Voltage 
vs.  Frequency  for  One  Amp. 
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PISTON  VELOCITY  VS.  FREQUENCY 


Figure  IV-21 .  Computer  Model  Output  for  NIB  with  Helium- 

Xenon  (B1  =  15,  Q  »  20) — Piston  Velocity 
vs.  Frequency  for  One  Amp. 
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ACOUSTIC  POWER  VS.  FR 


Figure  IV-22.  Computer  Model  Output  for  NIB  with  Helium- 

Xenon  (B1  =  15,  Q  *  20) --Acoustic  Power 
vs.  Frequency  for  One  Amp. 
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ELECTROACOUSTIC  EFFICIENCY  VS.  FREQUENCY 


Pi*ure  IV-23 .  Computer  Model  Output  for  NIB  with  Helium- 

Xenon  (B1  »  15,  O  *  20) — Electroacoustic 
Efficiency  vs.  Frequency. 
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TABLE?  IV-S.  COMPUTER  MODEL  RESULTS  FOR  NIB  DRIVER- 
PEAK  ELECTROACOUSTIC  EFFICIENCY  (%) 


a 

H«-X«  (Bl»15) 

H«-Xe  (Bl-19) 

s 

5.1 

5.9 

10 

71.3 

70.4 

9.2 

10.7 

15 

71.6 

74.3 

13.1 

15.1 

20 

70.4 

76.2 

16.6 

19.1 

25 

19.6 

22.7 

30 

22.8 

26.0 

v.  SlAfiL  ORIVSR  gQNSTRVSTIQN 


A.  INTRODUCTtON 

This  chapter  describes  the  construction  of  the  STAR 
driver.  A  large  part  of  this  thesis  involved  the  design  of 
the  parts  necessary  to  integrate  the  NIB  driver  with  the 
resonator  and  GAS  cannister,  as  illustrated  in  Figure  1-4. 

The  eaior  design  work  was  for  the  driver  housing.  The 
requirements  for  the  driver  housing  are: 

1  -  lightweight, 

2  -  strong  (pressuri2ed  to  210  psi), 

3  -  interface  to  GAS  canister  lid  (12  10-32  bolt  holes 

centered  at  9 . 5  inches  diameter ) , 

4  -  interface  with  resonator  (preferably  interchangeable 

with  earlier  test  resonators) , 

5  -  helium  leak-tight  feed-through  ports  for  electrical 

leads  (driver  electric  power,  accelerometer  and 

microphone  outputs ) ,  and 

6  -  hold  NIB  magnet  structure  securely. 

Figure  V-l  shows  a  cross-sectional  drawing  of  the  driver 
housing  (not  including  any  of  the  feed-through  ports,  which 
will  be  discussed  later).  Distinction  between  parts  is  shown 
by  different  hatching.  Part  A  is  the  NIB  magnet  structure 
which  was  discussed  in  Chapter  IV.  Part  B  is  the  main  body 
of  the  driver  housing,  machined  out  of  6061 -T6  aluminum. 
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Figure  V-i. 


Cross-sectional  Schematic  of  the  STAR  Driver 
Housing  Showing  (A)  NIB  Driver,  (B)  Driver 
Housing  Main  Body,  (C)  Aluminum  Pusher  Plate, 
(D)  Pressure  Lid,  (E)  Leak  Test  Lid, 

(F)  Resonator  Vacuum  Housing,  and  (G)  Bolt 
Holes  for  Interface  with  Resonator. 


Part  C  la  a  6061-T6  aluminum  pusher  plate  used  to  hold  the 
aatnet  structure  securely  in  the  body  using  six  10/32  bolts. 
Part  0  is  the  pressure  lid,  mealed  a  1/I6th  inch  diameter 
lead  o-ring  and  held  in  place  with  six  10/32  bolts.  Part  E 
represents  part  of  the  leak  test  lid  (see  Figure  C-l).  When 
completed  this  will  be  replaced  by  the  GAS  can  lid  provided 
by  NASA.  The  driver  housing  must  match  up  with  the  GAS  can 
lid's  9.5  inch  diameter  10-32  bolt  circle.  Part  F  is  the 
resonator  section  housing.  The  arrows  labeled  G  point  to  the 
bolt  holes  which  fora  the  driver-resonator  interface. 

B .  DESIGN  CONSIDERATIONS 

The  driver-resonator  system  will  be  pressurized  to  ten 
atmospheres  with  a  helium-xenon  mixture  (12.5%  xenon). 
NASA's  regulations  say  that  it  must  be  pressure  tested  to 
1.5  times  the  operating  pressure  (NASA,  1986a),  or  15  atmo¬ 
spheres  (*1. 52X10*  pa  «  220  psi) .  An  added  complication  of 
working  with  helium  (or  helium  mixtures)  is  that  it  diffuses 
through  almost  all  materials  except  metals.  As  a  result,  we 
are  not  able  to  use  rubber  o-rings  to  seal  metal  parts 
together,  but  must  use  lead  o-rings  instead. 

l-  Electrical  Feed-Throughs 

NASA's  regulations  prohibit  us  from  leaking  gases 
into  space,  primarily  to  avoid  contamination  of  other  exper¬ 
iments  and  for  the  safety  of  the  Shuttle  and  astronauts. 
This  regulation  is  important  not  only  for  the  major  o-ring 
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•••la,  but  also  for  the  Microphone  and  accelerometer  outputs 
to  tha  frequency  control  alactronics  (saa  Figure  1-5),  and 
connacting  tha  audio  aaplifiar  to  the  driver.  There  are 
three  individual  wires  and  three  co-axial  connectors  that 
must  penetrate  tha  driver  housing  without  allowing  leakage. 
The  other  housing  penetrations  are  a  DC  pressure  gage  for 
the  driver  housing  and  a  fill  line  from  a  helium-xenon 
bottle  which  will  stabilize  the  pressure  in  the  system 
against  temperature  changes  in  the  environment  and  those 
generated  by  the  refrigerator. 

Based  on  prior  experience  in  other  applications,  we 
decided  that  Emerson  &  Cuming’s  STYCAST<tw>  2850FT  epoxy 
(Canton,  MA  02021)  would  be  the  best  choice  for  providing 
the  electrically  insulating  seal  for  the  wires  in  the 
housing.  STYCAST«tm>  2850PT  is  a  quartz-filled  epoxy  formu¬ 
lation  with  remarkably  good  over-all  general  properties.  In 
addition  to  having  excellent  electrical  grade  insulation 
properties,  good  adhesion,  and  excellent  resistance  to  chem¬ 
icals  and  solvents,  2850FT  has  unusually  high  thermal  con¬ 
ductivity  and  low  thermal  expansion.  This  combination  of 
properties  makes  this  epoxy  valuable  for  applications  where 
electrical  insulat^n  and  mechanical  protection  must  be 
maintained  ur.ile  ,  jping  with  heat  transfer  considerations. 

One  definition  for  a  sealed  container  is  a  helium 
leak  rate  of  less  th  •  10-**  standard  cc's  per  second  (NASA, 
1986b) .  To  test  epoxy  we  designed  and  built  a  test 
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fixture  to  be  used  with  the  Veeco  MS-9  helium  leak  detector. 
A  cross-sectional  diagram  of  the  test  fixture  is  shown  in 
Figure  v-2.  Part  A  is  a  copper  tube  with  1-1/8"  outer  diam¬ 
eter  to  interface  with  the  leak  detector.  Soldered  to  the 
copper  tube  is  a  p.ece  of  solid  brass  stock  about  one  inch 
thick  (B)  with  a  0.052  inch  diameter  hole  bored  through  the 
center.  Epoxied  into  the  hole  is  a  0.025  inch  diameter  bare 
copper  wire  (C).  This  design  is  motivated  by  the  desire  to 
provide  as  long  a  travel  path  and  as  small  a  cross-section 
as  possible  for  the  helium  gas  molecules.  Part  D  is  a  hollow 
aluminum  housing  designed  to  hold  pressurized  helium.  We 
decided  that  a  rubber  o-ring  should  hold  the  helium  suffi¬ 
ciently  well  for  a  short  term  leak  test.  Part  E  is  a  1/8" 
NPT  (National  Pipe  Thread)  hole  for  connecting  to  a  helium 
fill  line.  A  photograph  of  the  test  fixture  with  and  without 
the  housing  is  shown  in  Figure  V-3. 

The  copper  tube  and  brass  stock  were  soldered  toge¬ 
ther.  We  boiled  the  assembly  to  remove  excess  particulates 
and  rinsed  it  with  alcohol.  We  then  etched  the  hole  in  the 
brass  and  the  copper  wire  with  nitric  acid  diluted  with 
distilled  water  (about  20%  acid)  to  provide  a  clean,  rough 
surface  for  the  epoxy.  Before  etching,  cellophane  tape  was 
placed  on  the  portion  of  the  brass  stock  where  the  o-ring 
would  be  sealing  (after  making  sure  that  the  tape  would  not 
dissolve  in  the  acid).  We  then  rinsed  the  assembly  with 
distilled  water  and  dried  it,  mixed  the  epoxy,  and  inserted 
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Figure  V-2 .  Cross-sectional  Schematic  of  the  Helium 
Leak  Detector  Test  Fixture  Showing 

(A)  Copper  Tube  for  Leak  Detector  Interface. 

(B)  Solid  Brass  Stock.  (C)  Bare  Copper  Wire. 
(D)  Aluminum  Housing,  and  (E)  Connection  for 
Helium  Fill  Line. 


Figure  V-3 .  Photograph  of  the  Helium  Leak  Detector 
Test  Fixture. 
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the  wire  into  ths  hole  in  ths  brass  stock,  making  surs  ths 
•poxy  fillsd  ths  tap. 

Ths  following  procsdurs  was  used  to  run  ths  leak 
dstsctor  tssts:  first  zero  ths  sstsr,  than  do  a  tost  run  on 
ths  ssnsitivity  calibrator.  Ths  ssnsitivity  calibrator  ussd 
was  a  Vsscs-typs  sc-4  serial  number  3179  with  a  leak  rate  of 
6.5X10***  std  cc/ssc  (standard  cubic  centisstsrs  per  second). 
We  would  then  reaove  the  calibrator,  rechsck  the  aster  for 
sero,  attach  ths  test  fixture,  read  the  aeter,  reaove  the 
fixture,  and  recheck  for  zero. 

The  teat  fixture  was  initially  pressurized  to  115 
psi  with  heliua.  We  tried  to  test  it  inmediately  but  found 
that  the  detector  would  not  draw  enough  of  a  vacuus  to  check 
for  heliua.  We  thought  that  there  Bight  be  outgaasing  from 
flux  or  soaething  in  the  tube.  So  we  pulled  a  vacuua  on  the 
tube  overnight  and  the  next  day  it  worked  fine.  We  tested 
the  fixture  on  the  first,  third,  and  sixth  days  after  it  had 
been  pressurized.  The  numerical  results  are  listed  in  Table 
V  -i .  Even  after  we  had  waited  six  days  for  the  heliua  to 
diffuse  through  the  epoxy,  the  fixture  leakage  rate  was  an 
order  of  aagnitude  less  than  the  calibrator.  The  leakage 
rate  was  therefore  four  orders  of  aagnitude  less  than  that 
allowed  by  NASA  to  be  considered  leak-tight  ( 1X10— *  std 
cc/sec) .  After  the  third  leak  test  we  checked  the  pressure 
in  the  aluminum  chamber  and  found  it  to  be  within  a  few 
percent  of  the  original  pressure.  We  therefore  concluded 
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that  the  bar*  wira-apoxy  aathod  with  a  length/diameter  ratio 
of  about  20  could  ba  usad  for  tha  drivar  housing  panatrating 

alactrical  f eed- through* . 

TABLE  V-i.  HELIUM  LEAK  DETECTOR  TEST  RESULTS 
ON  COPPER  WIRE-EPOXY  FIXTURE. 


Number  days  after 

null  aster 

calibrator 

copper  wira  feed 

pressurization 

raadins 

ruding 

-through  ending 

1 

0 

l.l 

0.02 

3 

0 

1.2 

0.08 

6 

0 

1.1 

0.10 

Figura  V-4  shows  a  closa-up  of  tha  feed-through  plug 
for  tha  wira  leads.  It  will  housa  a  total  of  five  wires 
(with  space  available  for  a  sixth  if  necessary ) .  Two  of  tha 
leads  will  ba  for  the  voltage  supply  to  the  voice  coil. 
These  leads  will  consist  of  bare  copper  wires  epoxied  into 
tha  snail  holes  seen  in  tha  figura.  Each  wire  will  then  be 
soldered  to  a  6/32  phillip-head  screw  with  a  1-72  blind- 
tapped  hole  in  it.  The  screw  will  ba  glued  into  an  3/32 
hole,  forming  a  mechanically  stable,  demountable,  electrical 
terminal.  The  accelerometer  and  microphone  leads  will  also 
have  bare  copper  wire  epoxied  into  the  small  hole,  and  will 
then  be  connected  to  Microdot «▼«>  coaxial  fittings  (Malco 
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Figure  V-4.  Driver  Housing  Feed-Through  Plug 
for  the  Electrical  Leads. 
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part  number  031-0001-0001,  South  Passadena,  CA  91030). 
Figure  V-5  shows  an  example  of  each  of  these  leads. 

2.  Measurement  Devices 

The  microphone  and  accelerometer  used  by  STAR  are 
similar  to  those  used  by  Hofler  (1986  and  1988).  In  fact, 
the  resonator  end  of  the  driver  housing  was  designed  to 
accept  Hofler 's  resonator,  calibration  volume,  and  instru¬ 
mentation  for  testing.  The  microphone  used  to  measure  the 
acoustic  pressure  in  the  resonator  is  a  small  quartz-crystal 
(Valpey-Fisher ,  Hopkinton,  MA  01748)  or  PZT  dynamic  pressure 
transducer  (Piezoelectric  Products  lot  number  4-9-84, 
F013184 ,  Metuchen,  NJ) .  It  is  mounted  in  the  driver  housing 
where  the  housing  connects  to  the  resonator  flange.  The 
microphone  will  be  calibrated  by  comparison  to  a  standard 
transducer,  which  will  have  been  previously  calibrated  via  a 
reciprocity  technique.  The  microphone  will  be  connected  to  a 
miniature  FET  follower  circuit.  The  FET  buffer  (Eletec  model 
304  impedance  converter,  Daytona  Beach,  FL  32020)  converts 
the  microphone’s  high  impedance  output  signal  to  a  low 
impedance  to  reduce  sensitivity  to  electrical  noise  and 
maintain  the  output  voltage  and  phase  while  driving  the 
cables  necessary  to  reach  the  signal  processing  electronics. 
The  buffer's  DC  power  supply  and  output  cables  are  complete¬ 
ly  shielded  through  the  driver  housing  to  the  buffer  (see 
Figures  C-3  and  C-7  for  location  of  microphone  and  FET  in 
driver  housing) . 


(b) 

Figure  V-5.  Example  of  Wire  Lead  Connectors 

for  Feed-Through  Plug:  (a)  Single 
Wire  and  (b)  Co-axial  Cable. 


140 


Tha  Endevco  nodal  22  Picomin<TM>  accelerometer  is 
attached  to  tha  side  of  the  bellows  stiffener  facing  the 

driver.  The  integrated  accelerometer  signal  will  be  used  as 
a  measure  of  the  bellows  velocity,  which  will  be  calibrated 
by  driving  a  small  cavity  of  calculaole  acoustic  impedance. 
Thus  the  product  of  the  outputs  from  the  microphone  and 
integrated  acceleroaeter  can  be  combined  to  yield  the  acous¬ 
tic  power  and  their  relative  phase  can  be  used  to  keep  the 
driver  frequency  at  the  tube  resonance.  The  operation  of  the 
refrigerator  at  resonance  is  essential  both  to  optimize  the 
acoustic  power  delivered  to  the  "stack"  and  to  guarantee 
that  the  "stack"  is  always  located  at  the  proper  position  in 
the  standing  wave  as  discussed  in  Sub-section  II. A. 4. 

The  DC  pressure  transducer  is  an  Omega  PX-80 

designed  for  high  pressure  applications  requiring  long-term 
stability.  Figure  V-6  shows  a  drawing  of  the  pressure  trans¬ 
ducer  and  its  driver  housing  fitting. 

The  helium  port  in  the  housing  is  provided  to 
empty  and  fill  the  driver  and  resonator  and  so  that  a  stable 
pressure  may  be  maintained  in  the  driver  and  resonator 
through  the  ballast  volume  provided  by  the  reservoir  (see 
Figure  1-4). 

There  is  a  1/32  inch  diameter  hole  filled  with  a  ten 
thousandths  inch  inner  diameter  half  inch  long  capillary  in 
the  driver  housing  where  it  connects  to  the  resonator.  The 
function  of  the  capillary  is  to  equilibriate  the  static 


141 


pressure  between  the  driver  and  resonator  during  filling  and 
non-operating  tines  while  maintaining  a  pressure  seal  at 
acoustic  operating  frequencies . 

A  complete  set  of  dimensioned  drawings  for  the 
driver  housing  and  its  auxiliary  equipment  is  given  in 
Appendix  C.  A  list  of  materials  for  the  driver  portion  of 
the  STAR  is  given  in  Appendix  D. 
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VI.  CONCLUSION  AND  RECOMMENDATIONS 


The  objective  of  the  STAR  project  is  to  test  and  space 
qualify  a  new  continuous  cycle  cryogenic  refrigeration 
systes  for  cooling  of  sensors  and  electronics  which  is  based 
upon  the  newly-discovered  thersoacoustic  heat  pusping 
effect.  The  new  refrigerator  has  no  sliding  seals,  a  cycle 
frequency  of  about  300  Hz,  and  uses  acoustic  resonance  to 
enhance  overall  efficiency.  This  thesis  is  concerned  speci¬ 
fically  with  the  design  and  testing  of  the  electrodynanic 
transducer  which  is  responsible  for  the  electro-acoustic 
power  conversion.  A  computer  model  of  the  driver /resonator 
system  was  presented  along  with  the  techniques  for  measure¬ 
ment  of  the  electrical  and  mechanical  parameters  used  as 
input  for  the  model. 

A  final  driver  design  (including  dimensional  drawings) 
utilizing  a  modified  JBL  2450J  neodymium-iron-boron 
compression  driver  and  associated  leak-tight  electrical 
feed-throughs,  microphone,  accelerometer,  pressure  gage, 
pressure  housing,  and  resonator  interface  is  provided. 

Since  we  only  received  the  NIB  driver  in  late  December 
1987,  there  wasn't  time  to  perform  enough  parameter 
measurements  to  obtain  accurate  values  for  the  suspension 
stiffness  (k)  and  transduction  coefficient  (Bl).  The  next 
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group  of  students  should  find  these  values  for  use  in  the 
■odel  end  efficiency  calculations. 

The  resonator  experisents  (not  reported  here)  indicate 
that  Helium-Xenon  may  be  the  preferred  cab  mixture  for  the 
STAR.  If  this  is  the  case,  the  NIB  driver  resonant  frequency 
will  have  to  be  lowered  for  it  to  match  the  operatinc 
(resonant)  frequency  of  the  resonator.  This  could  be 
accomplished  by  adding  movinc  mass  to  the  voice  coil- 
reducer-bellows  assembly ,  but  that  would  decrease  the 
velocity  and  lower  the  electroacoustic  efficiency.  The 
preferred  option  for  driver  resonance  adjustment  is  to 
decrease  the  suspension  stiffness  by  cuttinc  away  part  of 
the  surround  (see  Section  III.H).  Since  the  NIB  voice  coil 
has  a  titanium  surround  and  dlaphracm  this  operation  will 
have  to  be  performed  carefully.  The  titnnium  diaphragm 
(dome)  will  also  have  to  be  removed.  (The  dome’s  function  is 
performed  by  the  reducer  and  bellows.) 

The  task  for  the  next  croup  of  students  is  to  assemble 
the  NIB  driver  with  its  housinc  and  accessories  ( feed¬ 
throughs  ,  microphone,  accelerometer,  etc.)  and  perform 
efficiency  measurements  on  the  cround  similar  to  those  that 
will  be  performed  in  space.  This  series  of  measurements  will 
closely  resemble  those  performed  on  the  Altec  290- 16K  driver 
as  described  in  Chapter  III. 
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COMPUTER  MODEL  CALCULATION  PBQSRAM 


PROGRAM  FRIDGE 

C  INTRODUCTION: 

C 

C  THE  PURPOSE  OF  THIS  PROGRAM  WRITTEN  BY  LT.  MICHAEL  P. 
C  SUSALLA  IS  TO  CALCULATE  THE  DATA  POINTS  AND  PLOT  THE 
C  PERFORMANCE  CHARACTERISTICS  OF  A  CONCEPTUAL  MODEL  OF  A 
C  THERMOACOUSTIC  HEAT  ENGINE.  THE  COMMENTS  INCORPORATED  INTO 
C  THIS  PROGRAM  AND  THE  COMMENTS  INCORPORATED  IN  THE  PLOTTING 
C  PROGRAM  "GRAFRI  FORTRAN  Al*\  WILL  ALLOW  THE  USER  TO  PRO- 
C  DUCE  FOUR  GRAPHICAL  PLOTS  WHICH  WILL  CLOSELY  APPROXIMATE 
C  THE  PERFORMANCE  OF  AN  ACTUAL  THERMOACOUSTIC  HEAT  ENGINE. 

C  THE  GOAL  OF  THE  PROGRAM/SIMULATION  IS  TO  ALLOW  THE 
C  OPERATOR  TO  VARY  DIFFERENT  PARAMETERS  TO  ACHIEVE  WHAT  IS 
C  BELIEVED  TO  BE  AN  OPTIMUM  CONDITION  WHERE  THE  ENGINES 
C  DRIVER  RESONANCE  IS  MATCHED  WITH  THE  TUBE  RESONANCE. 

C 

C  GENERAL  INSTRUCTIONS: 

C 

C  THIS  PROGRAM  IS  SET  UP  TO  PROVIDE  FOUR  PLOTS:  VOLTAGE 
C  VERSUS  FREQUENCY,  PISTON  VELOCITY  VERSUS  FREQUENCY.  ACOUS- 
C  TIC  POWER  VERSUS  FREQUENCY,  AND  ELECTROACOUSTIC  EFFICIENCY 
C  VERSUS  FREQUENCY.  IN  GENERAL  THE  PROGRAM  HAS  HARD  INPUT 
C  DATA  SUCH  THAT  ALL  YOU  MUST  DO  TO  GET  A  SET  OF  PLOTS  IS 
C  COMPILE  "GRAFRI  FORTRAN  Al"  IN  FORTVS  AND  THEN  RUN  "FRIDGE 
C  FORTRAN  Al"  IN  WF77  FROM  FLIST  ON  THE  NPS  MAINFRAME  USING 
C  A  TEKTRONICS  618  TERMINAL.  THIS  ACTION  WILL  GIVE  YOU  A  SET 

C  OF  PLOTS  FOR  THE  INPUT  DATA.  TO  TRY  OTHER  DATA  YOU  MUST  GO 

C  TO  XEDIT  ON  THE  "FRIDGE  FORTRAN  Al"  PROGRAM  AND  VARY  THE 

C  INPUT  DATA.  AS  LONG  AS  YOUR  CHANGES  ARE  RESTRICTED  TO  THE 
C  "FRIDGE"  PROGRAM  THE  ONLY  ACTION  NECESSARY  TO  PRODUCE  NEW 
C  PLOTS  IS  TO  RUN  THE  PROGRAM  IN  WF77.  IF  YOU  FIND  IT  NECES- 
C  SARY  TO  CHANGE  THE  PLOTTING  PROGRAM,  "GRAFRI",  FOR  WHAT- 
C  EVER  REASON.  IT  MUST  BE  RECOMPILED  IN  FORTVS  PRIOR  TO  THE 

C  NEXT  ATTEMPT  TO  PLOT  OR  THE  CHANGES  WILL  NOT  BE  SEEN  BY 

C  THE  MACHINE. 

C 

C  INPUT  SECTION: 

C  A.  DEFINITION  OF  VARIABLES: 

C  ALL  OF  THE  VARIABLES  HAVE  BEEN  ESTABLISHED  AS  REAL 

C  OR  COMPLEX  QUANTITIES  BASED  ON  THE  REQUIREMENTS  OF  THE 
C  ARITHMETIC  ALGOR I THOM  OF  THE  CONCEPTUAL  MODEL.  EACH  OF 
C  THE  VARIABLES  WILL  NOW  BE  DEFINED. 

C  LCOIL-THE  IMAGINARY  PART  OF  THE  IMPEDANCE  OF  THE  VOICE 
C  COIL  OF  THE  DRIVER.  AND  IS  CALCULATED  TO  BE  W*LE. 

C  K«THE  PROPAGATION  CONSTANT  OF  THE  WAVES  IN  THE  TUBE 
C  AND  IS  GIVEN  BY  W/C . 
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C  IMZM*THE  IMAGINARY  PART  OF  THE  MECHANICAL  IMPEDANCE  OF 
C  THE  TUBE/DRIVER  COMBINATION.  XT  IS  GIVEN  BY  A  RATHER  LONG 
C  AND  COMPLEX  EQUATION  IN  THE  CALCULATION  SECTION  OF  THE 

0  program . 

C  MAGU-THE  MAGNITUDE  OF  THE  COMPLEX  QUANTITY  "U” ,  IN 

C  METERS/SEC. 

C  IMZT-THE  IMAGINARY  PART  OF  THE  TUBE  MECHANICAL  IMPE- 

C  DANCE  AND  IS  GIVEN  BY  SUBTRACTING  THE  DRIVER  PORTIONS  FROM 
C  XMZM. 

C  M-THE  MOVING  MASS  OF  THE  SYSTEM,  AND  IS  THE  SUM  OF  THE 
C  VOICE  COIL  AND  THE  DIAPHRAGM  MASSES  OF  THE  DRIVER,  IN 
C  KILOGRAMS. 

C  L-THE  LENGTH  OF  THE  TUBE  IN  METERS. 

C  LE«THE  INDUCTANCE  OF  THE  VOICE  COIL,  IN  HENRYS. 

C  I*THS  INPUT  CURRENT  FROM  THE  SOURCE,  IN  AMPS. 

C  MAGV-THE  MAGNITUDE  OF  THE  COMPLEX  QUANTITY  "V",  IN 

C  VOLTS. 

C  ZM-THE  MECHANICAL  IMPEDANCE  OF  THE  TUBE/DRIVER  SYSTEM, 
C  AND  IT  IS  MADE  UP  OF  THE  REAL  PART.  REZM,  AND  THE  'IMAGI- 
C  NARY  PART,  IMZM. 

C  ZE-THE  ELECTRICAL  EQUIVALENT  OF  THE  MECHANICAL  IMPE- 

C  DANCE  OF  THE  TUBE/DRIVER  SYSTEM.  IT  IS  GIVEN  BY  TRANSFOR- 
C  MING  ZM  ACROSS  A  TRANSFORMER  USING  THE  BL  COEFFICIENT  OF 
C  THE  DRIVER. 

C  ZCOIL-THE  COMPLEX  ELECTRICAL  IMPEDANCE  OF  THE  VOICE 

C  COIL. 

C  ZTOT-THE  TOTAL  ELECTRICAL  IMPEDANCE  OF  THE  SYSTEM. 

C  COMPLEX. 

C  V-THE  COMPLEX  VOLTAGE  OF  THE  SOURCE  NECESSSARY  TO 

C  PRODUCE  THE  SELECTED  INPUT  CURRENT. 

C  U-THE  VELOCITY  OF  THE  PISTON, (DIAPHRAGM) ,  OF  THE 

C  DRIVER  IN  M/S. 

C  ZTUBE*THE  COMPLEX  MECHANICAL  IMPEDANCE  OF  THE  TUBE 

C  ITSELF. 

C  P»THE  COMPLEX  PRESSURE  OF  THE  SYSTEM,  IN  NEWTONS /M* *2 . 
C  PI«A  CONSTANT  FOR  CONVERSIONS  TO  AND  FROM  ANGULAR 

C  UNITS. 

C  RM=THE  MECHANICAL  RESISTANCE  OF  THE  MOVING  MASS  OF  THE 
C  SYSTEM. 

C  SSUP-THE  MECHANICAL  STIFFNESS  OF  THE  SYSTEM  DUE  TO  THE 
C  DRIVER  SUSPENSION. 

C  RHO-THE  DENSITY  OF  THE  GAS  IN  THE  TUBE  AT  A  SPECIFIED 
C  MEAN  TEMPERATURE. 

C  C»THE  SPEED  OF  SOUND  FOR  THE  SPECIFIED  GAS  AND 

C  TEMPERATURE. 

C  S-THE  CROSS  SECTIONAL  AREA  OF  THE  TUBE. 

C  BL*A  COEFFICIENT  OF  THE  DRIVER  THAT  IS  BASED  ON  THE 

C  STEADY  AIR  GAP  FLUX  DENSITY  OF  THE  DRIVER  AND  THE  LENGTH 
C  OF  ITS  VOICE  COIL. 

C  Q-THE  QUALITY  FACTOR  OF  THE  TUBE. 

C  FO-THE  RESONANCE  FREQUENCY  OF  THE  TUBE. 

C  A-THE  CROSS  SECTIONAL  AREA  OF  THE  BACK  VOLUME  OF  THE 
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C  DRIVER , 

C  VOL* THE  VOLUME  OF  THE  BACK  VOLUME  OF  THE  DRIVER. 

C  RO-THX  INTERNAL  ELECTRICAL  RESISTANCE  OF  THE  POWER 

C  SOURCE. 

C  RE "THE  ELECTRICAL  RESISTANCE  OF  THE  VOICE  COIL. 

C  SVOL-THB  MECHANICAL  STIFFNESS  OF  THE  SYSTEM  DUE  THE 
C  BACK  VOLUME  OF  THE  DRIVER. 

C  WO* THE  RESONANCE  FREQUENCY  OF  THE  TUBE  IN  ANGULAR 

C  UNITS,  RAD/SEC. 

C  ALPHA* A  CHARACTERISTIC  QUANTITY  OF  THE  TUBE  THAT  IS  A 

C  FUNCTION  OF  THE  RESONANT  FREQUENCY,  THE  Q  AND  THE  SPEED  OF 
C  SOUND. 

C  F* FREQUENCY,  IN  HERTZ. 

C  W-THE  FREQUENCY  CONVERTED  TO  RAD/SXC. 

C  FORCE-THE  FORCE  ON  THE  MOVING  MASS  DUE  TO  THE  INTER- 
C  ACTION  OF  THE  VOICE  COIL  MAGNETIC  FIELD  AND  THE  DRIVER 
C  MAGNET. 

C  REZT-THE  REAL  PA  1.1  OF  THE  COMPLEX  MECHANICAL  TUBE 

C  IMPEDANCE. 

C  PU-THE  PRESSURE  TIMES  THE  VOLUME  VELOCITY.  ACOUSTIC 

C  POWER. 

C  PWR-THE  ELECTRICAL  POWER  NECESSARY  TO  DRIVE  THE  SYSTEM 
C  AT  THE  SPECIFIED  INPUT  CURRENT. 

C  PUPWR-A  POWER  RATIO,  ACOUSTIC  TO  ELECTRICAL.  AN  EFFI- 
C  CIENCY. 

C 

C  B . DATA  DESCRIPTION. 

C  FOR  THE  SELECTED  INPUT  DATA  THAT  IS  HARD  INPUT  THE 

C  FOLLOWING  VALUES  WERE  EMPIRICALLY  DETERMINED:  M,  L,  Q,  FO, 
C  A,  VOL,  BL.  S,  SSUP,  AND  RE.  RHO  AND  C  WERE  SELECTED  FOR 
C  HELIUM  GAS  AT  10  ATMOSPHERES  AND  A  MEAN  TEMPERATURE  OF  300 
C  DEGREES  KELVIN.  RM  WAS  CALCULATED  FROM  A  THEORETICAL 
C  EQUATION  USING  RHO.C  S  AND  THE  DRIVER  RESONANT  FREQUENCY. 
C  RG  WAS  SET  TO  ZERO  TO  MODEL  AN  IDEAL  POWER  SOURCE.  LE  WAS 
C  SET  TO  ZERO  FOR  EASE  OF  CALCULATION,  BUT  WAS  TESTED  DURING 
C  RUNNING  TO  VERIFY  THAT  IT  HAD  LITTLE  IF  ANY  EFFECT  WHEN 
C  INCLUDED.  A  CURRENT,  I.  OF  ONE  AMPERE  WAS  SELECTED  FOR 
C  EASE  OF  ANALYSIS. 

C  SINCE  YOU  MAY  WANT  TO  CHANGE  MY  INPUT  DATA  I  WILL  PUT 
C  A  COPY  IN  THE  COMMENTS  SO  THAT  YOU  CAN  RETURN  TO  IT  IF 
C  DESIRED.  MY  ANALYSIS  HAS  SHOWN  THE  BEST  BACK  VOLUME  TO  BE 
C  9.50E-8,  AND  THE  BEST  RM  TO  BE  1.J35  KG/S.  THIS  VALUE  FOR 
C  RM  WAS  DETERMINED  AFTER  REVIEW  OF  PLOTS  SHOWED  HIGHER  THAN 
C  EXPECTED  VOLTAGES  AND  IT  WAS  DETERMINED  BY  MEASUREMENT  IN 
C  THE  LABORATORY. 

C 

C 

C  DATA  RM, M, SSUP, RHO, C/3. 53E-3, 0.0045, 2. 6E4, 1.68, 1000/ 

C  DATA  S, BL ,Q, FO, A/1 . 521 E-3 , 8 , 15,800 , 7 . 13E-5/ 

C  DATA  VOL , RG , LS , RE , I , L/2 . 715E-6 , 0 , 0 , 8 , 1 , 0 . 625/ 
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REAL  LCOIL . K , IMZM , MAOU , ZMZT , M , L , LE , I , MAGV 
COMPLEX  ZM , ZE , ZCOIL , ZTOT , V , 0 , ZTUBE , P 
PI«3 . 141S926S4 

DATA  RM,M,SSUP,RHO, C/3. 531-3,0.0045, 2. 6E4, 1.68,1  "0/ 
DATA  S ,  BL ,  Q ,  PO , A/1 . 521E-3 ,8,15, 800 , 7 . 13K-5/ 

DATA  VOL , RG , LE , RE , I , L/2 . 71SE-6 , 0 , 0 , 8 , 1 , 0 . 625/ 

SVOL- ( A* A*RHO*C*C ) /VOL 
WO»2*PI*FO 
ALPHA-MO/ (2*Q*C) 


C  CALCULATION  SECTION: 

C  THE  NEXT  SIX  LINES  CREATE  AN  OUTPUT  FILE  FOR  THE  DATA 
C  THAT  IS  NECESSARY  TO  PRODUCE  THE  DESIRED  PLOTS,  AND  PLACES 
C  HEADINGS  ON  THE  DATA  AS  APPROPRIATE. 


OPEN (UNIT«69 , FILE*  *  FRIDGE  OUTPUT  A1 * , IOSTAT-IOS) 

WRITE (69, * ) '  MASS  VOLUME* 

WRITE ( 69 , 1 ) M , VOL 

WRITE ( 69 , * ) *  FREQUENCY  PIST  VEL  PRESS*PIST  VEL  PWR  RATIO 
*  VOLTAGE’ 

1  FORMAT (K13. 4, El 3. 4) 


C  THE  FOLLOWING  DO  LOOP  PROVIDES  FOR  THE  CALCULATION  OF 
C  THREE  HUNDRED  DATA  POINTS  FOR  EACH  OF  THE  DESIRED  PLOTS, 
C  USING  A  FREQUENCY  RANGE  OF  ONE  QUARTER  TUBE  RESONANCE  TO 
C  THREE  HALVES  TUBE  RESONANCE. 


DO  10  J»l,300 

F= ( FO/4 )  +  <  <  J- 1 ) *  < FO/240 ) ) 


C  THE  NEXT  THREE  LINES  CALCULATE  VARIOUS  QUANTITIES, 
C  EXPLAINED  ABOVE,  THAT  ARE  FUNCTIONS  OF  FREQUENCY,  FOR  USE 
C  IN  FUTURE  CALCULATIONS. 


W»2*PI«F 

K-W/C 

LCOIL=W«LE 


C  THE  NEXT  NINE  LINES  CALCULATE  THE  REAL  AND  IMAGINARY 

C  PARTS  OF  THE  TUBE/DRIVER  MECHANICAL  IMPEDANCE  AND  COMBINE 
C  THEM  INTO  A  COMPLEX  QUANTITY,  ZM. 


REZM=RM+ ( ( (RHO*C*S* ALPHA *COS(K*L) *SIN(K*L) ) + (K*RHO*C*S 
*  *SINH( ALPHA*L) *COSH( ALPHA* L) ) ) / (K* ( 1+ < ALPHA/K) * *2) * ( ( ( ( 
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*  SIN(K*L) ) **2) *( (COSH ( ALPHA* L) ) **2) >♦( < (COS(K«L) ) «*2) * ( ( 

*  SXNH(ALPHA*L) ) **2) ) ) ) > 

IMZMa  <  W  «M )-( SSUP/M )-(  SVOL/M  )♦(((  RHO  *C  *  S«  ALPHA  *  SINH  ( ALPHA  •  L ) 

*  *COSH( ALPHA *L) ) -(K*RHO*C*S*COS(K*L) *SIN(K*L) ) )/(K* ( 1+ (ALPHA 

*  /K) **2) *( ( ( (SIN(K*L) ) **2)  *( (CO$H( ALPHA* L) ) **2) )♦( ( (COS(K*L) ) * 

*  **2) * ( (SINH( ALPHA* L) )**2))))) 

ZM-CMPLX ( RIZM , IMZM ) 


C  THE  NEXT  THREE  LINES  CONVERT  ZM  TO  AN  ELECTRICAL 
C  QUANTITY  AND  THEN  COMBINE  ALL  OF  THE  ELECTRICAL  IMPEDANCES 
C  TO  GET  A  TOTAL  IMPEDANCE. 


ZE«(BL**2)/ZM 
ZCOIL-CMPLX ( RE , LCOIL ) 
ZTOT-RO+ZCOIL+ZS 


C  THE  NEXT  SIXTEEN  LINES  CALCULATE  THE  VALUES  TO  BE 
C  PLOTTED,  AND  PUT  THEM  IN  THE  OUTPUT  FILE,  AND  THEN  CLOSE 
C  THE  FILE. 


V«I*ZTOT 

FORCE-BL*! 

U-FORCE/ZM 

MAGU-CABS(U) 

REZT-REZM-RM 

IMZT-IMZM- ( W*M )  ■►  ( SSUP/W) ♦ ( SVOL/W ) 
ZTUBE-CMPLX { REZT , IMZT ) 

P»(ZTUBE*U)/S 

PU«0.5*(REAL(P*CONJG(U) )*S) 
PWR»0.5*(REAL(I*CONJG(V) ) ) 
PUPWR-PU/PWR 
MAGV-CABS(V) 

WRITE ( 69 , 100 ) F , MAGU , PU , PUPWR , MAGV 
100  FORMAT (E13 . 5.E13. S , E13 . 5 , El3. 5, E13 .5) 

10  CONTINUE 
CLOSE (69) 


C  PLOTTING  SECTION: 

C  THE  NEXT  TWO  LINES  INTERFACE  WITH  MY  PLOTTING  PROGRAM, 
C  "GRAFRI  FROTRAN  Al'\  AND  THE  DISSPLA  ROUTINE  TO  PLOT  THE 
C  DATA. 


PAUSE  ’PRESS  ENTER  TO  START  PLOTTING  ROUTINE’ 
CALL  EXCMS ( ’ EXEC  DISSPLA  GRAFRI’) 

STOP 

END 
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APPENDIX  B.  COMPUTER  MODEL  GRAPH  PROGRAM 


PROGRAM  GRAFRI 

C  THIS  IS  LT.  MICHAEL  P.  SUSALLA’S  PLOTTING  PROGRAM  FOR 
C  USE  IN  CONJUNCTION  WITH  "FRIDGE  FORTRAN  Al"  IN  “THE  THER- 
C  MOACOUSTIC  HEAT  ENGINE  SIMULALTION" .  THIS  PROGRAM  MUST  BE 
C  COMPILED  IN  FORTVS  PRIOR  TO  RUNNING  "FRIDGE  FORTRAN  Al"  IN 
C  WF77 .  THIS  PROGRAM  IS  FOR  USE  IN  DISSPLA  AND  CALLS  VARIOUS 
C  DISSPLA  SUBROUTINES.  THESE  VARIOUS  SUBROUTINES  ARE  USED 
C  MORE  THAN  ONCE  IN  THE  PROGRAM  BUT  WILL  ONLY  BE  DESCRIBED 
C  ONCE  EACH.  IN  ORDER  TO  USE  DISSPLA  YOU  MUST  HAVE  ARRAYS 
C  OF  DATA  POINTS  TO  BE  PLOTTED.  MY  ARRAYS  ARE  ALL  DIMEN- 

C  SIONED  TO  300  POINTS  SINCE  THAT  IS  THE  NUMBER  CALCULATED 
C  BY  "FRIDGE  FORTRAN  Al".  YOU  SHOULD  BE  AWARE  THAT  IF  YOU 
C  CHANGE  THE  NUMBER  OF  POINTS  TO  BE  CALCULATED,  AND  THE 
C  NUMBER  IS  GREATER  THAN  300,  YOU  MUST  REDIMENSION  THE 
C  ARRAYS.  TWO  CHARACTER  STRINGS  ARE  USED  IN  THIS  PROGRAM. 
C  NEW IS  IS  USED  IN  RENAMING  FILES  TO  THE  FORMAT  NECESSARY 

C  FOR  DISSPLA.  AND  DUMMYS  IS  USED  TO  READ  THE  HEADINGS  OUT 

C  OF  THE  DATA  FILE.  ALL  OF  THE  DATA  TO  BE  PLOTTED  IS  IN  THE 
C  FORM  OF  REAL  NUMBERS,  SO  ALL  OF  THE  VARIABLES  ARE  MADE 
*  C  REAL.  D  IS  USED  FOR  THE  FREQUENCY  DATA.  G  IS  USED  FOR  THE 

C  PISTON  VELOCITY  DATA,  H  IS  USED  FOR  THE  ACOUSTIC  POWER 
C  DATA,  O  IS  USED  FOR  THE  ELECTROACOUSTIC  EFFICIENCY  DATA 
C  AND.  T  IS  USED  FOR  THE  VOLTAGE  DATA.  THE  XMAX  AND  YMAX 
C  VARIABLES  ARE  INITIALIZED  TO  ZERO  AND  ARE  USED  TO  CALCU- 
C  LATE  THE  MAX  VALUE  OF  EACH  OF  THE  PARAMETERS  AS  THEY  ARE 

C  READ  INTO  THE  ARRAYS.  THIS  COMPLETES  THE  INTRODUCTION  AND 

C  EXPLAINS  THE  FIRST  FIVE  LINES  OF  THE  PROGRAM. 


DIMENSION  D(300) ,G(300) ,H(300) .0(300) ,T(300) 
CHARACTER *80  NEW1$, DUMMYS 
REAL  D,G,H,0,T,M 

DATA  XMAX . YMAX1 , YMAX2 , YMAX3 , YMAX4/0 .,0.,0.,0.,0./ 


C  THE  NEXT  SEVEN  LINES  OF  THE  PROGRAM  RENAME  MY  DATA 
C  FILE  TO  THE  FORMAT  NECESSARY  FOR  USE  IN  DISSPLA,  OPEN  THE 
C  FILE,  READ  THE  MASS  AND  VOLUME  OUT  OF  THE  FILE  AND,  READ 
C  THE  HEADINGS  OUT  OF  THE  FILE. 


NEW1$*’ RENAME  '//'  FRIDGE  OUTPUT  Al'//'  FILE  XYZ1  Al * 
CALL  EXCMS ( NEW 1 $ ) 

OPEN  ( 69 , FILEs  *  XYZ1 ’ ) 

READ  (69,’<A80)’)  DUMMYS 
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MAD  (69, 1)M, VOL 
READ  (69,’(A80)’)  DUMMY* 
1  FORMAT  (113.4,113.4) 


C  TH1  FOLLOWING  DO  LOOP  READS  THE  PLOTTING  DATA  INTO  THE 
C  ABOVE  SPECIFIED  ARRAYS  FOR  USE  IN  THE  DISSPLA  SUBROUTINES, 

C  AND  ALSO  CALCULATES  THE  MAXIMUM  VALUE  OF  EACH  OF  THE 
C  ARRAYS  FOR  SETTING  UP  THE  SCALES  OF  THE  PLOT  AXES.  THEN 
C  THE  DATA  FILE  IS  CLOSED.  THESE  VALUES  AM  LEFT  AS  REAL 
C  NUMBERS  VICE  INTEGERS  BECAUSE  THE  DATA  FROM  THE  PROGRAM 
C  CAN  VARY  SIGNIFICANTLY  WITH  DIFFERENT  INPUT  VALUES,  AND  IF 
C  THE  MAX  VALUE  IS  LESS  THAN  10  AN  INTEGER  REQUIREMENT  GIVES 
C  A  PLOTTING  ERROR. 


DO  20  J*1 ,300 

READ  ( 69 , IS , END*21 ) D ( J ) ,0(J> ,H<  J) ,0(J> ,T( J) 
IS  FORMAT ( 113 . 5 , 113 . S , E13 . 5 , 113 . 5 , 113 . S) 

IF(D(  J)  . GT . XMAX ) THEN 
XMAX-D(J) 

ENDIF 

IF(G( J) , GT . YMAX1 ) THEN 
YMAXl-Q(J) 

ENDIF 

IF(H( J) . GT . YMAX2 ) THEN 
YMAX2-H(J) 

ENDIF 

IF(0( J) . GT . YMAX3 ) THEN 
YKAX3-0(J) 

ENDIF 

IF(T( J) . GT . YMAX4 ) THEN 
YMAX4-T(J) 

ENDIF 

20  CONTINUE 

21  CL0SK69) 


C  THE  NEXT  THREE  LINES  TARE  CARE  OF  ANOTHER  ODDITY  OF 
C  THE  DATA.  IN  THE  CASE  OF  THE  VOLTAGE  PLOT.  DEPENDING  ON 
C  THE  INPUT  DATA  THE  VOLTAGE  WILL  LOCK  IN  AT  8  VOLTS,  AS 
C  EXPECTED.  BUT  THIS  IS  ALSO  THE  MAX  VALUE.  SO  YOU  GET  A 
C  STRAIGHT  LINE  AT  THE  TOP  OF  THE  PLOT  MARGIN.  THE  NEXT 
C  THREE  LINES  INCREASE  THE  MAX  VALUE  TO  9  VOLTS  TO  MAKE  THIS 
C  PLOTTING  CASE  MOM  CLEAR. 


IF ( YMAX4 . LE .9.0) THEN 

YMAX4-9.0 

ENDIF 
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C  THE  NEXT  EIGHT  LINES  CALCULATE  THE  INCREMENTAL  STEPS 
C  FOR  EACH  AXIS  OF  THE  VARIOUS  PLOTS.  THE  X  AXIS  IS  THE 
C  SAME  FOR  ALL,  AND  IS  LEFT  IN  INTEGER  FORMAT  FOR  CLARITY. 
C  EACH  AXIS  IS  SET  UP  FOR  TEN  INCREMENTS.  THE  DATA  FILE  IS 
C  THEN  RENAMED  TO  THE  DESIRED  FORMAT  FOR  MY  STORAGE. 


INC*INT(XMAX/10. ) 

YINC1-YMAX1/10. 

YINC2-YMAX2/10. 

YINC3-YMAX3/10 . 

YINC4*YMAX4/10. 

MAX=INT (XMAX ) 

NEW1$= *  RENAME  '//’  FILE  XYZ1  Al’//'  FRIDGE  OUTPUT  A1 ’ 
CALL  EXCMS(NEWIS) 


C  THIS  IS,  THE  FIRST  PLOTTING  SECTION.  IT  DOES  THE 
C  INITIAL  SET  UP  OF  ALL  OF  THE  CALLS  THAT  CAN  BE  COMMON  AND 
C  PLOTS  THE  VOLTAGE  VERSUS  FREQUENCY  PLOT.  CALL  TEK618  SETS 
C  UP  THE  REQUIRED  MACHINE  FOR  THE  PLOTTING  AND  LINKS  YOU  TO 
C  THE  PLOTTER.  CALL  BLOWUP  IS  USED  TO  MAGNIFY  THE  PLOT  FOR 
C  A  BETTER  PICTURE.  CALL  YAXANG  PUTS  THE  AXIS  NUMBERS 
C  PERPINDICULAR  TO  THE  PLOT  AXIS  TO  PREVENT  RUNNING  INTO 
C  EACH  OTHER.  CALL  PAGE  SETS  UP  THE  PAGE  SIZE  IN  INCHES. 
C  CALL  AREA2D  SETS  UP  A  TWO  DIMENSIONAL  PLOT  AREA  OF  THE 
C  SPECIFIED  SIZE  IN  INCHES.  CALL  SCMPLX  SELECTS  THE  LET- 
C  TERING  TYPE  FOR  THE  HEADING.  CALL  MESSAG  AND  CALL  REALNO 
C  ALLOW  ME  TO  PLACE  A  MESSAGE  IN  THE  PLOT  AREA  AT  A  SPECIFIC 
C  POINT  SPECIFIED  IN  INCHES  IN  THE  CALL.  CALL  XNAME  AND 
C  YNAME  ALLOW  ME  TO  LABEL  THE  AXES.  CALL  HEADIN  ALLOWS  ME 
C  TO  PLACE  A  TITLE  ON  THE  PLOT.  CALL  THKFRM  AND  FRAME  ALLOW 
C  ME  TO  PUT  FRAMES  AROUND  THE  PLOT  AREA  AND  THE  WHOLE  PLOT 
C  OF  A  SET  THICKNESS.  CALL  GRAF  SETS  UP  THE  PHYSICAL  ORIGINS 
C  AND  MAX  VALUES  AND  THE  INCREMENTS  FOR  EACH  AXIS.  CALL 
C  CURVE  TELLS  THE  MACHINE  TO  PLOT  THE  DESIRED  ARRAYS  OF  DATA 
C  POINTS.  CALL  ENDPL  ENDS  THAT  SPECIFIC  PLOT  AND  GIVES  A 
C  NEW  PAGE  FOR  THE  NEXT  PLOT. 


CALL  TEK618 
CALL  BLOWUP (1.5) 

CALL  YAXANG ( 0 ) 

CALL  PAGE (11,8.5) 

CALL  AREA2D ( 8 . ,6. ) 

CALL  SCMPLX 

CALL  MESSAG ( ' MASS=$  * ,100,0.1,5.5) 

CALL  REALNO (M, -2 , ’ ABUT ' , ’ ABUT ’ ) 

CALL  MESSAG ( * KG$ ' , 100 ,  ’ ABUT  * , ’ ABUT ’ ) 
CALL  MESSAG ( ’VOLUME=$’ ,100,0.1,5. ) 

CALL  REALNO ( VOL , -2 , ' ABUT ' , ’ ABUT ' ) 

CALL  MESSAG ( ’M**3$' ,100, ’ABUT' , ' ABUT ’ ) 
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CALL  XNAMfi ( ’ FREQUENCY ( IN  HERTZ ) $  \  1 00 ) 

CALL  YNAMt  ( '  VOLTAGE  { IN  VOLTS )  $  MOO ) 

CALL  HEADIN ( ' VOLTAGE  VS .  FREQUENCY* ' ,100,2. ,1) 
CALL  THKFRM < .02) 

CALL  FRAME 

CALL  GRAF (0.0, INC , MAX ,7.5, YINC4 , YMAX4 ) 

CALL  CURVE(D,T, 300,0) 

CALL  ENDPL ( 0 ) 


C  THIS  fECTION  PLOTS  PISTON  VELOCITY  VERSUS  FREQUENCY. 


CALL  AREA2D<8. ,6.) 

CALL  MESSAG( 'MASSES' ,100,0.1,5.5) 

CALL  RSALNO ( M , -2 , * ABUT * , ' ABUT • ) 

CALL  MESSAG ( ’ KGS ' ,100, 1  ABUT ' , • ABUT  * ) 

CALL  MESSAG ( 'VOLUME-* ',100,0.1,5.) 

CALL  REALNO ( VOL , -2 , ' ABUT ' , ' ABUT ' ) 

CALL  MESSAG ( 'M**3*' ,100, 'ABUT' , 'ABUT') 

CALL  XNAME( 'FREQUENCY (IN  HERTZ)* *, 100) 

CALL  YNAME( 'PISTON  VELOCITYdN  METERS/SEC)* ’,100) 

CALL  HEADIN { 'PISTON  VELOCITY  VS.  FREQUENCY*’ , 100,2. , 1) 
CALL  THKFRM ( . 02 ) 

CALL  FRAME 

CALL  GRAF (0.0, INC , MAX ,0.0, YINC1 , YMAX1 ) 

CALL  CURVE ( D , G , 300 , 0 ) 

CALL  ENDPL ( 0 ) 


C  THIS  SECTION  PLOTS  ACOUSTIC  POWER  VERSUS  FREQUENCY. 


CALL  AREA2D(8. ,6. ) 

CALL  MESSAG( 'MASS-*' ,100,0.1,5.5) 

CALL  REALNO (M, -2, 'ABUT' , 'ABUT* ) 

CALL  MESSAG ( * KG* ' , 100 , ' ABUT ’ , ' ABUT ' ) 

CALL  MESSAG ( ' VOLUME-* ’,100,0.1,5.) 

CALL  REALNO ( VOL , -2 , ’ ABUT ' , ' ABUT ' ) 

CALL  MESSAG ( *M« *3* ' , 100 , ’ ABUT ' , 'ABUT' ) 

CALL  XNAME ( ' FREQUENCY ( IN  HERTZ ) $ ' , 1 00 ) 

CALL  YNAME( 'ACOUSTIC  POWER  (IN  N-M/S)*’ , 100) 

CALL  HEADIN ( ' ACOUSTIC  POWER  VS .  FREQUENCY* ' , 100 , 2 . , 1 ) 
CALL  THKFRM (.02) 

CALL  FRAME 

CALL  GRAF ( 0 . 0 , INC , MAX ,0.0, YINC2 , YMAX2 ) 

CALL  CURVE(D,H, 300,0) 

CALL  ENDPL ( 0 ) 


C  THIS  SECTION  PLOTS  POWER  RATIO  VERSUS  FREQUENCY. 
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CALL  AREA2D ( 8 . ,6. ) 

CALL  MESSAG( 'MASS**' ,100,0.1,5.5) 

CALL  REALNO(M,-2, 'ABUT' , 'ABUT' ) 

CALL  *«SSAG< 'K0$' ,100, 'ABUT' , 'ABUT' ) 

CALL  NSSSAG ( ' VOLUME** ’ , 100„ 0 . 1 , 5 . ) 

CALL  REALNO(VOL,-2, 'ABUT' , 'ABUT' ) 

CALL  MESSAG('M**3*' ,100, 'ABUT' , 'ABUT' ) 

CALL  XNAME ( ' FREQUENCY ( IN  HERTZ) $'  ,100) 

CALL  YftAMEt  'ACOUSTIC  PUR/  ELECTRIC  PURS ',100) 

CALL  HEADXN ( ' POUER  RATIO  VS .  FREQUENCYS ',100,2.,!) 

call  thkfrm(.02) 

CALL  FRAME 

CALL  GRAF (0.0, INC , MAX ,0.0, YINC3 , YMAX3 ) 

CALL  CURVE(D, 0,300,0) 

CALL  ENDPL(O) 


CALL  DONEPL  ENDS  THE  PLOTTING  SEQUENCE  AND  SENDS 
BACK  TO  CMS. 


CALL  DONEPL 

STOP 

END 


YOU 


DRIVER  HOUSING  -  RESONATOR  VIEW 


Figure  C 


(Du 


in  inches  ) 


2.  Driver  Housing — Resonator  View 


VIEW  AA 


ELECTRONICS  PLUG 

(Oleenaiona  in  inchaa  I 


Feed-Through  for  Electrical  Leads 
Shown  in  Figure  v-5. 


c 


* 


Figure  C-5.  Driver  Housin* — Side  View  Section  BB 
Showing  Pressure  Gage  and  Helium 
Fill  Ports. 
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Fill  parts . 


VIEW  AA 


DC  PRESSURE  TRANSDUCER 

(Olmnaiona  in  inches) 


Figure  C-6 .  DC  Pressure  Gage  Feed 
Through  Plug  Details. 


DRIVER  HOUSING  -  LID  VIEW 


Figure  C- 


Driver  Housing — Lid  View  Showing 
Capillary,  Microphone,  end  Follower, 
and  Bolt  Circlet  for  Resonator 
Attachment,  Pusher  Plate,  Pressure 
Lid  and  GAS  Can  Lid. 


ote :  Plot*  is  3/10'  wide  end 

end  slopes  10°  at  outer 
edge .  ) 


PUSHER  PLATE  (LID  VIEW) 


(Dimensions  in  inches) 


Fisure  C-8.  Aluainus  Pusher  Plate. 


(Not*:  Q-nno  Qroova  0.018*  wide 
and  . 036 ’  daap : 

ltd  la  1 /2 *  wide) 


PRESSURE  LID  (LID  VIEW) 

(Dimensions  in  inches) 


Figure  C-9.  Pressure  Lid  with  Lead 
O-ring  Groove. 


PARKER  2-174  O-RING 
fWoEFIH  bM  WIDTH 


Figure  C-10.  Pressure  and  Heliue  Leak  Test  Lid. 
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DRIVER  TEST  LID 


APPENDIX  D.  SIM  emis.  ULSI  2£  MATERIALS 


The  driver  consists  of  a  modified  JBL  model  2450J  neo¬ 
dymium-iron-boron  magnet  with  iron  pole  pieces  on  either 
side.  The  diaphragm  is  titanium  and  the  voice  coil  is 
aluminum. 

The  housing,  pusher  plate,  pressure  lid,  and  reducer 
cone  will  all  be  machined  from  6061-T6  aluminum.  All  open¬ 
ings  in  the  driver  housing  will  be  sealed  with  1/16”  or 
1/32"  diameter  lead  o-rings.  The  electrical  feed-throughs 
will  be  sealed  with  Emerson  &  Cuming’s  STYCAST<TW»  2S50FT 
epoxy.  All  wire  leads  are  copper.  The  steel  bolts  connecting 
the  resonator  to  the  driver  are  size  1/4-20.  All  other  steel 
bolts  are  size  10-32. 
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